1.

Write a MATLAB Program for find SUM, SUB, MUL, DIVISION

of given two Matrices

a = input('Enter the First Matrix \n');

b = input('Enter the Second Matrix \n');
c = a+b;
d = a-b;
e = a*b;
f = a/b;
disp('A Matrix');
disp(a);
disp('B Matrx');
disp(b);
disp('Matrix Addition: ');
disp(c);
disp('Matrix Subtraction:');
disp(d);
disp('Matrix Multilication: ')
disp(e);
disp('Matrix Division: ');
disp(f);

OUTPUT

Enter the First Matrix

[3 4 2; 2 1 5; 5 3 2]

Enter the Second Matrix

[7 6 5; 2 3 4; 5 2 1]

A Matrix

3 4 2

2 1 5

5 3 2

B Matrx

7 6 5
 2 3 4

5 2 1

Matrix Addition:

10 10 7

4 4 9

10 5 3

Matrix Subtraction:

-4 -2 -3

0 -2 1

0 1 1

Matrix Multilication:

39 34 33

41 25 19

51 43 39

Matrix Division:

1.9444 -1.5556 -1.5000

-2.6111 3.8889 2.5000

1.3889-0.1111 0.5000
2. Write a MATLAB Program for find element by Element matrix operations of
MUL,DIVISION,EXPONENT and Transpose of of given two Matrices

a = input('Enter the First Matrix \n');

b = input('Enter the Second Matrix \n');

c = a .*+ b;

d = a ./ b;

e = a .^ b;

f = a';

f2 = b';

disp('A Matrix');

disp(a);

disp('B Matrx');

disp(b);

disp('Matrix Multiplication Element by Element: ');

disp(c);

disp('Matrix Division Element by Element:');

disp(d);

disp('Matrix Exponentiation Element by Element: ')

disp(e);

disp('Transpose of A Matrix: ');

disp(f);

disp(f2);
OUTPUT

Enter the First Matrix

[3 4 5; 6 5 7; 1 2 3]

Enter the Second Matrix

[6 7 8; 4 3 6; 9 8 7]

A Matrix

3 4 5

6 5 7

1 2 3

B Matrx

6 7 8

4 3 6

9 8 7

Matrix Multiplication Element by Element:

18 28 40

24 15 42

9 16 21

Matrix Division Element by Element:

0.5000 0.5714 0.6250

1.5000 1.6667 1.1667

 0.1111 0.2500 0.4286

Matrix Exponentiation Element by Element:

729 16384 390625

1296 125 117649

1 256 2187

Transpose of A Matrix:

3 6 1

4 5 2

5 7 3

Transpose of B Matrix:

6 4 9

7 3 8

8 6 7
3. Write a MATLAB Program for finding the Area and Circumference of a Circle?

r = input('Enter the Radius');

A = pi * r^2;
C= 2*pi*r;
fprintf('Radius of Circle = %d\n',r);
fprintf('Area of Circle = %f\n',A);
fprintf('Circumference of Circle = %d\n',C);

OUTPUT
Enter the Radius5

Radius of Circle = 5

Area of Circle = 78.539816

Circumference of Circle = 3.141593e+001

4. Write a MATLAB Program for finding Sum and Sum of Squares up to given numbers?

i = 1;
sum = 0;
sumsq = 0;
n = input('Enter How many Numbers');
for i = 1 : n
sum = sum + i;
sumsq= sumsq + i*i;
fprintf('Number = %d\t Square Number = %d\n',i,i*i);
end
fprintf('Sum of %d numbers = %d\n',n,sum);
fprintf('Sum of squares of %d numbers = %d\n',n,sumsq);

OUTPUT
1

Enter How many Numbers5

Number = 1 Square Number = 1

Number = 2 Square Number = 4

Number = 3 Square Number = 9

Number = 4 Square Number = 16

Number = 5 Square Number = 25

Sum of 5 numbers = 15

Sum of squares of 5 numbers = 55

Enter How many Numbers10

Number = 1 Square Number = 1

Number = 2 Square Number = 4

Number = 3 Square Number = 9

Number = 4 Square Number = 16

Number = 5 Square Number = 25

Number = 6 Square Number = 36

Number = 7 Square Number = 49

Number = 8 Square Number = 64

Number = 9 Square Number = 81

Number = 10 Square Number = 100

