0% found this document useful (0 votes)

26 views155 pages

Pre Solutions of Old Papers 1573360064915 1 1631425497745 1656220859798

MATLAB

Uploaded by

ChandraprakashAgrawal

We take content rights seriously. If you suspect this is your content, claim it here.

Available Formats

Download as PDF, TXT or read online on Scribd

0% found this document useful (0 votes)

26 views155 pages

Pre Solutions of Old Papers 1573360064915 1 1631425497745 1656220859798

MATLAB

Uploaded by

ChandraprakashAgrawal

We take content rights seriously. If you suspect this is your content, claim it here.

Available Formats

Download as PDF, TXT or read online on Scribd

You are on page 1/ 155

TA 164 Regular Mid Sem Test-Soultions.

docx
TA ZC 164 compre 17-18 Solutions.docx
TA ZC 164 EC-3 17-18 Solutions.docx
TA ZC 164 EC-3 17-18 Solutions.pdf
TA ZC164 EC-3R FIRST SEM 2018-2019-Revised solutions.docx
TA ZC164 EC-3R FIRST SEM 2018-2019-Revised Solutions.pdf
TA ZC164 EC-3R SECOND SEM 2017-2018-Solutions.docx
TA ZC164 EC-3R SECOND Semester 2018-19- Solutions.docx
TA ZC164 EC-II 2017-2018-Solutions.docx
CUMMINS-COMPRE SOLUTIONS.docx
MATLAB Old Paper+Solution.doc
Solutions-EC3M.docx
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2018-2019

Mid-Semester Test
(EC-2 Make-up)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 30% No. of Pages =2
Duration : 2 Hours No. of Questions =
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.

Q.1 Given x= [2, 7, 4, 8, 6, 7, 10]

What will be the result of y for the following relations? [3]

Solution

S.No Relation Result

(i) y=x>5 y=[ 0 1 0 1 1 1 1]
(ii) y=x(find(x>5)) y= 7 8 6 7 10
(iii) y=length(find(x>5)) y=5

𝟏. 𝟐 −𝟑. 𝟒 𝟑. 𝟔
Q.2 The array d has the values [𝟎. 𝟖 𝟏. 𝟒 −𝟎. 𝟖] [3]
𝟏. 𝟓 𝟎. 𝟔 𝟓. 𝟎
Solution
(i) What is the size of array d and which command will you d= 3 x 3 , size(d)
use in MATLAB to determine it
(ii) What is the value of d(3,2) ? 0.6
(iii) What subscripts does d (?,?)=0.6 will contain ? 3,2

Q.3. If m= [1,2,3,4;5,6,7,8;9,10,11,12] [3]

What elements will the following expression contains

(i) m(2:end,2:2:end) ?
ans =
6 8
10 12

(ii) How will the following command m(1:2,[1 2])=[20,21;22 23] ?

m=
20 21 3 4
22 23 7 8
9 10 11 12

The elements 1,2; 3,4 will be replace by 20,21;22,23

(iii) What is result of following m (2,:) ? ans = 22 23 7 8, It will print second row

Q.4 What will be the result of the following commands in MATLAB [3]
(i) a=eye(3) It will print an Identity Matrix of 3 *3
(ii) b=zeros(2,3) It will print a 2*3 matrix consisting of all zero elements
(iii) clear and clc clear will clear workspace variables and clc will clear the screen
Q.5 Write a MATLAB program that prints the areas of a cylindrical pressure vessel used in chemical industry. The radius varies from 4 cm to 4.5
cm with 0.1 cm increments. Determine the areas of such varying radii with a constant height of 5 m. The Area of the cylinder is 𝐴𝑟𝑒𝑎 = 2𝜋𝑟ℎ +
𝜋𝑟 2 [3]

h=5;
%area=zeros(1,5)
r=4.0:0.1:4.5
for i=1:5

area(i)=2*pi*r(i)*h+pi*r(i)^2

end
fprintf('%The area of cylinders with varying r is %f\n',area(i));

Q.6 Evaluate the following expression [3]

(i) 1<= 2 || ‘x’ = = ‘y’ Any condition is true result is true ?

ans =
logical
1

(ii) xor(3 >= 2,‘x’ = = ‘y’) ?

Any one condition is false result is 1
ans =
logical
1

(iii) (3 < 12) – 2 ?

( 3<12) condn is true means value=1 that is 1-2= -1
Q.7 Write if else statements for the following scenarios [4]
(a) Input the taxable income of a salaried employee
(b) If taxable income is less than Rs. 5,00,000 print zero tax. Tax rate=0%
(c) If taxable income is more than Rs. 5,00,000 and less than Rs.1000000. Tax rate=10% above 5,00000
(d) If taxable income is more than Rs. 10,00,000 and less than Rs.2000000. Tax rate=20% above 10,00000 and 10% above 500000-
1000000 range

% Consider x as Income

x=input ('Enter Income of the Employee');

if (x<500000)
tax1=0
elseif (x>500000 && x<1000000)
tax1=(x-500000)*0.10
elseif (x>1000000 && x<2000000)
tax1=(x-1000000)*0.20+500000*0.10
fprintf('The Tax for x is %.2f', x,tax1);
end
Q.8 The continuity equation in fluid dynamics for steady fluid flow through a stream tube equates the product of the density, velocity,
and area at two points that have varying cross-sectional areas. For incompressible flow, the densities are constant so the equation
𝐴1
is A1V1 = A2V2. If the areas and V1 are known, V2 can be found as 𝐴2 𝑉1.
Therefore, whether the velocity at the second point increases or decreases
depends on the areas at the two points. Write a script that will prompt the user
for the two areas in feet squared, and will print whether the velocity at the second
point will increase, decrease, or remain the same as at the first point [4]

% Continuity Equation Formula

A1=input(‘Enter area at Point 1 in ft2 \n’);

A2=input(‘Enter area at Point 2 in ft2 \n ’);
V1=input(‘Enter velocity in ft/s \n’);
V2= (A1/A2)*V1;
disp(V2)

% COMMENTS--- V2 will increase if A1 is greater that A2 that is V2 is more than V1 , Similarly viceversa, If Areas are same velocities
are same

Q.9. Water near the surface of a tropical ocean has a temperature of 298.2K (250oC), whereas water 700m beneath the surface has a
temperature of 280.2K (7.0oC). It has been proposed that the warm water be used as the hot reservoir and the cool water as the
cold reservoir of a heat engine. Find the maximum possible efficiency for such an engine. Analytical solution is give. Convert
the same into the necessary generalized function [4]

Solution:
Given data,
Th = 298.2 K
TI = 280.2 K
The maximum possible efficiency is the efficiency that a Carnot engine would have operating between these two temperatures. Substitute the value in the
corresponding formula,
η = (Th−TI)/Th
η = (298.2−280)/280)
η = 0.060 (6.0%)
Therefore, the maximum possible efficiency is 6 %

function eff=carnot(Th,Ti)
eff=(Th-Ti)/ Th * 100;
end
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Comprehensive Examination Test (EC-3 Regular)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book No. of Pages =3
Weightage : 40%
Duration : 3 Hours No. of Questions = 8
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer
book.
2. All parts of a question should be answered consecutively. Each answer should start from
a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

All Question Carries 5 Marks – 8*5=40 Marks

Q.1 Evaluate the following expressions in MATLAB

(a) (11/5)+6 (b) round(-11/5)+6
(b) (3^2)^3 (d) ceil(-11/5)+6
(c) 11/(5+6)

Solutions
(a) 8.2
(b) 4 Round to nearest decimal or integer
(c) 729
(d) 4 Round toward positive infinity
(e) 1

Q.2 Using the zeros, ones, and eye commands create the following arrays
0 0 1 1
0 0 1 1
a) [ ]
0 0 0 0
1 1 1 1
Solution (a)

Step 1 : Create a=zeros(4) will create a 4*4 matrix of 0

Step2 : a(:,3:4)=ones will replace contents of 3rd and 4th Column by ones
Step3: a(4,:)=ones will replace contents of 4th row with ones
Step4: a(3,:)=zeros will replace the 3rd row by zero

1 1 0 0 1
1 1 0 0 0
b) [ ]
1 1 0 0 0
1 1 0 0 0

Solution (b)

Step1 : z=ones(4) will create 4*4 matrix

Step2: z(:,3:5)=zeros will replace elements of 3 to 5th column with zeros
Step3: z(1,5)=ones will replace 5th column of 1st row by ones

Q.3 A juice company manufactures one-gallon bottles of three types of juice blends using
orange, pineapple, and mango juice. The blends have the following compositions:
1 gallon orange blend: 3 quarts of orange juice, 0.75 quart of pineapple juice,
0.25 quart of mango juice.
1 gallon pineapple blend: 1 quart of orange juice, 2.5 quarts of pineapple
juice, 0.5 quart of mango juice. 1 gallon mango blend: 0.5 quart of orange juice, 0.5
quart of pineapple juice,3 quarts of mango juice.

How many gallons of each blend can be manufactured if 7,600 gallons of orange juice,
4,900 gallons of pineapple juice, and 3,500 gallons mango juice are available? Write a
system of linear equations and solve in MATLAB only steps

Solution

Consider orange=x , pineapple=y, mango= z vectors

1 gallon Orange Blend: X

1 gallon Pineapple Blend: Y
1 gallon mango Blend: Z

Set up Equations
X=3*x+0.75*y+0.25*z
Y=x+2.5*y+0.5*z
Z=0.5*x+0.5*y+3*z

Now Put X=7600, Y=4900,Z=3500

MATLAB

a=[3,0.75,0.25;1,2.5,0.5;0.5,0.5,3]
b=[7600;4900;3500]

c=inv(a)* b

/* Answers not require only steps give full marks*/

Ans [ x= 2246 gallons orange juice, y=933 gallons pineapple, z=636 gallons mango juice]

Q.4 Viscosity, μ, is a property of gases and fluids that characterizes their resistance to flow.
For most materials viscosity is highly sensitive to temperature. Below is a table that
gives the viscosity of SAE 10W oil at different temperatures (Data from B.R. Munson,
D.F. Young, and T.H. Okiishi, Fundamentals of Fluid Mechanics, 4th ed., John Wiley
and Sons, 2002). Determine an equation that can be fitted to data.

T (0C) -20 0 20 40 60 80 100 120

µ (N 4 0.38 0.095 0.032 0.015 0.0078 0.0045 0.0032
s/m2)
x 10-5

The straight line does not fit the point , Therefore the equation can be fitted as
ln(𝜇) = 𝑎2 𝑇 2 + 𝑎1 𝑇+𝑎0 .

The function can be fitted to the data using polyfit(x,y,2) second degree
polynomial, T is independent variable in 0K (C+273) and ln µ. Write the
Program that best fits the program with calculation of coefficients and display
the data points and function.

Hint:Use semilog function

T=[-20:20:120];
mu=[4 0.38 0.095 0.032 0.015 0.0078 0.0045 0.0032];
TK=T+273;
p=polyfit(TK,log(mu),2)
Tplot=273+[-20:120];
muplot = exp(p(1)*Tplot.^2 + p(2)*Tplot + p(3));
semilogy(TK,mu,'o',Tplot,muplot)

Q.5 Create a data in some other program on your computer. (% use a separate file
say data) such as excel, word and so forth. Copy the data set to the clipboard
using the windows and then use function uiimport to load the data set into
MATLAB

Step1 : Create the Data in Excel file

Step2 : Say the file name as testexam.xlsx
Step3: The data in file is as follows

name Age marks

Samir 40 54
Ramesh 36 67
ravi 47 72

Step4: Opening in matlab

s=uiimport('C:\Users\BITS7\Documents\MATLAB\testexam.xlsx')

both test data and numeric data will be loaded

Q.6. Use ODE45 function to solve the differential equation

𝑡 2 𝑦 ′ = 𝑦 + 3𝑡
y(1)= -2

over the interval 1< t <4

>> ode

Q.7 Plot the function 𝑦 = 𝑒 −𝑥 𝑠𝑖𝑛𝑥 for x between 0 and 2 in steps of 0 1 . Create
the following plot types (a) stem plot (b) stair plot (c) bar plot . Be sure to
include titles and axis labels on all plots

x=0:0.1:2;
y=exp(-x)*sin(x);
stem(x,y);
bar(x,y);
title(‘stem and bar Plot’);
xaxis(‘x coordinate’)
yaxis(‘function’);

Q.8 Fluid flows in pipe networks can be analyzed in a manner similar to that used
for electric resistance networks. The Figure below shows a network with three
pipes. The volume flow rates in the pipes are q1, q2, and q3 . The pressures at the
pipe ends are pa, pb and pc, . The pressures at the junction is p1. Under certain
condition, the pressure flow rate relation in a pipe has the same form as the
voltage-current relation in a resistor.Thus, for the three pipes, we have

Where the Ri are the pipe resistances. From conservation of mass

q1= q2+ q3

a. Set up these equations in a matrix form Ax=b suitable for solving for the
three flow rates q1, q2, q3, and the pressure p1 given the values of pressures
pa, pb and pc and the values of resistances R1, R2, and R3. Find expressions
of A and B
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Comprehensive Examination Test (EC-3 Regular)

Course No. : TA ZC164

All Question Carries 5 Marks – 8*5=40 Marks

Q.1 Evaluate the following expressions in MATLAB

(a) (11/5)+6 (b) round(-11/5)+6
(b) (3^2)^3 (d) ceil(-11/5)+6
(c) 11/(5+6)

Solutions
(a) 8.2
(b) 4 Round to nearest decimal or integer
(c) 729
(d) 4 Round toward positive infinity
(e) 1

Q.2 Using the zeros, ones, and eye commands create the following arrays
0 0 1 1
0 0 1 1
a) [ ]
0 0 0 0
1 1 1 1
Solution (a)

Step 1 : Create a=zeros(4) will create a 4*4 matrix of 0

Step2 : a(:,3:4)=ones will replace contents of 3rd and 4th Column by ones
Step3: a(4,:)=ones will replace contents of 4th row with ones
Step4: a(3,:)=zeros will replace the 3rd row by zero

1 1 0 0 1
1 1 0 0 0
b) [ ]
1 1 0 0 0
1 1 0 0 0

Solution (b)

Step1 : z=ones(4) will create 4*4 matrix

Step2: z(:,3:5)=zeros will replace elements of 3 to 5th column with zeros
Step3: z(1,5)=ones will replace 5th column of 1st row by ones

Solution

Consider orange=x , pineapple=y, mango= z vectors

1 gallon Orange Blend: X

1 gallon Pineapple Blend: Y
1 gallon mango Blend: Z

Set up Equations
X=3*x+0.75*y+0.25*z
Y=x+2.5*y+0.5*z
Z=0.5*x+0.5*y+3*z

Now Put X=7600, Y=4900,Z=3500

MATLAB

a=[3,0.75,0.25;1,2.5,0.5;0.5,0.5,3]
b=[7600;4900;3500]

c=inv(a)* b

/* Answers not require only steps give full marks*/

Ans [ x= 2246 gallons orange juice, y=933 gallons pineapple, z=636 gallons mango juice]

T (0C) -20 0 20 40 60 80 100 120

µ (N 4 0.38 0.095 0.032 0.015 0.0078 0.0045 0.0032
s/m2)
x 10-5

The straight line does not fit the point , Therefore the equation can be fitted as
ln(𝜇) = 𝑎2 𝑇 2 + 𝑎1 𝑇+𝑎0 .

Hint:Use semilog function

Step1 : Create the Data in Excel file

Step2 : Say the file name as testexam.xlsx
Step3: The data in file is as follows

name Age marks

Samir 40 54
Ramesh 36 67
ravi 47 72

Step4: Opening in matlab

s=uiimport('C:\Users\BITS7\Documents\MATLAB\testexam.xlsx')

both test data and numeric data will be loaded

Q.6. Use ODE45 function to solve the differential equation

𝑡 2 𝑦 ′ = 𝑦 + 3𝑡
y(1)= -2

over the interval 1< t <4

>> ode