Sum of 10 numbers = 55

Sum of squares of 10 numbers = 385

5. Write a MATLAB Program for finding Area of Circle ,Triangle and Rectangle using
MENU options?

ch = menu('Area','Circle','Triangle','Rectangle');

if ch == 1

r = input('Enter Radius');

A1 = pi*r^2;

fprintf('Radius of Circle = %d \n',r);

fprintf('Area of Circle = %f \n',A1);

elseif ch == 2

a = input('Enter a Value');

b = input('Enter b Value');

c = input('Enter c Value');

s = (a+b+c)/2;

A= sqrt(s*(s-a)*(s-b)*(s-c));

fprintf('Side of a Triangle = %f \n',s);

fprintf('Area of Triangle = %d \f',A);

else

l = input('Enter Length');

br = input('Enter breadth');

A2=l*br;

fprintf('Area of Rectangle = %d \n',A2);

end
OUTPUT

Enter Radius 6

Radius of Circle = 6

Area of Circle = 113.097336

Enter a Value 3

Enter b Value 4

Enter c Value 6

Side of a Triangle = 6.500000

Area of Triangle = 5.332682e+000

Enter Length 10

Enter breadth 25

Area of Rectangle = 250

6. Write a MATLAB Program for finding Arithmetic operations of Two given Numbers
using SWITCH?

a = input('Enter a Value');

b = input('Enter b Value');

disp('1.Addtion');

disp('2.Subtraction');

disp('3.Multiplication');

disp('4.Division');

ch = input('Enter Your Choice');

switch(ch)

case 1

c = a+b;

fprintf('Addition of Two Numbers = %d\n',c);

case 2

c = a-b;

fprintf('Subtraction of Two Numbers = %d\n',c);

case 3

c = a*b;

fprintf('Multiplication of Two Numbers = %d\n',c);

case 4

c = a/b;

fprintf('Division of Two Numbers = %d\n',c);

otherwise

disp('Invalid Choice');

end
OUTPUT

Enter a Value 8

Enter b Value 4

1.Addtion

2.Subtraction

3.Multiplication

4.Division

Enter Your Choice 1

Addition of Two Numbers = 12

>>

Enter a Value 56

Enter b Value 43

1.Addtion

2.Subtraction

3.Multiplication

4.Division

Enter Your Choice 4

Division of Two Numbers = 1.302326e+000

>>

Enter a Value 45

Enter b Value 65

1.Addtion

2.Subtraction
3.Multiplication

4.Division

Enter Your Choice 3

Multiplication of Two Numbers = 2925

>>

Enter a Value 67

Enter b Value 32

1.Addtion

2.Subtraction

3.Multiplication

4.Division

Enter Your Choice2

Subtraction of Two Numbers = 35

7. (a)Compute the Integral

𝟑/𝟐 𝟐
I.∫𝟏/𝟐 𝒆−𝒙 𝒅𝒙
𝟏 𝟐
II. ∫−𝟏 ∫𝟎 𝟏 − 𝟔𝒙𝟐 𝒚 𝒅𝒙𝒅𝒚
function y = int1(x);
y=exp(-x.^2);
end

>> y=quad('int1',1/2,3/2)

y=

0.3949

>> F = inline('1-6*x.^2*y');

>> I=dblquad(F,0,2,-1,1)

I=

4.0000
(b).Solve the First order Linear Differential Equation

I. dx/dt = x+t
With the initial condition x0=0
II.Solve the dy/dx = (x3-2y)/x
With the initial condition 1<x<3 y0=4.2

function xdot = simpode(t,x);

xdot = x+t;
end

>> tspan=[0 2];

>> x0=0;

>> [t,x]=ode23('simpode',tspan,x0)

>> plot(t,x)

>> xlabel('t')

>> ylabel('x')

Ode23 graph
II.

>> f=@(x,y) (x^3-2*y)/x

>> [x,y]=ode45(f,[1:0.01:3],4.2)

>>plot(x,y)

Ode45 graph
 8. Write a Matlab program for Linear, Quadratic and Polynomial curve