x=0:0.1:2;
y=exp(-x)*sin(x);
stem(x,y);
bar(x,y);
title(‘stem and bar Plot’);
xaxis(‘x coordinate’)
yaxis(‘function’);

Where the Ri are the pipe resistances. From conservation of mass

q1= q2+ q3

Course No. : TA ZC164

All Question Carries 5 Marks – 8*5=40 Marks

Q.1 Evaluate the following expressions in MATLAB

(a) (11/5)+6 (b) round(-11/5)+6
(b) (3^2)^3 (d) ceil(-11/5)+6
(c) 11/(5+6)

Solutions
(a) 8.2
(b) 4 Round to nearest decimal or integer
(c) 729
(d) 4 Round toward positive infinity
(e) 1

Q.2 Using the zeros, ones, and eye commands create the following arrays
0 0 1 1
0 0 1 1
a) [ ]
0 0 0 0
1 1 1 1
Solution (a)

Step 1 : Create a=zeros(4) will create a 4*4 matrix of 0

Step2 : a(:,3:4)=ones will replace contents of 3rd and 4th Column by ones
Step3: a(4,:)=ones will replace contents of 4th row with ones
Step4: a(3,:)=zeros will replace the 3rd row by zero

1 1 0 0 1
1 1 0 0 0
b) [ ]
1 1 0 0 0
1 1 0 0 0

Solution (b)

Step1 : z=ones(4) will create 4*4 matrix

Step2: z(:,3:5)=zeros will replace elements of 3 to 5th column with zeros
Step3: z(1,5)=ones will replace 5th column of 1st row by ones

Solution

Consider orange=x , pineapple=y, mango= z vectors

1 gallon Orange Blend: X

1 gallon Pineapple Blend: Y
1 gallon mango Blend: Z

Set up Equations
X=3*x+0.75*y+0.25*z
Y=x+2.5*y+0.5*z
Z=0.5*x+0.5*y+3*z

Now Put X=7600, Y=4900,Z=3500

MATLAB

a=[3,0.75,0.25;1,2.5,0.5;0.5,0.5,3]
b=[7600;4900;3500]

c=inv(a)* b

/* Answers not require only steps give full marks*/

Ans [ x= 2246 gallons orange juice, y=933 gallons pineapple, z=636 gallons mango juice]

T (0C) -20 0 20 40 60 80 100 120

µ (N 4 0.38 0.095 0.032 0.015 0.0078 0.0045 0.0032
s/m2)
x 10-5

The straight line does not fit the point , Therefore the equation can be fitted as
ln(𝜇) = 𝑎2 𝑇 2 + 𝑎1 𝑇+𝑎0 .

Hint:Use semilog function

Step1 : Create the Data in Excel file

Step2 : Say the file name as testexam.xlsx
Step3: The data in file is as follows

name Age marks

Samir 40 54
Ramesh 36 67
ravi 47 72

Step4: Opening in matlab

s=uiimport('C:\Users\BITS7\Documents\MATLAB\testexam.xlsx')

both test data and numeric data will be loaded

Q.6. Use ODE45 function to solve the differential equation

𝑡 2 𝑦 ′ = 𝑦 + 3𝑡
y(1)= -2

over the interval 1< t <4

>> ode

x=0:0.1:2;
y=exp(-x)*sin(x);
stem(x,y);
bar(x,y);
title(‘stem and bar Plot’);
xaxis(‘x coordinate’)
yaxis(‘function’);

Where the Ri are the pipe resistances. From conservation of mass

q1= q2+ q3

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 40%
No. of Pages =2
Duration : 3 Hours
No. of Questions = 7
Date of Exam : 25/11/2018 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1 Short Questions [2 x 10= 20]

a) Explain the terms in the following statement
fid=fopen (‘compre.dat’,’r’) {Terms: fid, fopen, compre, and r}

fid:File Identifier
fopen: File opening [1]
compre.dat: File name in data type
‘r’: opened for reading [1]

b) What is the use of fscanf and fprintf in file handling? Give 1 example each

fscanf( fileID , formatSpec ) reads data from an open text file into column vector A and interprets
values in the file according to the format specified by [1]

fprintf(fileID,formatSpec,A1,...,An) applies the formatSpec to all elements of arrays A1,...An in column

order, and writes the data to a text file. Fprintf uses the encoding scheme specified in the call to fopen
[1]

c) If a is [1,2;3,4;5,6] what will the command min(a) print in MATLAB, explain

min(a) will print the minimum values of the columns. It will print 1 and 2 in column
[1]

1 2 [1]

d) if a is [1,2;3,4;5,6] what will the command min(a,[ ],2) print in MATLAB,explain

min(a,[],2) will print the minimum values of the columns. It will print 1 and 2 in row
[1]

1
3 [1]
5

e) randi (2,3,[10 20]). How will this command work in MATLAB?

This command should be written as randi([10 20],2,3] . Anyways if
explanation is ok, give marks it will print a matrix of 2 by 3 with numbers
between 10 and 20 [1]

Example
a=
18 11 16
19 20 11 [1]

f) >> ppm = [ 9 8.5 7 6 10 5.5];

>> index = [3 5 1 6 4 1];
>> ppm(index) What will it print?

ans =
7.0000 10.0000 9.0000 5.5000 6.0000 9.0000

g) >>a=eye(3);
>>sparse(a) which elements will it print ?

Sparse will print non zero elements ------------------

It will print [1]

ans =
(1,1) 1
(2,2) 1
(3,3) 1 [1]

h) Construct an equivalent statement using single if statement in MATLAB for the following
block If x<y
If z<10
W=x*y*z

If (x<y & z<10) W=xyz; [2]

i) Accepts marks of the students as p and if p is greater than 50 print pass else print fail

p=input(‘Enter marks of the student’)

if (p>50)
disp(‘pass’)
else
disp(‘fail’)
end [2]

j) Which command will you use to find the solution of simultaneous equation and how?

if a is the vector of coefficient of simultaneous equation we will use inv(a) *

coefficient of right side to give solution of equation
X=inv(a)*Y; [2]
Q.2 Atmospheric properties such as temperature, air density, and air pressure are important
in aviation. Create a file that stores temperatures in degrees Kelvin at various altitudes.
The altitudes are in the first column and the temperatures in the second. For example, it
may look like this:

1000 288
2000 281
3000 269
5000 256
10000 223

Write a script that will load this data into a matrix, separate it into vectors, and then plot
the data with appropriate axis labels and a title. [4]

Solution

A=[1000,2000,3000,5000,10000];
T=[288, 281,269,256,223];
save atm.dat A T –ascii
load atm.dat
plot(A,T)
legend(‘Altitudes versus Temperature’)
xlabel(‘Altitude’)
ylabel(‘Temperature’)
TA ZC164 (EC-3 Regular) First Semester 2018-2019 Page 2

Q.3 For the polynomial f(x)=𝑥 4 − 7𝑥 3 + 59𝑥 2 − 75𝑥 + 35 [4]

Find the following using commands

(a) Roots of polynomial
(b) Calculate f(9)
Plot the polynomial for -1.5<=x<=6.7

Solution

>> P=[1 -7 59 -75 35]

>> roots(p) [1.5]

ans =
2.7947 + 6.5277i
2.7947 - 6.5277i (No marks , this is just output, only command is sufficient)
0.7053 + 0.4436i
0.7053 - 0.4436i

>> syms x
>> f=(@)x 𝑥 4 − 7𝑥 3 + 59𝑥 2 − 75𝑥 + 35
>> f(9) [1.5]

f(9)
ans = (No marks , this is just output, only command is sufficient)
5597

>> x=linspace(-1.5,6.7,20)
>> f=@(x) x.^4-7*x.^3+59*x.^2-75.*x+35 [1]
>> plot(x,f(x))

(No marks , this is just output

Q.4 Quenching is the process of immersing a hot metal object in a bath for a specified time
to obtain certain properties such as hardness. A copper sphere 25 mm in diameter,
initially at 300 0C, is immersed in a bath at 0 0C. The following table gives measurements
of the sphere’s temperature versus time. Find a functional description of these data. Plot
the function and the data on the same plot. [4]
Time 0 1 2 3 4 5 6
Temperature 300 150 75 35 12 5 2

% Script to find functional description of the data

time=[0 1 2 3 4 5 6];
temperature=[300 150 75 35 12 5 2];
polyfit(time,temperature,1)
%% It does not describe correct description as at time=0
temperature should be 300

polyfit(time,temperature,2)
%% It does not describe correct description as at time=0
temperature should be 300
y=polyfit(time,temperature,3)
plot(time,temperature)
% now checking
time=[0 1 2 3 4 5 6];
y= y(1)*time.^3+y(2)*time.^2+y(3)*time+y(4)

time=[0 1 2 3 4 5 6];
y= y(1)*time.^3+y(2)*time.^2+y(3)*time+y(4)
plot(time,temperature,'b',time,y,'r')

Now Both the same so functional description is right

For command 1 Marls

For Procedure 1 marks
For complete procedure and correct commands give full marks
Q.5 Create a Script PLOT to plot the following functions in MATLAB
x=cos(t) y=sint(t) and z=√𝑥 2 + 𝑦 2 over the range 0<=t<=10. Plot 3 dimensional
figures using plot3 and surf commands [4]

% Script to Plot 3 Dimensional Figure

time=linspace(0,10*pi,50)
x=cos(time);
y=sin(time)
z=sqrt(x.^2+y.^2)
figure(1)
plot3(x,y,z);
figure(2)
[X Y]=meshgrid(x,y);
surf(X Y z);
Give full marks only if steps are followed

Q.6 Write a function that will receive the temperature in deg oC and convert the temperature
to dego F and deg0 K. Use the standard conversion formula for deg C to deg F [4]

% function to convert deg C to degf

function degf =tempconv (degc)

degf=1.8*degc+32;
end
>> tempconv(100)

212 Give full marks only if steps are followed

*********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division

Assignment :22 Marks Last Date: 13 Dec 2018

Optimization Problems –Refer Tutorial 1

1. A firm manufactures two types of products P1 and P2 and sells them on a profit of Rs. 3 on type P1 and Rs. 4 on type P2. Each product is processed on two machines A
and B. Type P1 requires 2 minutes of processing time on A and one minute on B; type P2 requires 3 minutes on A and 2 minutes on B. The machine A is available for not
more than 7 hours 30 minutes while machine B is available for 12 hours during any working day. Formulate the problem as an LPP.
2 Marks

ANSWER

clear
clc
c=[-3;-4;0;0];
A=[2 3 1 0;1 2 0 1;0 0 1 0;0 0 0 1];
b=[450;720;0;0];
Aeq=zeros(4,4);
Beq=zeros(4,1);
lb=zeros(4,1);
ub=[inf,inf,inf,inf];
x=linprog(c,A,b,Aeq,Beq,lb,ub)
2. A manufacturer produces three types of plastic fixtures. The time required for molding, trimming, and packaging is given in Table .1. (Times are given in hours per dozen
fixtures.) 2 Marks

Process Type A Type B Type C Total time available

Molding 1 2 12,000

Trimming 1 4,600

Packaging 2,400

Profit 110 160 150

How many dozen of each type of fixture should be produced to obtain a maximum profit? 2 marks

ANSWER

clear
clc
c=[-110;-160;-150;0;0;0];
A=[1 2 1.5 1 0 0;2/3 2/3 1 0 1 0; 0.5 1/3 0.5 0 0 1;0 0 0 1 0 0;0 0 0 0 1 0;0 0 0 0 0 1];
b=[12000;4600; 2400; 0;0;0];
Aeq=zeros(6,6);
Beq=zeros(6,1);
lb=zeros(6,1);
ub=[inf,inf,inf,inf,inf,inf];
x=linprog(c,A,b,Aeq,Beq,lb,ub)
Tutorial 2: Efficiency Problems:Use Functions
3. Calculate the maximum theoretical thermal efficiency of a coal-fired power station that heats steam to 500°C and cools it in a condenser at 35°C

The temperature of the gases in a car engine during combustion is 1700°C. The exhaust is expelled at 75°C. Calculate the maximum theoretical thermal efficiency of the
engine. 2 Marks

ANSWER
clear
clc
t1=1700 +273; % Temperature of steam in Kelvin
t2=75 + 273; % Temperature of the condenser in Kelvin
eff=1-(t2/t1); % maximum efficiency
efficiency=eff*100;
fprintf('Highest efficiency in percentage is %8.2f\n', efficiency)
4. Create a function in MATLAB to calculate the ideal air standard cycle efficiency based on the Otto cycle for a gas engine with a cylinder bore of 50 mm, a stroke of 75 mm
3
and a clearance volume of 21.3 cm .Take the inputs as Cylinder bore, clearance volume and get the outputs of compression ratio and efficiency. The calculations are done
as follows in the snapshot 2 Marks

Objective: To find compression ratio and efficiency

Analytical Solution


As we know swept volume =  502  75  147200 m3 =147.2 cm3
4
3
Therefore the total cylinder volume = 147 + 21.3 = 168 cm

168.5
Compression ratio , rv = = 7.914 / 1
21.3

1
Then use the equation   1  1
rv
1
= 1-
7.9140.4
= 56.3 %

Solve using MATLAB function :It should be generalized program, any values inputed should give result

ANSWER
clc
% Inputs
a=input('Cylinder bore in mm');
b=input('Clearance Volume in cm');
% Calculations of cylinder volume
vol=(pi/4)*a^2*75/1000+b;
% Calculations of Compression ratio
CR=vol/b;
eff=(1-(1/CR^0.4))*100;
fprintf('compression ratio is %8.2f\n', CR);
fprintf('Efficiency is %8.2f\n', eff);
Tutorial 3: Gauss Siedel Iteration

5. Use Gauss Siedel iteration to find the values in x1, x2, x3 2 marks

4x1 + x2 - x3 = 3
2x1 + 7 x2 + x3 = 19
x1 - 3 x2 +12 x3 = 31

ANSWER
% find the values of given equation using gauss seidel iteration method
clear
clc;
f1=@(x,y,z) (1/4)*(3-y+z);
f2=@(x,y,z) (1/7)*(19-2*x-z);
f3=@(x,y,z) (1/12)*(31-x+3*y);
x=0;
y=0;
z=0;
for i=1:5
f1(x,y,z);
f2(x,y,z);
f3(x,y,z);
x= f1(x,y,z);
y= f2(x,y,z);
z= f3(x,y,z);
end
disp(x)
disp(y)
disp(z)
o o
6. A long column with thermal conductivity k = 1 W/m K is maintained at 500 K on three surfaces while the remaining surface is exposed to a convective environment with h =
10 W/m  K and fluid temperature T. The cross sectional area of the column is 1 m by 1 m. Using a grid spacing x = y = 0.25 m, determine the steady-state temperature
2o

distribution in the column and the heat flow to the fluid per unit length of the column.

0
Refer Slide of Lecture 13 for Detailed explanation, run the code in matlab and Plot the temperature profile by changing the temperature on three surface to 510 k
2 Marks
Tutorial 4: Polynomials Solving

7. The following points are given

x 1 2.2 3.7 6.4 9 11.5 14.2 17.8 20.5 23.2

Y 12 9 6.6 5.5 7.2 9.2 9.6 8.5 6.5 2.2

(a) Fit the data with a first-order polynomial. Make a plot of the points and the polynomial.
(b) Fit the data with a second-order polynomial. Make a plot of the points and the polynomial.
(c) Fit the data with a third-order polynomial. Make a plot of the points and the polynomial.
(d) Fit the data with a fifth-order polynomial. Make a plot of the points and the polynomial.

ANSWER
clear
clc
x=[1 2.2 3.7 6.4 9 11.5 14.2 17.8 20.5 23.2];
y=[12 9 6.6 5.5 7.2 9.2 9.6 8.5 6.5 2.2];
% Solution of (a)
p1=Polyfit(x,y,1);
plot(x,y)
title('first-order polynomial')
% Solution of (b)
P2=Polyfit(x,y,2);
plot(x,y)
title('second-order polynomial')
% Solution of (c)
P3=Polyfit(x,y,3);
plot(x,y)
title('Third-order polynomial')
% Solution of (d)
P=Polyfit(x,y,4);
plot(x,y)
title('fifth-order polynomial')
8. The temperature (in degrees Fahrenheit) was recorded every three hours for a day at a particular location. Using a 24-hour clock where midnight is
0, for example,
the data might be

Time 0 3 6 9 12 15 18 21
Temp 55.5 52.4 52.6 55.7 75.6 77.7 70.3 66.6

a. Create vectors for the data.

b. Plot the data.
c. Find a curve that fits the data.
d. At what time(s) was it 60 degrees? 65 degrees? 2 Marks

ANSWER
clear
clc
% Solution for a
x=[0:3:21];
y=[55.5 52.4 52.6 55.7 75.6 77.7 70.3 66.6];
p=polyfit(x,y,1);
curve=polyval(p,x);
xlabel('Time')
ylabel('Temperature')
plot(x,y, 'ro',y,curve)
title('Time Vs Temp');
% Values at 60 & 65.
polyval(p,60)
polyval(p,65)
Tutorial 5: Control Problems

9. If the transfer function is given below.Plot he Step response and Impulse response for the transfer function.

( )

2 Marks
ANSWER

% Step response
G=tf([2],[1 3 2]);
step(G)
hold on;
% Impulse response
% considering zeta = 0,0.2,0.4,0.6,…..1.8
t=0:0.2:10;
zeta=[0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8];
for n=1:10;
num=[0 0 1];
den=[1 3*zeta(n) 2];
[y(1:51,n),x,t]=step(num,den,t);
End
Plot(t,y)
grid;
title('Imuplse response');
ODE 45 Questions

10. Use ode45 to solve the initial value problem for a first order differential equation:

2 Marks

ANSWER

function f=myfun1(t,y)
f=-t*y/sqrt(2-y^2);
[tv1 f1]=ode45('myfun1',[0 5],1);
Plot(tv1,f)
11. Given dy/dx= 1+y2 , where y=0 when x=0 find y(0.2), y(0.4) and y(0.6) compare Runga Kutta 4th order with euler and modified euler method
using MATLAB 2 Marks

ANSWER

%Euler's Method
x=zeros(4,1);
y=zeros(4,1);
x(1)=0;
y(1)=1;
for i= 1:3
x(i+1)=x(i)+.2;
y(i+1)=y(i)+.2.*((y(i).^2)+1);
end
Eulers_Method_Result1=[x,y]
plot(x,y,'-*r');
hold on
% Moified Euler's Method
x=zeros(4,1);
y=zeros(4,1);
f=@(y,x)(1+y^2);
y(1)=1;
x(1)=0;
h=.2;
for i=1:3
x(i+1)=x(i)+h;
y(i+1)=y(i)+h*f(y(i),x(i));
y(i+1)=y(i)+((f(y(i),x(i))+f(y(i+1),x(i+1)))*h/2);
end
Modified_Eulers_Method_Results=[x,y]
plot(x,y,'-*b');
%Runge kutta 4th order method
x=zeros(4,1);
y=zeros(4,1);
x(1)=0;
y(1)=1;
h=.2;
f=@(x,y)y^2+1;
for i=1:3
x(i+1)=x(i)+h;
k1 = f(x(i),y(i));
k2 = f(x(i)+0.5*h,y(i)+0.5*h*k1);
k3 = f((x(i)+0.5*h),(y(i)+0.5*h*k2));
k4 = f((x(i)+h),(y(i)+k3*h));
y(i+1) = y(i) + (1/6)*(k1+2*k2+2*k3+k4)*h;
end
Runge_Kutta_4th_Order_Method_Result3=[x,y]
plot(x,y,'-*g');
xlabel('X')
ylabel('Y')
title('Comparison of RK4,Eulers and Modified Euler Methods')
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018

Comprehensive Examination
(EC-3 REGULAR)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 40%
No. of Pages =3
Duration : 3 Hours
No. of Questions = 10
Date of Exam : 22/04/2018 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1. a = 2; b = 3; c = 4; d = 5;
y1 = b/a*c + b^a - d;
y2 = b/a *(b -a) ^( b -a )*( b -a ); [4]

Solution:

y1 = b/a*c + b^a - d;
(a) y1= 3/2*4+ 3^2-5

Precedence from left to right so [/ & * same precedence]

y1=3/2=1.5*4=6 + 9-5=10

(b) y2 = b/a (b -a) ^( b -a )( b -a );

Precedence from left to right so

y2= 3/2 *(3-2)^(3-2)*(3-2)
y2=1.5*1^(1*1)
y2=1.5

Q.2. If the function f(x)=𝑥 3 − 12𝑥 + 40.25𝑥 − 36, Use the fminbnd with the interval 2 ≤ 𝑥 ≤
7. The format of fmind is x = fminbnd (function,xl,x2) , the value x is where the function
has minimum value and x1 and x2 is the interval. [4]

f=@(x) x^3-12*x+40.25*x-36
x=fminbnd(f,2,7)

{ x = 2.0000 in MATLAB}
Q.3. Write a MATLAB script that doubles the positive elements are divisible by 2 or 3 else if
the elements are negative assign them as zero. Take the same vector as a=[4,6,-9,1,3,2,1]
[4]
% Function script to double positive elements and assign zero to negative elements

function divisible=vector(a)

for i=1:7
if a(i) > 0 && a(i)%2==0 ||a(i)%3==0
divisible(i)=a(i)*a(i);
else
divisible(i)=0;
end
end

[output optional in all prgrams]

Output

vector([1,2,4,-8,4,2,3])

ans =

1 4 16 0 16 4 9

Q.4. An experiment is conducted to find the value of universal gas constant R for 1 mol of gas.
The gas follows ideal gas law PV=nRT, V in Volume in liters, P-atm, T is in 0K and n is
the number of mol of gas. Use the data below to determine R, by plotting V versus T/P and
fit to the linear equation. { Write only commands } [4]
V 0.8 0.6 0.4 0.2
T 28 39 57 69.5
P 1.7 1.98 2.5 4

n=1;
V=[0.8 0.6 0.4 0.2];
T=[28 39 57 69.5];
P=[1.7 1.98 2.5 3 4];
TdP=(T+273)./P;
p=polyfit(TdP,V,1);
R=p(1)/n
TdPplot=linspace(200,70,50)
Vplot=p(1)*TdPplot+p(2);
plot(TdP,V,'o',TdPplot,Vplot)
axis([50 220 0 1])
xlabel('T/P (K/atm)')
ylabel('V (L)’)

0.006396799854418
Q.5. The resistance R in ohms of a conductor is given by R= E/I, where E is the potential in volts
and I is the current in amperes. Write a script that will (a) Call a function to prompt the
user for the potential and the current. And (b) Call a function that will print the resistance;
this will call a sub function to calculate and return the resistance. [4]

function r=vector (n1,n2)

r= n1/n2;
disp(r);
end

(ii)
% subfunction
n1=input(‘value of voltage’);
n2=input(‘value of current’);
fprintf(‘The resistance value is %.2f’,vector(n1,n2));

Q.6. The Wind Chill Factor (WCF) measures how cold it feels with a given air temperature (T,
in degrees Fahrenheit) and wind speed (V, in miles per hour). One formula for the WCF is

WCF = 35.7 + 0.6𝑇 − 35.7 ( 𝑉 0.16 ) + 0.43𝑇(𝑉 0.16 )

Create a table showing wind chill factors for temperatures ranging from –20 to 55 in steps
of 5, and wind speeds ranging from 0 to 55 in steps of 5. Write this to a file ‘wcftable.dat’.
[4]
%Script to print table of wind chill factors
% Print column headers
fprintf('%45s\n ', 'Wind Speeds')
for v = 0:5:55
fprintf('%7d', v)
end
fprintf('\nTemp\n')
for t = -20:5:55
fprintf('%3d', t)
for v = 0:5:55
fprintf('%7.1f',wcf(t,v));
end
fprintf('\n')
end

function WCF=wcf(t,v)
WCF=35.7+0.6*t-35.7*v.^0.16+0.43*v.^0.16;
end

TA ZC164 (EC-3 Regular) Second Semester 2017-2018 Page 2

Q.7. In product design, it is useful to gauge how important different features of the product would
be to potential customers. One method of determining which features are most important is
a survey in which people are asked for different features: “Is this feature important to you?”
The number of potential customers who responded Yes is then tallied. For example, a
company conducted such a survey for 10 different features; 200 people took part in the
survey. The data was collected into a file that might look like this:

1 2 3 4 5 6 7 8 9 10
30 83 167 21 45 56 55 129 69 55

A Pareto chart is a bar chart in which the bars are arranged in decreasing values. The bars
on the left in a Pareto chart indicate which are the most important features. Create a data
file, and then a subplot to display the data with a bar chart organized by question on the
left and a Pareto chart on the right [4]

Create a data file in command window

>> Questions_respondents=[1:10;30 83 167 21 45 56 55 129 69 55]

>> save pareto.dat Questions_respondents –ascii

>> questions=pareto(:,1)
>> responses=pareto(2,:)
>> bar(res)

Q.8. Find out the specific speed of a turbine of 10 MW capacity working under a head of 500m
and having the normal working speed of 300 RPM.

Write a function in MATLAB where if you enter the above three parameters the specific
speed should calculate. Suppose the 10 MW is replaced by 100 MW and so on the program
should be generalized for all inputs [4]

function specific_speed=turbinespeed(h,P,N)
specificspeed=N*sqrt(P)/h^(5/4);
end

h=input(‘Head’);
P=input(‘Power input’);
N=input(‘Specific Speed’);
fprintf(‘the specific speed of turbine is %f’,turbinespeed(h,P,N);

Q.9. Create a mesh, surface plot and contour plot of the function z=𝑒 𝑥+𝑖𝑦 for the interval -1≤
𝑥 ≤ 1 𝑎𝑛𝑑 − 2𝝅 ≤ 𝑦 ≤ 2𝝅. In each case plot the real part of z versus x and y [4]

close all
clear all
clc

% Number of gridpoints in y- and x-direction

n = 60;
m = 100;

% Create vectors with x- and y-values

x = linspace(-1,1,m);
y = linspace(-2*pi,2*pi,n);

% Product a grid of all possible

% combination of x- and y-values.

[X,Y] = meshgrid(x,y);
% compute function

Z = real(exp(X + sqrt(-1)*Y));
Clf

surf(X,Y,Z);
shading interp
view(3)
xlabel(’x’);
ylabel(’y’);
label(’z’);

figure(2)
mesh(X,Y,Z+4,Z);

contour(X,Y,Z)

Q.10. A long column with thermal conductivity k = 1.2 W/moK is maintained at 500oK on three
surfaces while the remaining surface is exposed to a convective environment with h = 10
W/m2oK and fluid temperature T. The cross sectional area of the column is 1 m by 1 m.
Using a grid spacing x = y = 0.25 m, determine the steady-state temperature distribution
in the column and the heat flow to the fluid per unit length of the column.

Understand the Process by energy balance and Determine the Node temperatures using
MATLAB code [4]
The equation is divided by k(xy1) and simplified to
***********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2018-2019

Comprehensive Examination
(EC-3 REGULAR)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 40%
No. of Pages =3
Duration : 3 Hours
No. of Questions = 10
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1. Evaluate the values of i and m through hand calculation for the following script [4]

v = [3 1 5];
i = 1;
for j = v
i = i + 1;
if i = = 3
i = i + 2;
m = i + j;
end
end

Solution:
Initially v is assign the value [3 1 5] that is 1 row and 3 columns
i=1;
j will be assigned v vector
i=2
condition fails goes to loop
i=3 now
i=i+2 will make value of i=5
j will get 1 value as initially i=1
therefore m= then value of m will
m=5+1=6

Q.2. Create a multiplication table from 1 to 10 in MATLAB and write it to a spreadsheet. [4]

a = 1:10;
filename='multiplication.xlsx';
A = a .* a'
xlswrite(filename,A)
Please also check if somebody does alternate steps , give partial marks

Q.3. For the following data set x=[1 2 3 4 5 6 7 8] and y=[ 25 26 78 11 13 15 18 21] Create the
bar char, pie chart, stem chart. All the charts should appear in one figure window only. Use
subplot command [4]

x=[1 2 3 4 5 6 7 8]
y=[ 25 26 78 11 13 15 18 21]
x=[1 2 3 4 5 6 7 8]
y=[ 25 26 78 11 13 15 18 21]
subplot(2,2,1)
bar(x,y)
subplot(2,2,2)-/
pie(x,y)
subplot(2,2,3)
stem(x,y)

Q.4 The following program (writer.m) invites you to enter any number of names and marks from the
command line, and writes them as binary data to a file. To terminate the process just hit Enter for
the next name. Explain what task will each line do and overall output [4]
namelen = 10; % 10 bytes for name
fid = fopen(’marks.bin’, ’w’); % open for write only
str = ’?’; % not empty to start
while ~isempty(str)
str = input( ’Enter name: ’, ’s’ );
if ˜isempty(str)
if length(str) > namelen
name = str(1:namelen); %only first ten chars allowed
else
name = str;
name(length(str)+1:namelen) = ’ ’; %pad with blanks if.. too short
end
fwrite(fid, name);
mark = input( ’Enter mark: ’ );
fwrite(fid, mark, ’float’); % 4 bytes
for mark
end
end
fclose(fid);

Solution

namelen = 10; % 10 bytes for name

%*+/namelen is assigned length of 10
fid = fopen(’marks.bin’, ’w’); % open for write only
f
i
d
i
s
f
i
l
e
i
d
e
n
t
i
f
i
e
r
a
n
d
f
o
p
e
n
w
i
l
l
o
p
e
n
m
a
r
k
s
.
b
i
n
f
i
l
e
f
o
r
w
r
i
t
i
n
g
str = ’?’; % not empty to start str
is string to enter name of student
while ~isempty(str)
~ not operator to make sure string
str = input( ’Enter name: ’, ’s’ );
if it is empty and length is more
if ˜isempty(str)
it will assigned name else name=str(1:namelen)
length(str) > namelen
name = str(1:namelen); %only first ten chars allowed
else
name = str;
name(length(str)+1:namelen) = ’ ’; %pad with blanks if.. too short
end
fwrite(fid, name); will write
all the names to the file
mark = input( ’Enter mark: ’ ); will again
enter marks
fwrite(fid, mark, ’float’); % 4 bytes for mark
…then file will be closed
for mark
end
end
fclose(fid);

Q.5 Write a MATLAB script that doubles the positive elements are divisible by 2 or 3 else if the
elements are negative assign them as zero. Take the same vector as a=[4,6,-9,1,3,2,1]

[4]

function divisible=vector(a)

for i=1:7
if a(i) > 0 && a(i)%2==0 ||a(i)%3==0
divisible(i)=a(i)*a(i);
else
divisible(i)=0;
end
end

Q.6 If we define a structure as follows

meteo = struct('Site',{'Adelaide', 'Sydney'},'Time',{2.3 4},'Temperature',{24 19},'Pressure',{1023 1015})

(a) >> meteo(1)

(b) >> meteo(2)
(c) >> meteo(2).Temperature
(d) >> meteo(2).Pressure [4]

(a) Time: 2.3000

Temperature: 24
Pressure: 1023

(b) Site: 'Sydney'

Time: 4
Temperature: 19
Pressure: 1015
(c ) ans= 19

Q.7 The coefficient of variation is useful when comparing data sets that have quite different
means. The formula is CV = (standard deviation/mean) * 100%. A history course has two
different sections; their final exam scores are stored in two separate rows in a file. For
example

99 100 95 92 98 89 72 95 100 100

83 85 77 62 68 84 91 59 60

Create this data file, read the data into vectors, and then use the CV to compare the two
sections of this course. [4]

>> firstfile=[99 100 95 92 98 89 72 95 100 100]

>> secondfile=[83 85 77 62 68 84 91 59 60 ]
>> save pareto.dat firstfile secondfile –ascii
>> cv1=std(firstfile)/mean(firstfile)*100
>> cv2=std(secondfile)/mean(secondfile)*100

Q.8 Make a 3 D surface plot of the function z=0.5x2+0.5y2 in the domain −2 ≤ 𝑥 ≤ 2 𝑎𝑛𝑑
−2 ≤ 𝑦 ≤ 2 [4]

clear;

clc;
x=-2:0.1:2;
y=-2:0.1:2

z=0.5.*x^2+
0.5.*y^2;
plot3(x,y,z
)

title('Plot 3 Dimensional
figure');
xlabel('x axis');

ylabel('y axis');
zlabel('z axis');
grid on
Q.9 If there are two strings str1= ‘Non’ str2= ‘specific’. What will be the result of the following
command? [4]

(i) >> a=strcat(str1,str2)

(ii) >> a=strvcat(str1,str2)
(iii) >> strcmp(str1,str2)
(iv) >> strncmp(str1,str2,3)

Solution:

(i) a = 'Nonspecific
(ii) a=
2×8 char array 'Non '
'specific'
(iii) logical=0
(iv) logical=0

Q.10 Create a function in MATLAB to calculate the ideal air standard cycle efficiency based on the
Otto cycle for a gas engine with a cylinder bore of 40 mm, a stroke of 70 mm and a clearance
volume of 19.3 cm3.Take the inputs as Cylinder bore, clearance volume and get the outputs of
compression ratio and efficiency. The calculations are done as follows in the snapshot
𝜋
* Swept Volume= 4 × 402 × 70 *Total Cyl vol=Swept vol+clearance vol
* compression ratio is total cyl vol/clearance vol
* efficiency = 1- (1/comp ration ^0.4) *100

function Efficiency=ottocycle(d,S,VC)
% d is bore diameter ENTER IN MM
% S is stroke ENTER IN MM
%VC is clearance volume ENTER IN MM3
VS=(pi*d^2*S/4)*(1/1000) ;% conversion to cm
disp(VS);
r=(VC+VS)/VC; % no difference in unit as units are same
disp(r);
Efficiency=(1- (1/(r^0.4)))*100;
Best of Luck

***********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Mid-Semester Test (EC-2 Make-up)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 30% No. of Pages =2
Duration : 2 Hours No. of Questions = 7
Date of Exam : 18/03/2018 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.

Q.1. Short Questions [2 x 5 = 10]

(a) Create a 2 × 4 matrix of random integers, each in the range from 10 to 30

randi([10,30],2,4)

(b) If mat = [4 3 1; 2 5 6] what is the output of mat = [2:4; 3:5]

mat =

2 3 4
3 4 5

linspace(-3,-15,5)

(d) Create a 4 × 2 matrix of all zeros and store it in a variable. Then, replace the second
row in the matrix with a 3 and a 6.

a=zeros(4,2)
a(2,:)=[3,6]

(e) Write a command to find the Eigen value of a 2*2 matrix z.

z=[1,2;3,4]
eig(z)

Q.2. Evaluate the value of c if c=a.\b , a and b are matrices shown below [2]

 3 1   1 3
a  b 
  1 2  2 2

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 1 of 4

B vector will be divided by a

-0.3333 3.0000
-2.0000 1.000

Q.3. What is the output of the following output statement [2]

i) fprintf(‘%3.2f\n’, 1234.5678)

Solution:

1234.57 it will take only 2 decimal points

ii) for i=0:10:100

if mod(i,10) ==0
disp(i)
else
continue;
end
end

Solution:
It will print 0 to 100 in the gap of 10 such as 0,10,20,20,40,50,60,70,80,90,100

kx 2
Q.4. The potential energy stored in a spring is where k is spring constant and x is
2
compression in the spring. The force required to compress the spring is kx. The following
table gives the data for five springs
Spring No 1 2 3 4 5
Force (N) 11 7 8 10 9
Spring Constant k (N/m) 1000 600 900 1300 700

Use MATLAB to find (i) compression x in each spring and (ii) Potential energy stored in
each spring

Solution:

>> % force f
>>f=[11, 7, 8, 10, 9]
>>k=[1000, 600, 900,1300,700]
>>x=f./k;
>>pe=k.*x.^2/2;

Give Partial marks for steps

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2 of 4

[4]
TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2

Q.5. Define the vector v = [2 3 4 5 6]. Then use the vector in a mathematical expression to create
the following vectors: {use single command for each expression} [4]
(a) a = [4 6 8 10 12 ]
(b) b = [8 27 64 125 216 ]
(c) c = [22 33 44 55 66 ]
(d) d = [1 1.5 2 2.5 3]

Solution:

a= v.*2;
b=v.^3;
c=v.^v;
d=a.^0.25

Q.6. Clouds are generally classified as high, middle, or low level. The height of the cloud is the
determining factor, but the ranges vary depending on the temperature. For example, in
tropical regions the classifications may be based on the following height ranges (given in
feet):
low 0–6500
middle 6500–20000
high > 20000
Write a script that will prompt the user for the height of the cloud in feet, and print the
classification. [4]

Script

% Cloud Classification
x=input('Enter the height of the cloud in feet')
if (x>=0 & x<6500)
disp('Low cloud')
elseif (x>=6500 & x<20000)
disp('Medium cloud')
elseif (x>2000)
disp('high cloud')
end

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 3 of 4

Q.7. Write a function that will receive a matrix as an input argument, and will calculate and
return the overall average of all numbers in the matrix. Use loops, not built-in functions, to
calculate the average [4]

Solution

function average=overall_avg(x)
n=length(x);
sum=0;
for i=1:n
sum=sum+i;
average=sum/n;
end

>> overall_avg([1,2,3,4,5,6,7,8,9,10])

ans =

5.5000

********

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 4 of 4

Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Comprehensive Test (EC-2 Regular)

Course No. : MTU ZC221

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 50% No. of Pages =2
Duration : 3 Hours No. of Questions = 10
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

1. Evaluate the output of the following MATLAB expressions, if a is the input vector containing the
following elements a=[1,2,3;4,5,6]

a. Convert the Vector into a matrix of 2 rows and 3 columns using reshape
b. Multiply all elements of second row and second column by two
c. Add the third row which contains average of all the columns
d. Add the fourth column containing average of all rows
e. Find the min value for each row and each column

2. The demand for water during a fire is often most important factor in the design of distribution
storage tanks and pumps. For communities with population less than 1,00,000 the demaond Q in
gallons/min can be calculated by

𝑄 = 1020√𝑃 (1 − 0.01√𝑃) Where P is population in thousands. Set up vector P that starts at 10

and increments by 10 upto 200. Use element by element computations to determine Q for each
population in P

3. Using commands in MATLAB to plot a 3-D figure of the trajectory of a particle defined by the
following set of equations for x,y, and z.Use commands for proper labelling and title

x=t
y=t*cos(t)
z= e 0.2t

4. Consider the following steam table and determine the values using polyfit, polyval and interp1
commands using MATLAB
T, 0.01 5 10 15 20 25 30 35 40 45 50 55
deg C
P, bar 0.0061 0.0087 0.0123 0.0170 0.0234 0.0317 0.0425 0.0563 0.0738 0.0959 0.1.235 0.1576

a. Determine the pressure at 27 deg C using Interp1 (Refer Q.4)

b. Determine the temperature at 0.0465 bar
c. Determine the pressure at 70 deg C
d. Determine the temperature at 0.20 bar pressure
e. Determine the relationship by polynomial expression between T and expression where T is in 0K

5. The operations manager of a plant that manufactures tires wants to compare the actual inner
diameters of two grades of tires, each of which is expected to be 575 millimeters. A sample of five
tires of each grade was selected, and the results representing the inner diameters of the tires,
ranked from smallest to largest, are as follows:

Grade X 568 570 575 578 584

Grade Y 573 574 575 577 578

Create a data file tire.xls in Excel and calculate the mean, median and std. dev in MATLAB and
export it to excel or write the computed values to excel

6. Construct an equivalent statement using single if statement in MATLAB for the following block
(a) If x<y
If z<10
W=x*y*z

(b) Accepts marks of the students as p and if p is greater than 50 print pass else print fail

7. Given the system of equations

3x1  7x 2  13x3  76
x 1  5x 2  3x 3  28
12x1  3x 2 - 5x 3  1
find the solution using Gauss-Seidal method. Use x1 , x2 , x3   1 0 1 as the initial guess.
𝑑𝑥
8. Consider the example 𝑑𝑡 = 3𝑒 −𝑡 with initial condition x(0)=0, Solve the following equation using
ODE45 in MATLAB

9. Create a function in MATLAB to calculate the ideal air standard cycle efficiency based on the Otto
cycle for a gas engine with a cylinder bore of 50 mm, a stroke of 75 mm and a clearance volume
of 21.3 cm3.Take the inputs as Cylinder bore, clearance volume and get the outputs of
compression ratio and efficiency. The calculations are done as follows in the snapshot

10. Explain the significance of MATLAB in your organization with specific examples, can you think of
a problem { Process control, optimization anything} in your department and try to solve in
MATLAB. Present your case for solving any one problem with stepwise methodology
TA ZC164 EC-3R SECOND SEM 2017-2018-Solutions.docx
TA ZC164 EC-3M SECOND SEM 2016-2017_revised
SOLUTIONS.docx TA ZC164 EC-3 M I SEM 2017-2018.docx TA ZC164
EC-3 I 2016-2017.docx
TA ZC164 EC-2M SECOND SEM 2017-2018 solutions.docx
TA ZC164 EC-2M I 2017-2018-Corrected.docx
TA ZC164 EC-2M FIRSTSEM 2017-2018-CorrecteD
SOLUTIONS.docx TA ZC 164 EC-3 17-18-Makeup.docx TA ZC 164
EC-3 17-18 Solutions.docx
POWTP ZC164 EC-3R FIRSTSEM 2017-2018.docx
POW TP 164 MAKE UP PAPER.docx
CUMMINS-MID SEM-SOLUTIONS.docx
CUMMINS-COMPRE.docx
CUMMINS-COMPRE Solutions-Final.docx
CUMMINS-COMPRE SOLUTIONS.docx
TA ZC164 Mid Sem -first Semester 2018-19.docx
TA ZC164 EC-II 2017-2018-Solutions.docx
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018

Comprehensive Examination
(EC-3 REGULAR)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 40%
No. of Pages = 3
Duration : 3 Hours
No. of Questions = 10
Date of Exam : 22/04/2018 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1. a = 2; b = 3; c = 4; d = 5;
y1 = b/a*c + b^a - d;
y2 = b/a *(b -a) ^( b -a )*( b -a ); [4]

Solution:

y1 = b/a*c + b^a - d;
(a) y1= 3/2*4+ 3^2-5

Precedence from left to right so [/ & * same

precedence] y1=3/2=1.5*4=6 + 9-5=10

(b) y2 = b/a (b -a) ^( b -a )( b -a );

Precedence from left to right so

y2= 3/2 *(3-2)^(3-2)*(3-2)
y2=1.5*1^(1*1)
y2=1.5
Q.2. If the function f(x)= 3 − 12 + 40.25 − 36, Use the fminbnd with the interval 2 ≤
7. The format of fmind is x = fminbnd (function,xl,x2) , the value x is where the function
≤
has minimum value and x1 and x2 is the interval. [4]

f=@(x) x^3-12*x+40.25*x-36
x=fminbnd(f,2,7)

function divisible=vector(a)

for i=1:7
if a(i) > 0 && a(i)%2==0 ||a(i)%3==0
divisible(i)=a(i)*a(i);
else
divisible(i)=0;
end
end

[output optional in all prgrams]

Output

vector([1,2,4,-8,4,2,3])

ans =

1 4 16 0 16 4 9

Q.4. An experiment is conducted to find the value of universal gas constant R for 1 mol of gas. The
gas follows ideal gas law PV=nRT, V in Volume in liters, P-atm, T is in 0K and n is the
number of mol of gas. Use the data below to determine R, by plotting V versus T/P and
fit to the linear equation. { Write only commands } [4]
V 0.8 0.6 0.4 0.2
T 28 39 57 69.5
P 1.7 1.98 2.5 4

0.006396799854418
Q.5. The resistance R in ohms of a conductor is given by R= E/I, where E is the potential in volts
and I is the current in amperes. Write a script that will (a) Call a function to prompt the user
for the potential and the current. And (b) Call a function that will print the resistance;
this will call a sub function to calculate and return the resistance. [4]

function r=vector (n1,n2)

r= n1/n2;
disp(r);
end

(ii)
% subfunction
n1=input(‘value of voltage’);
n2=input(‘value of current’);
fprintf(‘The resistance value is %.2f’,vector(n1,n2));

Q.6. The Wind Chill Factor (WCF) measures how cold it feels with a given air temperature (T,
in degrees Fahrenheit) and wind speed (V, in miles per hour). One formula for the WCF is
WCF = 35.7 + 0.6 − 35.7 ( 0.16) + 0.43 ( 0.16)

Create a table showing wind chill factors for temperatures ranging from –20 to 55 in steps of
5, and wind speeds ranging from 0 to 55 in steps of 5. Write this to a file ‘wcftable.dat’.
[4]
%Script to print table of wind chill
factors % Print column headers
fprintf('%45s\n ', 'Wind Speeds')
for v = 0:5:55
fprintf('%7d', v)
end
fprintf('\nTemp\n' )
for t = -20:5:55
fprintf('%3d', t)
for v = 0:5:55
fprintf('%7.1f',wcf(t,v));
end
fprintf('\n')
end

function WCF=wcf(t,v) WCF=35.7+0.6*t-

35.7*v.^0.16+0.43*v.^0.16; end

TA ZC164 (EC-3 Regular) Second Semester 2017-2018 Page 2

Q.7. In product design, it is useful to gauge how important different features of the product
would be to potential customers. One method of determining which features are most
important is a survey in which people are asked for different features: “Is this feature
important to you?” The number of potential customers who responded Yes is then tallied.
For example, a company conducted such a survey for 10 different features; 200 people
took part in the survey. The data was collected into a file that might look like this:

1 2 3 4 5 6 7 8 9 10
30 83 167 21 45 56 55 129 69 55

A Pareto chart is a bar chart in which the bars are arranged in decreasing values. The bars on
the left in a Pareto chart indicate which are the most important features. Create a data file,
and then a subplot to display the data with a bar chart organized by question on the
left and a Pareto chart on the right [4]
Create a data file in command window

>> Questions_respondents=[1:10;30 83 167 21 45 56 55 129 69 55]

>> save pareto.dat Questions_respondents –ascii

>> questions=pareto(:,1)
>> responses=pareto(2,:)
>> bar(res)

Q.8. Find out the specific speed of a turbine of 10 MW capacity working under a head of 500m
and having the normal working speed of 300 RPM.

function specific_speed=turbinespeed(h,P,N)
specificspeed=N*sqrt(P)/h^(5/4);
end

h=input(‘Head’);
P=input(‘Power input’);
N=input(‘Specific Speed’);
fprintf(‘the specific speed of turbine is %f’,turbinespeed(h,P,N);

Q.9. Create a mesh, surface plot and contour plot of the function z= + for the interval -1≤
≤ 1− 2 ≤ ≤ 2 . In each case plot the real part of z versus x and y [4]

close all
clear all
clc

% Number of gridpoints in y- and x-direction

n = 60;
m = 100;

% Create vectors with x- and y-values

x = linspace(-1,1,m);
y = linspace(-2*pi,2*pi,n);
% Product a grid of all possible
% combination of x- and y-values.

[X,Y] = meshgrid(x,y);
% compute function

Z = real(exp(X + sqrt(-1)*Y));
Clf

surf(X,Y,Z);
shading interp
view(3)
xlabel(’x’);
ylabel(’y’);
label(’z’);

figure(2)
mesh(X,Y,Z+4,Z);

contour(X,Y,Z)

Comprehensive Examination
(EC-3 Makeup)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 40%
No. of Pages = 3
Duration : 3 Hours
No. of Questions = 9
Date of Exam : 23/04/2017 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1 Short Questions [1 x 5 = 5]

i. Complete the code to integrate the function, y(x) with respect to x.

syms x
y=x^2+6*x+5
a. int(y,x) b. int(x,y)
c. integrate(y,x) d. trapz(x,y)

Solution:

a. int(y,x)

ii. Which MATLAB command generates 50 uniform random numbers between -2 and
+2?
a. 10*rand(50)-2 c. 4*randn(50,1)-2
b. 2*rand(50,1)-2 d. 2*randn(50)-2

Solution:

b. 2*rand(50,1)-2
iii. Write the MATLAB commands to find a determinant for a 3*3 matrix

Solution:
>> a=[1,2,3;4,5,6;7,8,9]
>> det(a)

iv Which MATLAB command to find the covariance of x and y in MATLAB?

Solution:
>> x=[1,2,3]
>> y=[2,4,6]
>> z=cov(x,y)

v. What sign will be prefix for a function if you want to find a differential equation or a
integral for any function

Solution
>>@

Q.2 Determine the output of the following [2]

(i) >> grades = [85 70 100 95 80 91];
>> grade_index = [3 5 1 6 4 1];
>> grades(grade_index)

Solution

ans =

100 80 85 91 95 85

(ii) >>a=eye(3);
>>sparse(a)
ans =

(1,1) 1
It will print non zero elements
(2,2) 1 or
(3,3) 1
TA ZC164 (EC-3 Make-up) Second Semester 2016-2017 Page 2

Q.3 Examine the following (i) and (ii) for loops and determine the value of b at the end of each
loop and also the number of times each loop executes [4]

(i) (ii)
b=0; b=0;
for i:-10:10 for i=10:-2:4
b=b+1; for j=2:2:i
end if j==6
break
end
b=b+j;
end
End

Solution

(a) It will print

values of b from
1 to 20
b=1, b=2,---------
b=20

(b) (b) It will print

values of b from
2 to 24
b=2, b=4,---------
b=24

Q.4 Write a function called geomser that will receive values of r and n, and will calculate and
return the sum of the geometric series: 1 + r + r2 + r3 + r4 + ... + rn
The following examples of calls to this function illustrate what the result should be:
>> geomser(1,5)
ans = 6
>> disp(geomser(2,4))
31 [4]

Solution
function s=geomser(r,n)
v = [1,r*ones(1,n)];
p = cumprod(v);
s = sum(p);

end

Q.5 Create a mesh, surface plot and contour plot of the function
z= + for the interval -1≤ ≤ 1− 2 ≤ ≤ 2 . In each case plot the real part of

z versus x and y [4]

% Mesh Plot

close all
clear all
clc

% Number of gridpoints in y- and x-

direction n = 60;
m = 100;

% Create vectors with x- and y-values

x = linspace(-1,1,m);
y = linspace(-2*pi,2*pi,n);

% Product a grid of all possible

% combination of x- and y-values.

[X,Y] = meshgrid(x,y);
% compute function

Z = real(exp(X + sqrt(-

1)*Y)); Clf

surf(X,Y,Z);
shading interp
view(3)
xlabel(’x’);
ylabel(’y’);
label(’z’);

figure(2)
mesh(X,Y,Z+4,Z);

contour(X,Y,Z)

Q.6 The open loop transfer function G(s)H(s) of a control system is [4]
( )( )=
( + 0.5)( 2 + 0.5 + 8)
Plot the root loci for the system using MATLAB

Solution

% Enter coefficients of

Numerator Num=[1];

%Solve for Denominator

Den=[1,1,8.25,4,0]

Rootlocus(num,den)
Q.7 For a project, some biomedical engineering students are designing a device that will monitor a
person’s heart rate while on a treadmill. The device will let the subject know when the target
heart rate has been reached. A simple calculation of the target heart rate (THR) for a
moderately active person is THR=(220−A)*0.6 where A is the person’s age. Write a
function that will calculate and return the THR [4]

Solution

function thr= age(a)

thr=(220-a)*0.6;
end

TA ZC164 (EC-3 Make-up) Second Semester 2016-2017 Page 3

Q.8 Using this polynomial, predict the accumulation over the range of 20 - 30 hours. Finally,
calculate the time derivative of the accumulation over the period 0-10 hours. [4]
Mass of A accumulated as a function of time:

Mass of A accumulated (arbitrary units) Time (hours)

9 1
55 3
141 5
267 7
345 8
531 10

Only write commands and write what we can get as far as analysis

Solution:

>> mass=[9,55,141,267,345,531];
>> time=[1,3,5,7,8,10];
>> coeff = polyfit(time,mass,2)

/* coeff =
5.0000 3.0000 1.0000

/* This is not needed

So, Mass A = 5*(time)2 + 3 * (time) + 1.
The coefficients of the polynomial fit are now stored in the row vector coeff. 8?
To evaluate this polynomial over the range of 20-30 hours we define a n e w time vector (the
independent variable)

>> newtime = [20:30]

>> pred = polyval(coeff,newtime)

Q.9 Consider the system of reactions in a constant volume, constant temperature batch reactor.
A —- > D reaction 1
A + A —— > U reaction 2

where D is a desired product, U is undesired product and where k1 and k2 are the rate
constants for reactions 1 and 2. Let k1, k2 be given parameters, and the initial concentration
of A (ca0) be a design variable. The independent variable is time (tfin), and the dependent
variables are the concentration of the species, ca, cd, and cu. Note that in most design
situations k1 and k2 might be design parameters, adjusted via the temperature.

By writing the mass balance equations over the batch reactor, the system of differential
equations is the following: they are not linearly independent.
d(ca)/dt = -k1(ca) - k2(ca)2
d(cd)/dt = k1 (ca)
d(cu)/dt = k2 (ca)2
Solving this problem in Matlab involves two parts. First, write a function file that describes
the set of ODEs in terms of a single, combined matrix variable (the dependent variable).
Plot the concentration variables with respect to time

Next, in a second, (main) script invoke the ode solver, ode45.

Solution
function dC_dt = exampleode(t,C)
global k1 k2 % variables that we wish to share with the main script
dC_dt(1) = -k1*C(1) -
k2*C(1)*C(1); dC_dt(2) = k1*C(1);
dC_dt(3) = k2*C(1)*C(1);

global k1 k2 % Set parameters,

k1 = 2;
k2 = 1;
ca0 = 2; % ca0 = initial concentration of species A
tfin = 1/3; % all parameters in arbitrary units for purposes of demo % t is time,
tfin is final time, % c will be a matrix with three columns, one for each
species. % each row will be for another time point. % initial
conditions (c0) is a vector of three initial conds.
c0(1) = ca0 ;
c0(2) = 0 ;
c0(3) = 0;

[t,c] = ode45('exampleode',0,tfin,c0); % concentrations of the three species as function of time, each is a %

column vector
ca = c(:,1)
cd = c(:,2)
cu = c(:,3) % like to know the size of the result matrix
c last = size(c) % extract the final value of the species out of c
caf = c(last(1),1)
cdf = c(last(1),2)
cuf = c(last(1),3)
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Comprehensive Examination (EC-3 Make-up)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Open Book
Weightage : 40% No. of Pages = 2
Duration : 3 Hours No. of Questions = 8
Date of Exam : 19/11/2017(FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

Q.1. Evaluate the following expressions in MATLAB [5]

(a) (4\12)+4 (b) sign(-11/5)+6
(b) 3+4^2 (d) fix(-13/5)+5
(c) A=pi+0.1^2

Q.2. Write a MATLAB Program to evaluate the function: [5]

= ln 1−1

For any user specified value of x where x is a number < 1.0 . Use an if structure to
verify the value passed to the program is legal. If the value of x is legal, calculate
y(x), if not write suitable error message and quit.
Q.3. Create a matrix variable, mymat, which stores the following [5]

mymat= [2 5 8] and transform each of them into the following forms .

7 5 3

2 7
(i) 5 5 (iii) 8 5 2
8 3 3 5 7

(ii) 2 5 (iv) 2 5 8 2 5 8
7 8 7 5 3 7 5 3
5 3

Q.4. The number of bacteria Ns measured at different times t is given in the following table.
Determine an exponential function in the form [5]

NB= N e αt that best fits the data. Use the equation to estimate the number of bacteria
after 4.5 hr. Make a plot of the points and the equation.

t (hr) 1 2 3 4 5 6
NB 2000 4500 7500 15000 31000 64000

Q.5 A production facility is producing some nails that are supposed to have a diameter of
0.15 inch. At five different times, 10 sample nails were measured; their diameters were
stored in a file that has five lines and 10 diameters on each. First, create a data file to
simulate this data. Then, write a script to print the mean and standard deviation for each
of the five sets of sample nails. [5]
TA ZC164 (EC-3 Make-up) First Semester 2017-2018 Page 2

Q.6 Consider the chemical reaction equation [5]

where the initial values are y1(0) = 1; y2(0) = y3(0) = 0. This equation can be
regarded as stiff equation. Solve the problem with ode45.

Q.7 The following list gives the measured gas mileage in miles per gallon for 22 cars of the
same model. Plot the absolute frequency histogram and the relative frequency
histogram. [5]

23 25 26 25 27 25 24 22 23 25 26
26 24 24 22 25 26 24 24 24 27 23

Q.8. Tables of materials properties list density, in units of kg/m3, when the international
system of units (SI) is used and list specific weight, in units of lb/in3, when the U.S.
customary system of units are used. [5]

Write a user-defined MATLAB function that converts density to specific weight. For
the function name and arguments, use [sw] = DenToSw (den).

The input argument den is the density of a material in kg/m3, and the output argument
sw is the specific weight in lb/in3. Use the function in the Command Window to:

(a) Determine the specific weight of steel whose density is 7860 kg/m3
(b) Determine the specific weight of titanium whose density is 4730 kg/m3

*********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2016-2017

Comprehensive Examination
(EC-3 Makeup)

Course No. : TA ZC164

Q.1 Short Questions [1 x 5 = 5]

i. Complete the code to integrate the function, y(x) with respect to x.

syms x
y=x^2+6*x+5
a. int(y,x) b. int(x,y)
c. integrate(y,x) d. trapz(x,y)

ii. Which MATLAB command generates 50 uniform random numbers between -2 and
+2?
a. 10*rand(50)-2 c. 4*randn(50,1)-2
b. 2*rand(50,1)-2 d. 2*randn(50)-2

iii. Write the MATLAB commands to find a determinant for a 3*3 matrix
iv Which MATLAB command to find the covariance of x and y in MATLAB?
v. What sign will be prefix for a function if you want to find a differential equation or a
integral for any function

Q.2 Determine the output of the following [2]

(i) >> grades = [85 70 100 95 80 91];
>> grade_index = [3 5 1 6 4 1];
>> grades(grade_index)

(ii) >>a=eye(3);
>>sparse(a)
TA ZC164 (EC-3 Make-up) Second Semester 2016-2017 Page 2

Q.3 Examine the following (i) and (ii) for loops and determine the value of b at the end of each
loop and also the number of times each loop executes [4]

(i) (ii)
b=0 b=0
for i:-10:10 for i:10:-2:4
b=b+1 for j=2:2:i
end if j==6
break
end
b=b+j
end
end

Q.5 Create a mesh, surface plot and contour plot of the function
z= + for the interval -1≤ ≤ 1− 2 ≤ ≤ 2 . In each case plot the real part

of z versus x and y [4]

Q.6 The open loop transfer function G(s)H(s) of a control system is [4]
( )( )=
( + 0.5)( 2 + 0.5 + 8)
Plot the root loci for the system using MATLAB

Q.7 For a project, some biomedical engineering students are designing a device that will
monitor a person’s heart rate while on a treadmill. The device will let the subject know
when the target heart rate has been reached. A simple calculation of the target heart rate
(THR) for a moderately active person is THR=(220−A)*0.6 where A is the person’s age.
Write a function that will calculate and return the THR [4]
TA ZC164 (EC-3 Make-up) Second Semester 2016-2017 Page 3

Q.8Using this polynomial, predict the accumulation over the range of 20 - 30 hours. Finally,
calculate the time derivative of the accumulation over the period 0-10 hours. [4]
Mass of A accumulated as a function of time:

Mass of A accumulated (arbitrary units) Time (hours)

9 1
55 3
141 5
267 7
345 8
531 10

Only write commands and write what we can get as far as analysis

Q.9 Consider the system of reactions in a constant volume, constant temperature batch reactor.
A —- > D reaction 1
A + A —— > U reaction 2

where D is a desired product, U is undesired product and where k1 and k2 are the rate
constants for reactions 1 and 2. Let k1, k2 be given parameters, and the initial
concentration of A (ca0) be a design variable. The independent variable is time (tfin), and
the dependent variables are the concentration of the species, ca, cd, and cu. Note that in
most design situations k1 and k2 might be design parameters, adjusted via the
temperature.

By writing the mass balance equations over the batch reactor, the system of differential
equations is the following: they are not linearly independent.
d(ca)/dt = -k1(ca) - k2(ca)2
d(cd)/dt = k1 (ca)
d(cu)/dt = k2 (ca)2
Solving this problem in Matlab involves two parts. First, write a function file that
describes the set of ODEs in terms of a single, combined matrix variable (the dependent
variable). Plot the concentration variables with respect to time

Next, in a second, (main) script invoke the ode solver, ode45.

This script that might include other things like the initial conditions and other given
parameters. (The solution to a single ODE i s analogous, but the dependent variable is not
a matrix). [3 Marks- Algorithms 3 Marks- Code 3 Marks- Functions]
*********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Mid-Semester Test (EC-2 Make-up)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 30%
No. of Pages= 2
Duration : 2 Hours
No. of Questions = 7
Date of Exam : 18/03/2018(FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.

Q.1. Short Questions [2 x 5 = 10]

(a) Create a 2 × 4 matrix of random integers, each in the range from 10 to 30

Solution:
A=randi([10,30],2,4)

(b) If mat = [4 3 1; 2 5 6] what is the output of mat = [2:4; 3:5]

Solution:

mat =

2 3 4
3 4 5

Solution:

linspace(-3,-15,5)

(d) Create a 4 × 2 matrix of all zeros and store it in a variable. Then, replace the second
row in the matrix with a 3 and a 6.
Solution:

z=zeros(4,2)
z(2,:)=[6 3]

(e) Write a command to find the Eigen value of a 2*2 matrix z.

Solution:

Z=[1,2;3,4]
a=eig(Z)

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 1 of 4

Q.2. Evaluate the value of c if c=a.\b , a and b are matrices shown below [2]

3 1 1 3

a    b   

1 2  2
2

C= -1/3 3
-2 1

Q.3. What is the output of the following output statement [2]

ii) fprintf(‘%3.2f\n’, 1234.5678)

solution:

The output is 1234.57

ii) for i=0:10:100

if mod(i,10) ==0
disp(i)
else
continue;
end
end

Solution:

The output will be

0 10 20 30 40 50 60 70 80 90 100

Use MATLAB to find (i) compression x in each spring and (ii) Potential energy stored in
each spring [4]

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2 of 4

Solution

>> F=[11,7,8,10,9];
>> k=[1000,600,900,1300,700];
>> x=F./k;
>> e=0.5*k.*x.^2;

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2

Solution

a=v.^2
b=v.^3
c=v.^v
d=v.*(1/2)

Solution:

% Script to Classify Clouds

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 3 of 4

x=input('Enter the height of the cloud in feet')
if (x>=0 & x<6500)
disp('Low cloud')
elseif (x>=6500 & x<20000)
disp('Medium cloud')
elseif (x>20000)
disp('high cloud')
end

% Create the file calcaverag(a)

function avg= calcaverag(a)

n=length(a);
sum=0;
for i=1:n
sum=sum+i;
avg=sum/n;
end

********

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 4 of 4

Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Mid-Semester Test (EC-2 Makeup)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 30%
No. of Pages= 3
Duration : 2 Hours
No. of Questions = 10
Date of Exam : 15/10/2017(FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

Marks Short questions [2 x 6= 12]

Determine the max value and min value of each row and each column of Matrix A using
MATLAB. [2]

5 3 −3

A=[−2 2 3]

2 5 7

Q.2. Write the commands to solve the following and obtain Laplace transform respective for the
following equations
(a) t2e-t (b) 3 e -2t sin 4t [2]

Q.3. Create a variable, pounds, to store a weight in pounds. Convert this

to kilograms and assign the result to a variable kilos. The conversion factor
is 1 kilogram = 2.2 pounds. [2]

Q.4. Solve the following system of three linear equations using MATLAB [2]

3x+ 2y+5z=7.5
-4.5x+2y+ 3z=5.5
-5x+ y-2.5z=4.5

Q.5. Examine the following for loops and determine the value of ires at the end of each of the
loops, and also the number of times each loop executes [2]

ires=0;
for index1=10:-2:4
for index2=2:2:index1
if index2==6
break
end
ires=ires+index2
end
end
TA ZC164 (EC-2 Make-up) First Semester 2017-2018 Page 2

Q.6. Examine the following for loop and determine the value of b at the end of each of the loops, and
also the number of times each loop executes [2]

b=0;
for i:-10:10
b=b+1;
end

Q.7. Write a function sumsteps2 that calculates and returns the sum of 1 to n in steps of 2, where n is an
argument passed to the function. For example, if 11 is passed, it will return 1 + 3 + 5 + 7 + 9 + 11.
Do this using a for loop. Calling the function will look like this:

>> sumsteps2(11) [4]

ans =
36

Q.8. The steady-state heat conduction q from a cylindrical solid wall is determined by:
=2 ( 1 − 2)
1
( 2)

Where k is the thermal conductivity. 401 watts/0C/m of length L=300 cm with outer radius
r2= 5 cm and inner radius r1=3 cm. The external temperature is T2=20 0 C and the internal
temperature is T1=100 0C . Perform MATLAB command and calculate. [4]

Q.9. Create three row vectors [4]

a=[ 7 2 -3 1 0]
b=[-3 10 0 7 -2]
and c = [1 0 4 -6 5]

(a) Use the three vectors in a MATLAB command to create a 3 x 5 matrix in which the
rows are the vectors a,b and c
(b) Use the three vectors in a MATLAB command to create a 5 x3 matrix in which
the columns are the vectors a,b, and c
TA ZC164 (EC-2 Make-up) First Semester 2017-2018 Page 3

Q.10. The electrical circuit shown consists of resistors and voltage sources. Determine the current in each
resistor, using the mesh current method based on Kirchhoff’s voltage law. [6]

V1= 20 V, V2= 12 V, V3=40 V

R1= 18 Ω R2= 10 Ω R3= 16 Ω R4= 6 Ω R5= 15 Ω R6= 8 Ω R7= 12 Ω R8= 14 Ω

Kirchhoff’s voltage law states that the sum of the voltage around a closed circuit is zero. In the
mesh current method a current is first assigned for each mesh (i1, i2, i3, i4 in the figure). Then
Kirchhoff’s voltage law is applied for each mesh. This results in a system of linear equations for
the currents (in this solution gives the values of the mesh currents. The he current in a resistor that
belongs to two meshes is the sum of the currents in the corresponding meshes. It is convenient to
assume that all the currents are in the same direction (clockwise in this case). In the equation for
each mesh, the voltage source is positive if the current flows to the – pole, and the voltage of a
resistor is negative for current in the direction of the mesh current. For example the equation for
two meshes is shown. Write equations for other two meshes, arrange it in matrix form and write
the script to solve it for the current value

***********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Mid-Semester Test (EC-2 Makeup)

Course No. : TA ZC164

Marks Short questions [2 x 6= 12]

Determine the max value and min value of each row and each column of Matrix A using
MATLAB. [2]
5 3 −3 A=[−2 2 3 ]
2 5 7

Solutions

randi([10,30],2,4)

Q.2. Write the commands to solve the following and obtain Laplace transform respective for the
following equations
(a) t2e-t (b) 3 e -2t sin 4t [2]

(a) f=t^2*exp(-t)
laplace(f)

(b) f=3*exp(-2t)*sin(4*t)
laplace(f)

Q.3. Create a variable, pounds, to store a weight in pounds. Convert this to kilograms
and assign the result to a variable kilos. The conversion factor is 1 kilogram = 2.2 pounds.
[2]

>>pounds=input(‘weight in pound’);
>>kilos=2.2*pounds;
>>disp(kilos) 2 mar k
Q.4. Solve the following system of three linear equations using MATLAB [2]

3x+ 2y+5z=7.5
-4.5x+2y+ 3z=5.5
-5x+ y-2.5z=4.5

>> a=[3,2,5;-4.5,2,3;-5,1,-2.5]
2 mark
>> b=[7.5;5.5;4.5]
>> X=inv(a)* b

Q.5. Examine the following for loops and determine the value of ires at the end of each of the
loops, and also the number of times each loop executes [2]

ires=0;
for index1=10:-2:4
for index2=2:2:index1
if index2==6
break
end
ires=ires+index2
end
end

Initial value=0 then it will get into loop

Index1=10:-2:4 index1= 10
Then it will go to inner loop
2 mark
Index2 will be aassign to 2 index2=2
Condition is false inner loop
Ires=0+2 2,6,8, 12,14,18,20,24
The loop will terminate after executing 8 times
TA ZC164 (EC-2 Make-up) First Semester 2017-2018 Page 2

Q.6. Examine the following for loop and determine the value of b at the end of each of the loops, and
also the number of times each loop executes [2]

b=0;
for i:-10:10
b=b+1;
end

>> sumsteps2(11)
[4]
ans =
36

Q.8. The steady-state heat conduction q from a cylindrical solid wall is determined by:
=2 ( 1 − 2)
1
( 2)

Q.9. Create three row vectors

[4]

a=[ 7 2 -3 1 0]
b=[-3 10 0 7 -2]
and c = [1 0 4 -6 5]
(a) Use the three vectors in a MATLAB command to create a 3 x 5 matrix in which
the rows are the vectors a,b and c
(b) Use the three vectors in a MATLAB command to create a 5 x3 matrix in which
the columns are the vectors a,b, and c

TA ZC164 (EC-2 Make-up) First Semester 2017-2018 Page 3

Q.10. The electrical circuit shown consists of resistors and voltage sources. Determine the current in each
resistor, using the mesh current method based on Kirchhoff’s voltage law. [6]

V1= 20 V, V2= 12 V, V3=40 V

R1= 18 Ω R2= 10 Ω R3= 16 Ω R4= 6 Ω R5= 15 Ω R6= 8 Ω R7= 12 Ω R8= 14 Ω

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Open Book No. of Pages = 3
Weightage : 40%
Duration : 3 Hours No. of Questions = 8
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the
answer book.
2. All parts of a question should be answered consecutively. Each answer should start
from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

All Question Carries 5 Marks – 8*5=40 Marks

Q.1 Evaluate the following expressions in MATLAB

(a) (4\12)+4 (b) sign(-11/5)+6
(b) 3+4^2 (d) fix(-13/5)+5
(c) A=pi+0.1^2

Q.2 Write a MATLAB Program to evaluate the function:

= ln 1−1

For any user specified value of x where is a number < 1.0 . Use an if structure to
verify the value passed to the program is legal. If the value of x is legal, calculate
y(x), if not write suitable error message and quit.

Q.3 Create a matrix variable, mymat, which stores the following

mymat= [2 5 8] and transform each of them into the following forms .

7 5 3
2 7
(i) 5 5
8 3
(iii) 8 5 2
3 5 7
(ii) 2 5
7 8
5 3 (iv) 2 5 8 2 5 8
7 5 3 7 5 3

Q.4 The number of bacteria Ns measured at different times t is given in the

following table. Determine an exponential function in the form

NB= N e αt that best fits the data. Use the equation to estimate the number of
bacteria after 4.5 hr. Make a plot of the points and the equation.

t (hr) 1 2 3 4 5 6
NB 2000 4500 7500 15000 31000 64000

Q.5 A production facility is producing some nails that are supposed to have a
diameter of 0.15 inch. At five different times, 10 sample nails were measured;
their diameters were stored in a file that has five lines and 10 diameters on
each. First, create a data file to simulate this data. Then, write a script to print
the mean and standard deviation for each of the five sets of sample nails.

Q.6 Consider the chemical reaction equation

where the initial values are y1(0) = 1; y2(0) = y3(0) = 0. This equation can
be regarded as stiff equation. Solve the problem with ode45.

Q.7 The following list gives the measured gas mileage in miles per gallon for 22
cars of the same model. Plot the absolute frequency histogram and the relative
frequency histogram.

23 25 26 25 27 25 24 22 23 25 26
26 24 24 22 25 26 24 24 24 27 23
8. Tables of materials properties list density, in units of kg/m3, when the international
system of units (SI) is used and list specific weight, in units of lb/in3, when the
U.S. customary system of units are used.

Write a user-defmed MATLAB function that converts density to specific weight. For
the function name and arguments, use [sw] = DenToSw (den).

The input argument den is the density of a material in kg/m3, and the output argument
sw is the specific weight in lb/in3. Use the function in the Command Window to:

(a) Determine the specific weight of steel whose density is 7860 kg/m3
(b) Determine the specific weight of titanium whose density is 4730 kg/m3
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Comprehensive Examination Test (EC-3 Regular)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book No. of Pages = 3
Weightage : 40%
Duration : 3 Hours No. of Questions = 8
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the
answer book.
2. All parts of a question should be answered consecutively. Each answer should start
from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

All Question Carries 5 Marks – 8*5=40 Marks

Q.1 Evaluate the following expressions in MATLAB

(a) (11/5)+6 (b) round(-11/5)+6
(b) (3^2)^3 (d) ceil(-11/5)+6
(c) 11/(5+6)

Solutions
(a) 8.2
(b) 4 Round to nearest decimal or integer
(c) 729
(d) 4 Round toward positive infinity
(e) 1

Q.2 Using the zeros, ones, and eye commands create the following arrays
0 0 1 1

a) [0 0 1 1]

0 0 0 0

1 1 1 1
Solution (a)

Step 1 : Create a=zeros(4) will create a 4*4 matrix of 0

Step2 : a(:,3:4)=ones will replace contents of 3rd and 4th Column by ones
Step3: a(4,:)=ones will replace contents of 4th row with ones
Step4: a(3,:)=zeros will replace the 3rd row by zero

1 1 0 0 1
b) [1 1 0 0
0]

1 1 0 0 0

Solution (b)

Step1 : z=ones(4) will create 4*4 matrix

Step2: z(:,3:5)=zeros will replace elements of 3 to 5th column with zeros
Step3: z(1,5)=ones will replace 5th column of 1st row by ones

How many gallons of each blend can be manufactured if 7,600 gallons of orange
juice, 4,900 gallons of pineapple juice, and 3,500 gallons mango juice are available?
Write a system of linear equations and solve in MATLAB only steps

Solution

Consider orange=x , pineapple=y, mango= z vectors

1 gallon Orange Blend: X

1 gallon Pineapple Blend: Y
1 gallon mango Blend: Z

Set up Equations
X=3*x+0.75*y+0.25*z
Y=x+2.5*y+0.5*z
Z=0.5*x+0.5*y+3*z

Now Put X=7600, Y=4900,Z=3500

MATLAB

a=[3,0.75,0.25;1,2.5,0.5;0.5,0.5,3]
b=[7600;4900;3500]

c=inv(a)* b

/* Answers not require only steps give full marks*/

Ans [ x= 2246 gallons orange juice, y=933 gallons pineapple, z=636 gallons mango juice]

Q.4 Viscosity, μ, is a property of gases and fluids that characterizes their resistance to
flow. For most materials viscosity is highly sensitive to temperature. Below is a table
that gives the viscosity of SAE 10W oil at different temperatures (Data from B.R.
Munson, D.F. Young, and T.H. Okiishi, Fundamentals of Fluid Mechanics, 4th ed.,
John Wiley and Sons, 2002). Determine an equation that can be fitted to data.

T (0C) -20 0 20 40 60 80 100 120

µ (N 4 0.38 0.095 0.032 0.015 0.0078 0.0045 0.0032
s/m2)
x 10-5

The straight line does not fit the point , Therefore the equation can be fitted as
ln( ) = + 2+.
2 1
0

Hint:Use semilog function

Step1 : Create the Data in Excel file

Step2 : Say the file name as testexam.xlsx
Step3: The data in file is as follows

name Age marks

Samir 40 54
Ramesh 36 67
ravi 47 72

Step4: Opening in matlab

s=uiimport('C:\Users\BITS7\Documents\MATLAB\testexam.xlsx')

both test data and numeric data will be loaded

Q.6. Use ODE45 function to solve the differential equation

2′ =+3 y(1)= -2

over the interval 1< t <4

>> ode

Q.7 Plot the function = − for x between 0 and 2 in steps of 0 1 . Create the following plot types (a) stem plot (b) stair plot (c) bar plot . Be sure to
include titles and axis labels on all plots

x=0:0.1:2;
y=exp(-x)*sin(x);
stem(x,y);
bar(x,y);
title(‘stem and bar Plot’);
xaxis(‘x coordinate’)
yaxis(‘function’);

Q.8 Fluid flows in pipe networks can be analyzed in a manner similar to that used
for electric resistance networks. The Figure below shows a network with three
pipes. The volume flow rates in the pipes are q1, q2, and q3 . The pressures at
the pipe ends are pa, pb and pc, . The pressures at the junction is p1. Under certain
condition, the pressure flow rate relation in a pipe has the same form as the
voltage-current relation in a resistor.Thus, for the three pipes, we have

Where the Ri are the pipe resistances. From conservation of

mass q1= q2+ q3

a. Set up these equations in a matrix form Ax=b suitable for solving for the
three flow rates q1, q2, q3, and the pressure p1 given the values of pressures
pa, pb and pc and the values of resistances R1, R2, and R3. Find
expressions of A and B
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018
Comprehensive Test (EC-3 Regular)

Course No. : POWTP ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 50%
No. of Pages= 3
Duration : 3 Hours
No. of Questions = 10
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

Q.1. Calculate the efficiency of a power plant using function in MATLAB if the efficiencies of
the boiler, turbine and generator are 88, 40 and 98%, respectively.

Write the function as general so that for varying efficiency it should calculate the overall
efficiency. Formula to calculate the efficiency is

Epower plant = [Eboiler] [Eturbine] [Egenerator] [5]

Q.2 The relationship between two variables P and t is known to be:

= +

The following data points are given Determine the constants m and b by curve-fitting the
equation to the data points. Make a plot of P versus t. In the plot show the data points with
markers and the curve-fitted equation with a solid line. (The curve fitting can be done by
writing the reciprocal of the equation and using a first-order polynomial.) [5]

t 1 3 4 7 8 10
P 2.1 4.6 5.4 6.1 6.4 6.6

Q.3 Create a mesh, surface plot and contour plot of the function z= + for the interval -1≤

≤ 1− 2 ≤ ≤ 2 . In each case plot the real part of z versus x and y [5]

Q.4 Solve the equation using MATLAB (only program no calculations) using Gauss
technique [5]
2 +3 + =9
+2 +3 =63 ++2 =8
Q.5 The xlswritefunction can write the contents of a cell array to a spreadsheet.A manufacturer
stores information on the weights of some parts in a cell array. Each row stores the part
identifier code followed by weights of some sample parts.To simulate this, create the
following cell array:

>>parts = {‘A22’, 4.41 4.44 4.39 4.39‘Z29’, 8.88 8.95 8.84 8.92}

Then, write this to a spreadsheet file [5]

Q.6 Use ODE45 function to solve the differential equation

2 ′= +3

y(1)= -2

over the interval 1<t <4 [5]

Q.7 Using the zeros, ones, and eye commands create the following arrays

[5]
Q.8 The temperature dependence of vapor pressure p can be estimated by the
Antoine equation:

where ln is the natural logarithm, p is in mm Hg, T is in kelvins, and A, B, and C are material
constants. For toluene (C6H5CH3) in the temperature range
from 280 to 410 K the material constants are A= 16.0137 B = 3096.52, , and C = –53.67
Calculate the vapor pressure of toluene at 315 and 405 K in MATLAB [5]

Q.9 Write a script called prtemps that will prompt the user for a maximum Celsius value in the
range from –16 to 20; error-check to make sure it’s in that range. Then, print a table
showing degrees F and degrees C until this maximum is reached. The first value that
exceeds the maximum should not be printed. The table should start at 0 degrees F, and
increment by 5 degrees F until the max (in DEG C) is reached. Both temperatures should
be printed with a field width of 6 and one decimal place. The formula is C = 5/9* (F – 32).
For example, the execution of the script might look like this (the format should be exactly
like this): [5]
>> prtemps

When prompted, enter a temp in degrees C in range −16

to 20.
Enter a maximum temp: 30
Error! Enter a maximum temp: 9

F C

0.0 −17.8
5.0 −15.0
10.0 −12.2
15.0 −9.4
20.0 −6.7
25.0 −3.9
30.0 −1.1
35.0 1.7
40.0 4.4
45.0 7.2

Q.10 Write a function that will receive a name and department as separate strings and will
create and return a code consisting of the first two letters of the name and the last two
letters of the department. The code should be uppercase letters. For example,
>> namedept(‘Robert’,‘Mechanical’) [5]
ans =
ROAL

***********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
First Semester 2017-2018

Comprehensive Examination
(EC-3 Makeup)

Course No. : POW TP ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 50%
No. of Pages = 3
Duration : 3 Hours
No. of Questions = 9
Date of Exam : (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. Write the codes only, solve only if the question requires it

Q.1 Short Questions [1 x 6 = 6]

i. Complete the code to integrate the function, y(x) with respect to x.

syms x
y=x^2+6*x+5
a. int(y,x) b. int(x,y)
c. integrate(y,x) d. trapz(x,y)

ii. Which MATLAB command generates 50 uniform random numbers between -2 and
+2?
a. 10*rand(50)-2 c. 4*randn(50,1)-2
b. 2*rand(50,1)-2 d. 2*randn(50)-2

Q.2 Determine the output of the following [3]

(i) >> grades = [85 70 100 95 80 91];
>> grade_index = [3 5 1 6 4 1];
>> grades(grade_index)

(ii) >>a=eye(3);
>>sparse(a)

POWTP ZC164 (EC-3 Make-up) First Semester 2017-2018 Page 2

Q.3 Examine the following (i) and (ii) for loops and determine the value of b at the end of each
loop and also the number of times each loop executes [6]

(i) (ii)
b=0 b=0
for i:-10:10 for i:10:-2:4
b=b+1 for j=2:2:i
end if j==6
break
end
b=b+j
end
end

Q.5 Create a mesh, surface plot and contour plot of the function
z= + for the interval -1≤ ≤ 1− 2 ≤ ≤ 2 . In each case plot the real part of

z versus x and y [5]

Q.6 The open loop transfer function G(s)H(s) of a control system is [5]
( )( )=
( + 0.5)( 2 + 0.5 + 8)
Plot the root loci for the system using MATLAB

Q.7 For a project, some biomedical engineering students are designing a device that will monitor a
person’s heart rate while on a treadmill. The device will let the subject know when the target
heart rate has been reached. A simple calculation of the target heart rate (THR) for a
moderately active person is THR=(220−A)*0.6 where A is the person’s age. Write a
function that will calculate and return the THR [5]
POW TP ZC164 (EC-3 Make-up) First Semester 2017-2018 Page 3

Q.8 Using this polynomial, predict the accumulation over the range of 20 - 30 hours. Finally,
calculate the time derivative of the accumulation over the period 0-10 hours. [5]
Mass of A accumulated as a function of time:

Mass of A accumulated (arbitrary units) Time (hours)

9 1
55 3
141 5
267 7
345 8
531 10

Only write commands and write what we can get as far as analysis

Q.9 Consider the system of reactions in a constant volume, constant temperature batch reactor.
A —- > D reaction 1
A + A —— > U reaction 2

Next, in a second, (main) script invoke the ode solver, ode45.

This script that might include other things like the initial conditions and other given
parameters. (The solution to a single ODE i s analogous, but the dependent variable is not a
matrix). [10]
*********
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Comprehensive Test (EC-2 Regular)

Course No. : MTU ZC221

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 30%
No. of Pages= 2
Duration : 2 Hours
No. of Questions = 8
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

1. Using the zeros, ones, and eye commands create the following arrays [3]

Solutions:

(a)

>> a=zeros(2,4)
a=
0 0 0 0
0 0 0 0
>> a(:,4:5)=ones
a=
0 0 0 1 1
0 0 0 1 1
>> a(1,1)=ones
a=
1 0 0 1 1
0 0 0 1 1
>> a(2,2)=ones
a=
1 0 0 1 1
0 1 0 1 1
(b)

>> a=zeros(4)

0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

>> a(:,3:4)=ones

0 0 1 1
0 0 1 1
0 0 1 1
0 0 1 1

>> a(4,:)=ones

0 0 1 1
0 0 1 1
0 0 1 1
1 1 1 1

>> a(3,:)=zeros

0 0 1 1
0 0 1 1
0 0 0 0
1 1 1 1
(c) >> c=ones(4)
c=
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
>> c(:,3:5)=zeros
c=
1 1 0 0 0
1 1 0 0 0
1 1 0 0 0
1 1 0 0 0
>> c(1,5)=ones
c=
1 1 0 0 1
1 1 0 0 0
1 1 0 0 0
1 1 0 0 0

2. Evaluate the following expressions in MATLAB [3]

(a) (6/12)+4
6 will be divided by 12 which gives 0.5 and then added to 4 = 4.5

(b) sign(-11/5)+6

-11/5 gives -2.2 however sign will give +1 or –1ve only therefore -1+6=5
(c) fix(-13/5)+5
-13/5 gives -2.6 which truncated decimal point so -2+5=3

3. Tables of materials properties list density, in units of kg/m3, when the international
system of units (SI) is used and list specific weight, in units of lb/in3, when the U.S.
customary system of units are used.

Write a user-defined MATLAB function that converts density to specific weight. For the
function name and arguments, use [sw] = DenToSw (den).

The input argument den is the density of a material in kg/m3, and the output argument sw
is the specific weight in lb/in3. Use the function in the Command Window to:
(a) Determine the specific weight of steel whose density is 7860 kg/m3
(b) Determine the specific weight of titanium whose density is 4730 kg/m3 [4]

clear
clc
den=7860;
sw=dentosw(den);
disp('The specific weight of steel in lb/inch3
is'); disp(sw);

den=4730;
sw=dentosw(den);
disp('The specific weight of titanium in lb/inch3
is'); disp(sw);

% function file

function sw = dentosw(den)
sw=den*3.605/100000;
end

4. Write a MATLAB Program to evaluate the function: [4]

= ln 1−1

For any user specified value of x where x is a number < 1.0 . Use an if structure to verify
the value passed to the program is legal. If the value of x is legal, calculate y(x), if not
write suitable error message and quit.

x=input('The value of x');

if(x>=1.0)
disp('Error Input');
else
y=log(1.0/(1.0-x));
disp(y);
end
5. The following polynomials are given: [4]

p1(x) = x5 + 2x4 – 3x3 + 7x2 – 8x + 7

Use MATLAB functions with polynomial coefficient vectors to evaluate the expressions
at x = 2.

P=[1,2,-3,7,-8,7]
Polyval(P,2)

>> P=@(X) X^5+2X^4-3X^3+7*X^2-

8*X+7 P =
function_handle with value:
@(X)X^5+2*X^4-3*X^3+7*X^2-8*X+7
>> P(2)
ans =
59

6. Find the Laplace transforms of the following [4]

 1 
 
a)  2 
s (s + 4) 

dy
b) dt + 3y = e-2t y(0) = 2

>> syms s
>> f=(1/s*(s^2+4));
>> ilpalace(f)

>>syms s
>>f=(2s+5)/(s+2)*(s+4);
>>ilapace(f)
7. The following functions describe the oscillations in electric circuits and the vibrations of
machines and structures. Plot these functions on the same plot. Because they are similar,
decide how best to plot and label them to avoid confusion [4]
( ) = . − . sin(3t+2)
( ) = . − . cos (5t-3)

8. Write a function called geomser that will receive values of r and n, and will calculate and
return the sum of the geometric series: 1 + r + r2 + r3 + r4 + ... + rn [4]

The following examples of calls to this function illustrate what the result should be:
>> geomser(1,5)
ans = 6
>> disp(geomser(2,4))
31

function s=geomser(r,n)
v = [1,r*ones(1,n)];
p = cumprod(v);
s = sum(p);
end
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Comprehensive Test (EC-2 Regular)

Course No. : MTU ZC221

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 50%
No. of Pages= 2
Duration : 3 Hours
No. of Questions = 10
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

1. Evaluate the output of the following MATLAB expressions, if a is the input vector containing the
following elements a=[1,2,3;4,5,6]

= 1020√ (1 − 0.01√ ) Where P is population in thousands. Set up vector P that starts at 10 and increments by 10 upto 200. Use
element by element computations to determine Q for each population in P

3. Using commands in MATLAB to plot a 3-D figure of the trajectory of a particle defined by the
following set of equations for x,y, and z.Use commands for proper labelling and title

x=t
y=t*cos(t)
z= e 0.2t

a. Determine the pressure at 27 deg C using Interp1 (Refer Q.4)

5. The operations manager of a plant that manufactures tires wants to compare the actual inner
diameters of two grades of tires, each of which is expected to be 575 millimeters. A sample of
five tires of each grade was selected, and the results representing the inner diameters of the
tires, ranked from smallest to largest, are as follows:

Grade X 568 570 575 578 584

Grade Y 573 574 575 577 578

Create a data file tire.xls in Excel and calculate the mean, median and std. dev in MATLAB and
export it to excel or write the computed values to excel

6. Construct an equivalent statement using single if statement in MATLAB for the following block
(a) If x<y
If z<10
W=x*y*z

(b) Accepts marks of the students as p and if p is greater than 50 print pass else print fail

7. Given the system of equations

3x1  7x2  13x3  76

x1  5x2  3x3  28
12x1  3x2 - 5x3  1
find the solution using Gauss-Seidal method. Use x1 , x2 , x3  1 0 1 as the initial guess.
8. Consider the example = 3 − with initial condition x(0)=0, Solve the following equation using ODE45 in MATLAB

Course No. : MTU ZC221

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 50%
No. of Pages= 2
Duration : 3 Hours
No. of Questions = 10
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

1. Evaluate the output of the following MATLAB expressions, if a is the input vector containing the
following elements a=[1,2,3;4,5,6]

Solutions:

1a) a=[1,2,3;4,5,6]
The vector is 2 rows and three colums already so no command needed
1b) a(2,:)*2 for row and a (:,2)*2 for column

1c) a(3,:)=mean(a)

1d) a(:,4)=mean(a,2) also [ ] is ok in between