fits for the following data and also estimate sales in the year
2018.Compare these 3 fits.Which is better?
Year 2000 2002 2004 2006 2008 2010 2012
Sales(in 45 32 62 75 95 114 144
lackhs)
x=[2000:2:2012];
y=[45,32,62,75,95,114,144];
p=polyfit(x,y,1);
y_fit=polyval(p,x);
plot(x,y,'ro',x,y_fit);
title('Sales of KFC company');
Xlabel('years');
Ylabel('Sales (in Thousands)');

x=[2000:2:2012];
y=[45,32,62,75,95,114,144];
p=polyfit(x,y,2);
y_fit=polyval(p,x);
plot(x,y,'ro',x,y_fit);
title('Sales of KFC company');
Xlabel('years');
Ylabel('Sales (in Thousands)');

x=[2000:2:2012];
y=[45,32,62,75,95,114,144];
p=polyfit(x,y,3);
y_fit=polyval(p,x);
plot(x,y,'ro',x,y_fit);
title('Sales of KFC company');
Xlabel('years');
Ylabel('Sales (in Thousands)');

Output
 The estimated sales in the year 2018
 By using linear curve is 187
 By using Quadratic curve is 253
 By using 3rd degree polynomial is 176
Therefore Quadratic curve is best fit for the given data

9. Create Bar , Pie and Area Chart of the following Data

continent Asia Europe Africa N.America S.America

Population 3332 696 694 437 307

cont=char('Asia', 'Europe','Africa','N.America','S.America');
pop=[3332;696;694;437;307];
barh(pop)
for i=1:5,
gtext(cont(i,:))
end
xlabel('Population in millions')
title('world population(1992)','fontsize',18)

cont=char('Asia', 'Europe','Africa','N.America','S.America');
pop=[3332;696;694;437;307];
pie(pop)
for i=1:5,
gtext(cont(i,:))
end
title('world population(1992)','fontsize',18)

x=linspace(-3*pi,3*pi,100);
y=-sin(x)./x;
area(x,y)
xlabel('x'),ylabel('sin(x)./x')
hold on
x1=x(46:55);
y1=y(46:55);
area(x1,y1,'facecolor','y')
10. Write a MATLAB Program to generate an overlay plot with the Line Command for

y1=sin(t)

y2=t

(t-t3)
y3 = ---------

( 6+t5/120)

t=linspace(0,2*pi,100);
y1=sin(t);y2=t;
y3=t-(t.^3)/6+(t.^5)/120;
plot(t,y1)
line(t,y2,'linestyle','--')
line(t,y3,'marker','o', ...
'linestyle','none')
axis([0 5 -1 5])
xlabel('t')
ylabel(' Approximation of sin(t)')
title('Fun with sin(t)')
legend('sin(t)','Linear Approx','Fifth-order-Approx')

OUTPUT
 11. Write a MATLAB Program using Spline and Pchip Function

Time (sec) 0 10 20 30 40 50 60 70 80 90
Speed 45 32 0 0 7 12 20 15 29 55
Can taken 1 ½ min from this end of the road to that end of the road

%Speed in 1.5 minutes

>>t=0:10:90;

>> s=[45 32 0 0 7 12 20 15 29 55];

>>plot(t,s,’bo’); hold on;

%interpolation using spline

>>tI=0:90;

>>sI=spline(t,s,tI);

>>plot(tI,sI,’-k’);

%interpolation using PCHIP

sI_pc=pchip(t,s,tI);

plot(tI,sI_pc,’--r’);

OUTPUT
12. Write a MATLAB Program to find the Eigen values and Eigen vectors by using Gaussian
Elimination method?

5x =3y-2z+10

8y+4z=3x+20

2x+4y-9z=9

>> A=[5 -3 2; -3 8 4; 2 4 -9];

>> B=[10 ;20;9];

>>x=A\B

>> C=A*x %check the solution%

Gaussian Elimination method

>>C1=[A B] %form the augmented matrix%

>>Cr=rref(C1) %row reduce the augmented matrix%

Eigen values and Eigen Vectors

>> [V ,D]=eig(A)

Output
 3.4442
x =[3.1982]
1.1868

10
C =[20]
9

C1=

5 -3 2 10

-3 8 4 20

2 4 -9 9

Cr =

1 0 0 3.4442

0 1 0 3.1982

0 0 1 1.1868

V=

-0.1709 0.8729 0.4570

-0.2365 0.4139 -0.8791

0.9565 0.2583 -0.1357

D=

-10.3463 0 0

0 4.1693 0

1 0 10.1770