1e) min(a) for column and min(a,[],2) for row
2. The demand for water during a fire is often most important factor in the design of distribution
storage tanks and pumps. For communities with population less than 1,00,000 the demaond Q
in gallons/min can be calculated by

Solution

clear,
clc
P=10:10:200;
Q=1020*sqrt(P).*(1-.01*sqrt(P))

3. Using commands in MATLAB to plot a 3-D figure of the trajectory of a particle defined by the
following set of equations for x,y, and z.Use commands for proper labelling and title

x=t
y=t*cos(t)
z= e 0.2t

Solution
clear;
clc;
t=0:0.1:10*pi;
x=t;
y=t.*cos(t);
z=exp(0.2*t)
;
plot3(x,y,z)
title('Plot 3 Dimensional
figure'); xlabel('x axis');
ylabel('y axis');
zlabel('z axis');
grid on
4. Consider the following steam table and determine the values using polyfit, polyval and interp1
commands using MATLAB

T, 0.01 5 10 15 20 25 30 35 40 45 50 55
deg C
P, bar 0.0061 0.0087 0.0123 0.0170 0.0234 0.0317 0.0425 0.0563 0.0738 0.0959 0.1.235 0.1576

a. Determine the pressure at 27 deg C using Interp1 (Refer Q.4)

Solution

T=[0.01,5:5:55]
P=[0.061,0.087,0.0123,0.0170,0.0234,0.0317,0.0425,0.0563,0.0738,0.0959,0.1235,0.1576]

(a) q=interp1(T,P,27)
(b) q=interp1(P,T,0.0465)
(c) q=polyfit(T,P,2) m=polyval(q,70)
(d) q=polyfit(P,T,2) m=polval(q,0.2)
(e) TK=T+273 m=polyfit(TK,P,2)

5. The operations manager of a plant that manufactures tires wants to compare the actual inner
diameters of two grades of tires, each of which is expected to be 575 millimeters. A sample of
five tires of each grade was selected, and the results representing the inner diameters of the
tires, ranked from smallest to largest, are as follows:

Grade X 568 570 575 578 584

Grade Y 573 574 575 577 578

Create a data file tire.xls in Excel and calculate the mean, median and std. dev in MATLAB and
export it to excel or write the computed values to excel
Solution:

Create the following excel file in MATLAB folder

[s,x]=xlsread('Tire.xlsx');
s(:,6)=mean(s,2);
k=s(1:2,1:5);
s(:,7)=median(k,2);
m=s(1:2,1:5);
s(:,8)=std(m,0,2);
xlswrite('Tire.xlsx',s);

6. Construct an equivalent statement using single if statement in MATLAB for the following block
(a) If x<y
If z<10
W=x*y*z

(a) if (x<y) & z<10 w=xyz; end

(b) Accepts marks of the students as p and if p is greater than 50 print pass else print fail

% Accepts marks and Print fail or pass

p=input(‘Enter marks of students’);

if p>50
disp('Pass')
elseif disp('Fail')
end

7. Given the system of equations

3x1  7x2  13x3  76

x1  5x2  3x3  28
12x1  3x2 - 5x3  1
find the solution using Gauss-Seidal method. Use x1 , x2 , x3  1 0 1 as the initial guess.
x=[1 0 1];
for i=1:5
xsave=x;
x(1)=(1/3)*(76-13*x(3)-7*x(2));
x(2)=(1/5)*(28-x(1)-3*x(3));
x(3)=(1/5)*(12*x(1)+3*x(2)-1);
if abs(max(xsave-x))< 1e-10
break;
end
end

fprintf('x = %2.4f %2.4f %2.4f\n',x)

fprintf('# of iterations = %g\n',i)

fprintf('x = %11.5f %11.5f %11.5f\n',x)

8. Consider the example = 3 − with initial condition x(0)=0, Solve the following equation using ODE45 in MATLAB

tspan = [0
5]; x0 = 0;
[t,y] = ode45(@(t,x) 3*exp(-t), tspan,
x0); figure(1);
plot(t,y,'-o');

function Efficiency=ottocycle(d,S,VC)
% d is bore diameter ENTER IN MM
% S is stroke ENTER IN MM
%VC is clearance volume ENTER IN MM3
VS=(pi*d^2*S/4)*(1/1000) ;% conversion to
cm disp(VS);
r=(VC+VS)/VC; % no difference in unit as units are same
disp(r);
Efficiency=(1- (1/(r^0.4)))*100;
end

Apart from MATLAB Significance one plant example with step wise methodology was desired
Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Comprehensive Test (EC-2 Regular)

Course No. : MTU ZC221

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 50%
No. of Pages= 2
Duration : 3 Hours
No. of Questions = 10
Date of Exam : (AN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.
4. No need for calculations, only commands and logics is desired wherever necessary.

1. Evaluate the output of the following MATLAB expressions, if a is the input vector containing the
following elements a=[1,2,3;4,5,6]

3. Using commands in MATLAB to plot a 3-D figure of the trajectory of a particle defined by the
following set of equations for x,y, and z.Use commands for proper labelling and title

x=t
y=t*cos(t)
z= e 0.2t

a. Determine the pressure at 27 deg C using Interp1 (Refer Q.4)

5. The operations manager of a plant that manufactures tires wants to compare the actual inner
diameters of two grades of tires, each of which is expected to be 575 millimeters. A sample of
five tires of each grade was selected, and the results representing the inner diameters of the
tires, ranked from smallest to largest, are as follows:

Grade X 568 570 575 578 584

Grade Y 573 574 575 577 578

Create a data file tire.xls in Excel and calculate the mean, median and std. dev in MATLAB and
export it to excel or write the computed values to excel

6. Construct an equivalent statement using single if statement in MATLAB for the following block
(a) If x<y
If z<10
W=x*y*z

(b) Accepts marks of the students as p and if p is greater than 50 print pass else print fail

7. Given the system of equations

3x1  7x2  13x3  76

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING (MATLAB)
Nature of Exam : Closed Book
Weightage : 30%
No. of Pages= 2
Duration : 2 Hours
No. of Questions = 5
Date of Exam : 29/09/2018 (FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.

Q.1. Short Questions [2 x 5 = 10]

(a) Evaluate the following expressions in MATLAB

i. (6/12)+4
ii. sign(-11/5)+6
iii. fix(-13/5)+5
iv. vec=9:-2:1

(b) If mat = [4 3 1; 2 5 6] what is the output of mat = [1:2,2:3]

(c) Using the linspace function, create the following vectors: –4 –8 –12 –16 –20
(d) Evaluate the value of c if c=a.\b , a and b are matrices shown below

3 1 1 3

a    b   
2
1
2  2
(e) What is the output of the following output statement
(i) fprintf(‘%3.2f\n’,456.789)
(ii) >> mat=[2:4;3:5]
>> mat(1,2)=11

Q.2 Evaluate the output of the following MATLAB expressions, if a is the input vector
containing the following elements a=[1,2,3;4,5,6] [ 5 Marks]

TA ZC164 (EC-2 Regular) First Semester 2018-2019 Page 1 of 2

Q.3 Examine the following (i) and (ii) for loops and determine the value b at the end of each
of loop and also the number of times each loop executes [ 5 Marks]

(i) (ii)

b=0 b=0
for i=-10:10 for i=10:-2:4
b=b+1 for j=2:2:i
end if j==6
break
end
b=b+j
end
end

Q.4 The number of combinations C n,r of taking r objects out of n objects is given by [ 5 Marks]
!
, = !( − )!

(a) Determine how many combinations are possible in a lottery game for selecting 6
numbers that are drawn out of 49 ?
(b) Using the formula, determine the probability of guessing two out of six drawn numbers
6,2 43,4
49,6

Use Built in function factorial

Q.5 Write the commands to find the following from the data set [ 5 Marks]

mat=[4,5,6,7,8;10,11,12,13]

(a) max number (largest number) in both rows

(b) max number (largest number) in four columns
(c) mean in each row
(d) mean in each column
(e) standard deviation in each row

********

TA ZC164 (EC-2 Regular) First Semester 2018-2019 Page 2 of 2

Birla Institute of Technology & Science, Pilani
Work-Integrated Learning Programmes Division
Second Semester 2017-2018
Mid-Semester Test (EC-2 Make-up)

Course No. : TA ZC164

Course Title : COMPUTER PROGRAMMING
Nature of Exam : Closed Book
Weightage : 30%
No. of Pages= 2
Duration : 2 Hours
No. of Questions = 7
Date of Exam : 18/03/2018(FN)
Note:
1. Please follow all the Instructions to Candidates given on the cover page of the answer book.
2. All parts of a question should be answered consecutively. Each answer should start from a fresh page.
3. Assumptions made if any, should be stated clearly at the beginning of your answer.

Q.1. Short Questions [2 x 5 = 10]

(a) Create a 2 × 4 matrix of random integers, each in the range from 10 to 30

randi([10,30],2,4)

(b) If mat = [4 3 1; 2 5 6] what is the output of mat = [2:4; 3:5]

mat =

2 3 4
3 4 5

linspace(-3,-15,5)

(d) Create a 4 × 2 matrix of all zeros and store it in a variable. Then, replace the second
row in the matrix with a 3 and a 6.

a=zeros(4,2)
a(2,:)=[3,6]

(e) Write a command to find the Eigen value of a 2*2 matrix z.

z=[1,2;3,4]
eig(z)

Q.2. Evaluate the value of c if c=a.\b , a and b are matrices shown below [2]

3 1 1 3

a    b   
2
1 2  2

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 1 of 4

B vector will be divided by a

-0.3333 3.0000 -
2.0000 1.000

Q.3. What is the output of the following output statement [2]

i) fprintf(‘%3.2f\n’, 1234.5678)

Solution:

1234.57 it will take only 2 decimal points

ii) for i=0:10:100

if mod(i,10) ==0
disp(i)
else
continue;
end
end

Solution:
It will print 0 to 100 in the gap of 10 such as 0,10,20,20,40,50,60,70,80,90,100

Use MATLAB to find (i) compression x in each spring and (ii) Potential energy stored in
each spring

Solution:

>> % force f >>f=[11, 7, 8, 10,

9] >>k=[1000, 600,
900,1300,700] >>x=f./k;
>>pe=k.*x.^2/2;

Give Partial marks for steps

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2 of 4

[4]
TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 2

Solution:

a= v.*2;
b=v.^3;
c=v.^v;
d=a.^0.25

Script

% Cloud Classification
x=input('Enter the height of the cloud in feet')
if (x>=0 & x<6500)
disp('Low cloud')
elseif (x>=6500 & x<20000)
disp('Medium cloud')
elseif (x>2000)
disp('high cloud')
end

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 3 of 4

Solution

function average=overall_avg(x)
n=length(x);
sum=0;
for i=1:n
sum=sum+i;
average=sum/n;
end

>> overall_avg([1,2,3,4,5,6,7,8,9,10]

) ans =

5.5000

********

TA ZC164 (EC-2 Make-up) Second Semester 2017-2018 Page 4 of 4

Physics A Level Syllabusdec2024final
100% (2)
Physics A Level Syllabusdec2024final
56 pages
Solutions Manual Matlab A Practical Introduction Programming Problem Solving 2nd Edition Stormy Attaway
73% (41)
Solutions Manual Matlab A Practical Introduction Programming Problem Solving 2nd Edition Stormy Attaway
35 pages
Palm ch4
100% (1)
Palm ch4
36 pages
Assignment & Quiz (Matlab)
100% (1)
Assignment & Quiz (Matlab)
24 pages
IB Physics IA Ideas
38% (8)
IB Physics IA Ideas
1 page
20ec0454 Matlab Programming Question Bank
No ratings yet
20ec0454 Matlab Programming Question Bank
6 pages
Hydraulics
No ratings yet
Hydraulics
18 pages
Mechanical Engineering Fluid-Mechanics Differential-Equations Notes PDF
No ratings yet
Mechanical Engineering Fluid-Mechanics Differential-Equations Notes PDF
167 pages
Study Guide
100% (1)
Study Guide
20 pages
Matlab Marina For Loops Exercises
100% (1)
Matlab Marina For Loops Exercises
4 pages
Sol Matlab
No ratings yet
Sol Matlab
10 pages
DSP Lab 1 Fall 20.PDF NEW
No ratings yet
DSP Lab 1 Fall 20.PDF NEW
12 pages
Solved Problems: Problem 1
No ratings yet
Solved Problems: Problem 1
48 pages
Lab Mariano2023 en
No ratings yet
Lab Mariano2023 en
64 pages
Lecture1 - Modelling-Computers-Error Analysis - 19feb2024 - NK - v1
No ratings yet
Lecture1 - Modelling-Computers-Error Analysis - 19feb2024 - NK - v1
75 pages
Matlab Note: by Abdullah (Student of B S Maths 4 Semester GPG College Timergara Lower Dir) - Date: 9/5/2019
No ratings yet
Matlab Note: by Abdullah (Student of B S Maths 4 Semester GPG College Timergara Lower Dir) - Date: 9/5/2019
66 pages
Intoduction To Atlab
No ratings yet
Intoduction To Atlab
59 pages
DSP All Manuals
No ratings yet
DSP All Manuals
78 pages
Intoduction To Atlab: M. Subramanian
No ratings yet
Intoduction To Atlab: M. Subramanian
58 pages
Introduction To Matlab, Signal Processing & Speech Signal Processing
No ratings yet
Introduction To Matlab, Signal Processing & Speech Signal Processing
49 pages
Numerical Methods Lab: List of Experiments
No ratings yet
Numerical Methods Lab: List of Experiments
27 pages
Assignment Final
No ratings yet
Assignment Final
35 pages
Code and Answers-VBKM
No ratings yet
Code and Answers-VBKM
23 pages
Art's Process Spec Sheets
No ratings yet
Art's Process Spec Sheets
19 pages
Solutions To Homework 6o
No ratings yet
Solutions To Homework 6o
7 pages
Homework Solutions - MATLAB
No ratings yet
Homework Solutions - MATLAB
27 pages
4th Assignment Bangla College
No ratings yet
4th Assignment Bangla College
29 pages
Assignment 1 Answers PDF
No ratings yet
Assignment 1 Answers PDF
16 pages
Matlab: BY K. Kalyan Kumar
No ratings yet
Matlab: BY K. Kalyan Kumar
30 pages
Final Second Sugar Exam 1
No ratings yet
Final Second Sugar Exam 1
4 pages
Book - Kollicoat PDF
No ratings yet
Book - Kollicoat PDF
229 pages
Week - 2 Question Sheet
No ratings yet
Week - 2 Question Sheet
7 pages
Home Work 2
No ratings yet
Home Work 2
25 pages
Matlab Code 3
No ratings yet
Matlab Code 3
29 pages
University of Kirkuk Mechanical Engineering Department Second Stage - Morning Study
No ratings yet
University of Kirkuk Mechanical Engineering Department Second Stage - Morning Study
16 pages
Se - 1
No ratings yet
Se - 1
23 pages
4rth Semester Matlab Assignment Solved 70 Question
No ratings yet
4rth Semester Matlab Assignment Solved 70 Question
19 pages
Simulation Lab Unit I Solution
No ratings yet
Simulation Lab Unit I Solution
12 pages
PracticeFinal Fall2023
No ratings yet
PracticeFinal Fall2023
17 pages
ASTM D88-68 - Saybolt Viscosity
No ratings yet
ASTM D88-68 - Saybolt Viscosity
9 pages
Principles of Matlab (Fall 07) Workout #1
No ratings yet
Principles of Matlab (Fall 07) Workout #1
11 pages
C Se 200 Practice Midterm
No ratings yet
C Se 200 Practice Midterm
12 pages
ENEL2CMH1 - Applied Computer Methods
No ratings yet
ENEL2CMH1 - Applied Computer Methods
9 pages
244 Cheat Sheet
No ratings yet
244 Cheat Sheet
4 pages
1 2 Basics Input Output On Matlab Quiz
No ratings yet
1 2 Basics Input Output On Matlab Quiz
10 pages
Assignment-3 Tanmay
No ratings yet
Assignment-3 Tanmay
7 pages
Matlab Practicals
No ratings yet
Matlab Practicals
17 pages
251 Math Midterm PDF
No ratings yet
251 Math Midterm PDF
10 pages
MEE1006 - MATLAB For Engineers - SET-3 - Solutions
No ratings yet
MEE1006 - MATLAB For Engineers - SET-3 - Solutions
7 pages
DSP Lab 1
No ratings yet
DSP Lab 1
9 pages
nUMERICAL 01
No ratings yet
nUMERICAL 01
9 pages
Ibtasam Mateen (18-Arid-5157) Practical
No ratings yet
Ibtasam Mateen (18-Arid-5157) Practical
6 pages
Help Commands - Name I.E Gives Help On How To Use A Command
No ratings yet
Help Commands - Name I.E Gives Help On How To Use A Command
4 pages
EECS1560 Test 1 - Sample
No ratings yet
EECS1560 Test 1 - Sample
7 pages
Syllabus
No ratings yet
Syllabus
22 pages
GENN004 Spring2012 Midterm1
No ratings yet
GENN004 Spring2012 Midterm1
3 pages
Lab 1
No ratings yet
Lab 1
4 pages
Bhoomi Tiwari 2427030380 Assignment-3 MEE1006
No ratings yet
Bhoomi Tiwari 2427030380 Assignment-3 MEE1006
10 pages
FM & F Machines 2013
No ratings yet
FM & F Machines 2013
355 pages
Worksheet 1 - Recap of Algorithms: Advanced MATLAB
No ratings yet
Worksheet 1 - Recap of Algorithms: Advanced MATLAB
5 pages
Practice Questions
No ratings yet
Practice Questions
3 pages
Assignment-3 MEE1006
No ratings yet
Assignment-3 MEE1006
2 pages
Pe ZC164 Ec-2r First Sem 2022-2023
No ratings yet
Pe ZC164 Ec-2r First Sem 2022-2023
2 pages
Basics of Programing Languages nugpSecQKH
No ratings yet
Basics of Programing Languages nugpSecQKH
5 pages
VIN Solenoid Valve DSL103series
No ratings yet
VIN Solenoid Valve DSL103series
14 pages
The Effect of Flow Geometry On Shear Induced Particle Segregationa and Resuspension
No ratings yet
The Effect of Flow Geometry On Shear Induced Particle Segregationa and Resuspension
294 pages
Migration of Oil
No ratings yet
Migration of Oil
12 pages
Wave Propagation Theories &applications by Yi Zheng
No ratings yet
Wave Propagation Theories &applications by Yi Zheng
392 pages
Effects of Thermal Treatments On Five Different Waste Activated Sludge Samples Solubilisation, Physical Properties and Anaerobic Digestion
No ratings yet
Effects of Thermal Treatments On Five Different Waste Activated Sludge Samples Solubilisation, Physical Properties and Anaerobic Digestion
9 pages
Lubricantes para Cojinetes de Molinos - Bernardo Vianna
No ratings yet
Lubricantes para Cojinetes de Molinos - Bernardo Vianna
22 pages
InstruCalc Programs4
No ratings yet
InstruCalc Programs4
2 pages
Sanders On
No ratings yet
Sanders On
18 pages
504-C-002 Deethanizer HDS Rev 0
No ratings yet
504-C-002 Deethanizer HDS Rev 0
9 pages
10 Viscoelasticity 03 Rheological
No ratings yet
10 Viscoelasticity 03 Rheological
13 pages
Iso 5636 3 2013 en PDF
No ratings yet
Iso 5636 3 2013 en PDF
8 pages
Elasticity - Viscosity - Surface Tension
No ratings yet
Elasticity - Viscosity - Surface Tension
3 pages
Transport Phenomena
No ratings yet
Transport Phenomena
21 pages
MRSPTU B. Tech. Aeronautical Engg. (3rd - 4th Sem.) Syllabus 2018 Batch
No ratings yet
MRSPTU B. Tech. Aeronautical Engg. (3rd - 4th Sem.) Syllabus 2018 Batch
29 pages
Ebsray Pumps: V Series Model V40
No ratings yet
Ebsray Pumps: V Series Model V40
8 pages
Synthesis and Characterization of Two-Component Acrylic Base
No ratings yet
Synthesis and Characterization of Two-Component Acrylic Base
4 pages
Knight 2012
No ratings yet
Knight 2012
19 pages
Mecanismo de Consumo de Energía en Un Votator Trommelen 1971 PDF
No ratings yet
Mecanismo de Consumo de Energía en Un Votator Trommelen 1971 PDF
10 pages
Waukesha Cherry-Burrell: How To Read and Use Centrifugal Pump Curve Performance Data
No ratings yet
Waukesha Cherry-Burrell: How To Read and Use Centrifugal Pump Curve Performance Data
4 pages
PDS - EQUINAS 8 - 2023.v3
No ratings yet
PDS - EQUINAS 8 - 2023.v3
1 page
Assignment 1: Assigned Date: 31 January 2020 Due Date: 12 February 2020 Max Marks: 100
No ratings yet
Assignment 1: Assigned Date: 31 January 2020 Due Date: 12 February 2020 Max Marks: 100
2 pages