Matlab
Matlab
Department of Mechanical
Guided By
Prof. Dr. Veera Venkata Nagaraju
Rushikesh R Mankar(23MEM6R05)
Slide 12
Problem 6
syms x y z a
f1=2*x+3*y-4*z;
f2=x-2*y+5*z;
a=f2-f1
Solution-
a = 9*z - 5*y – x
%problem 7
syms a b c
f1=2*a+3*b;
f2=a+b;
c=f1*f2
Solution-
c = (a + b)*(2*a + 3*b)
%problem 8
syms x y
f1=3*x-2*y^2+4*x*y;
f2=2*x-5*y;
c=f1*f2
Solution-
c = (2*x - 5*y)*(- 2*y^2 + 4*x*y + 3*x)
%problem 1
syms x y z b a
b=(x^2*y^3*z)*(x^3*y*z^2)
a=simplify(b)
subs(a,[x y z],[0.5 2 3])
a = x^5*y^4*z^3
ans = 27/2
%problem 29
syms x y b a
factor(a*x-a*y+b*x-b*y)
ans =[x - y, a + b]
%problem 30
syms x y b a
factor(2*a*x-3*a*y+2*b*x-3*b*y)
ans = [2*x - 3*y, a + b]
%problem 31
syms x
factor(x^3+3*x^2-x-3)
ans = [x + 3, x - 1, x + 1]
%problem 3
syms x y
simplify((x^3+y^3)/(x+y))
ans = x^2 - x*y + y^2
Slide 13
%problem 1
syms x y a b
y=7*x^4;
a=diff(y,x,2);
b=diff(y,x,3);
a = 84*x^2
b = 168*x
%problem 2
syms x y a b
y=sqrt(x);
a=diff(y,x,2)
b=diff(y,x,3)
a = -1/(4*x^(3/2))
b = 3/(8*x^(5/2))
%problem 3
syms t y a b
y=sqrt(t^3);
a=diff(y,t,2)
b=diff(y,t,3)
a = (3*t)/(t^3)^(1/2) - (9*t^4)/(4*(t^3)^(3/2))
b = 3/(t^3)^(1/2) - (27*t^3)/(2*(t^3)^(3/2)) + (81*t^6)/(8*(t^3)^(5/2))
%problem 4
syms x y a b
y=6+(1/x^3)
a=diff(y,x,2)
b=diff(y,x,3)
%problem 5
syms x y a b
y=3*x-(1/sqrt(x))+(1/x)
a=diff(y,x,2)
b=diff(y,x,3)
%problem 6
syms x y a b
y=(5/x^2)-(1/sqrt(x^7))+2
a=diff(y,x,2)
b=diff(y,x,3)
%problem 7
syms t y a b
y=3*(t-2)^2
a=diff(y,t,2)
a=6
%problem 9
syms x y a b
y=(2/5)*x^3-(4/x^3)+4*sqrt(x^5)+7
a=diff(y,x)
a=
%problem 10
syms t y a b
y=5*t+(1/sqrt(t^3))
a=diff(y,t)
%problem 11
syms x y
y=((x+2)^2)/x
diff(y,x)
%problem 13
syms x y
y=7*sin(2*x)-3*cos(4*x)
diff(y,x)
Slide 14
%problem 1
syms x y a b
y=3*x^4+2*x^3-3*x+2
a=diff(y,x,2)
b=diff(y,x,3)
%problem 2
syms t y a
y=((2*t^2)/5)-(1/t^3)+(3/t)-sqrt(t)+1
a=diff(y,t,2)
subs(a,t,1)
%problem 3
syms t y1 y2 a b
y1=3*sin(2*t)+cos(t)
y2=2*log(4*t)
a=diff(y1,t,2)
b=diff(y2,t,2)
%problem 4
syms x y1 y2 a b
y1=2*cos(x)^2
y2=(2*x-3)^4
a=diff(y1,x,2)
b=diff(y2,x,2)
Slide 15
%problem 1
clear all
syms x y a b
y=(x^2-4*x)*(x^2-4*x+1)-20
a=expand(y)
b=y/(x^2-4*x+4)
c=factor(b)
d=solve(x^2-4*x+5)
e=subs(b,2-j)
Slide 16
%problem 3
syms x y z
exp1=x+3*y-z-2;
exp2=x-y+z-3;
exp3=3*x-5*y-4;
[x,y,z]=solve(exp1,exp2,exp3)
Solution-
x = 33/16
y = 7/16
z = 11/8
%problem 5(i)
syms x
y=exp(-x)
a=int(y,0,Inf)
%problem 5(ii)
syms x
y=exp(-x^2)
b=int(y,0,Inf)
%problem 5(iii)
syms x
y=exp(-x^2)*sin(x)
b=int(y,0,Inf)
a=1
b = -(pi^(1/2)*erf(1i/2)*exp(-1/4)*1i)/2
Slide 17
syms f(x)
syms x
f(x)=x^3+2*x^2-x+3
dfdt=diff(f,x,1)
Solution-
dfdt(x) = 3*x^2 + 4*x – 1
Slide 18
syms f(x)
syms x
f(x)=0.2+25*x-200*x^2+675*x^3-900*x^4+400*x^5
dfdt=diff(f,x,1)
Sol-
dfdt(x) = 2000*x^4 - 3600*x^3 + 2025*x^2 - 400*x + 25
%problem 1
syms x y t
y1=sin(3*x)+sin(x)
a=diff(y1,x)
%problem 2
y2=0.5*(sin(9*t)-sin(7*t))
b=diff(y2,t)
%problem 3
y3=cos(5*t)+cos(3*t)
c=diff(y3,t)
%problem 4
y4=(cos(5*t)-cos(t))/8
d=diff(y4,t)
Slide 22
%problem 1
clc
clear all
y1=[1 4 -32]
a=roots(y1)
%problem 2
y2=[1 0 -16]
b=roots(y2)
a =-8
4
b= 4
-4
%problem 3
syms x
solve((x+2)^2==16)
ans =-6
2
%problem 4
clc
clear all
y4=[2 -1 -3]
a=roots(y4)
%problem 5
y5=[6 -5 1]
b=roots(y5)
%problem 6
y6=[10 3 -4]
c=roots(y6)
%problem 7
y7=[1 -4 4]
d=roots(y7)
%problem 8
y8=[21 -25 -4]
e=roots(y8)
%problem 9
y9=[6 -5 -4]
f=roots(y9)
%problem 10
y10=[8 2 -15]
h=roots(y10)
Sol-
a = 1.5000
-1.0000
b = 0.5000
0.3333
c = -0.8000
0.5000
d =2
2
e = 1.3333
-0.1429
f = 1.3333
-0.5000
h = -1.5000
1.2500
%problem 2
syms x y
a=factor(x^2-6*x+9)
b=factor(4*x^2-25)
Sol-
a = [x - 3, x - 3]
b = [2*x - 5, 2*x + 5]
% problem 3
syms x
a=factor(4*x^2+8*x+3)
b=factor(15*x^2+2*x-8)
Sol-
a = [2*x + 3, 2*x + 1]
b = [5*x + 4, 3*x - 2]
SLIDE 23
clear all
syms x
a=solve(x^3-1==0)
b=[1 -1 -5 2]
c=roots(b)
syms x
a=[1 -2 -2 2];
b=round(roots(a),1)
c=[2 -1 -9.08 8.28];
d=round(roots(c),1)
e=[8 36 54 27];
f=round(roots(e),1)
Sol-
b = 2.5000
-1.2000
0.7000
d = -2.3000
1.7000
1.0000
f =-1.5000
-1.5000
-1.5000
SLIDE 26
clear all
syms x
p1=[1 8 5 4 3 2 1 1];
a=roots(p1)
p2=[1 -7 7 -10 -8 7 15];
b=roots(p2)
p3=[1 -8 10 12 8 -15];
c=roots(p3)
Sol-
clear all
syms x
p4=[1 7 12 -25 8];
a=roots(p4)
p5=[1 15 -23 105];
b=roots(p5)
p6=[1 -18 23];
c=roots(p6)
p7=[1 0 0 7];
d=roots(p7)
SLIDE 29
clear all
%problem 8(a)
syms x
a=solve(8*x+3==0)
% problem 8(b)
format long
b=round(solve(8*x+3==0),15)
%problem 8(c)
syms x p q
c=solve(x^3+p*x+q==0)
%problem 8(d)
d=solve(exp(x)-8*x+4)
SOL-
a = -3/8
b = -0.375
c = root(z^3 + p*z + q, z, 1)
root(z^3 + p*z + q, z, 2)
root(z^3 + p*z + q, z, 3)
%problem 9(b)
syms y x
y=sin(1/x^2);
ezplot(y,[-2,2])
grid on
xlabel('x')
ylabel('y')
title('i) using ezplot')
x=linspace(-2,2)
y=sin(1./x.^2);
plot(x,y)
grid on;
xlabel('x')
ylabel('y')
title('ii) Using plot')
%problem 9(c)
syms y x
y=tan(x/2);
ezplot(y,[-pi,pi,-10,10])
grid on
xlabel('x')
ylabel('y')
%problem 9(d)
x=linspace(-4,4);
y1=exp(-x.^2);
y2=x.^4-x.^2;
plot(y1,x,y2,x)
xlabel('x')
ylabel('y')
grid on
SLIDE 37
syms x t
f=4+(3/7)*x-6*x^2
I=int(f,x)
%problem 3(a)
f1=(2*x^3-3*x)/(4*x)
I1=int(f1,x)
%problem 3(b)
f2=(1-t)^2
I2=int(f2,t)
%problem 4
syms x t
f=(3/x^2)
I=int(f,x)
%problem 6
f1=(-5/9)*(t^(-3/4))
I1=int(f1,t)
%problem 8(a)
f2=4*cos(3*x)
I2=int(f2,x)
%problem 8(b)
f3=5*sin(2*t)
I2=int(f3,t)
%problem 10(a)
syms x t
f=5*exp(3*x)
I=int(f,x)
%problem 10(b)
f1=(2/3)*exp(-4*t)
I1=int(f1,t)
%problem 3(a)
syms x t
f=((3*x^2-5*x)/x)
I=int(f,x)
%problem 3(b)
f1=(2+t)^2
I1=int(f1,t)
Sol-
%problem 4(a)
syms x t
f=(4/3)*x^(-2)
I=int(f,x)
f = 4/(3*x^2)
I -4/(3*x)
SLIDE 38
clear all
%problem 19.2
y=inline('1-exp(-x)')
z=quad(y,0,4)
%problem 19.3
y1=@(x) 8+4*cos(x)
z1=quad(y1,0,pi/2)
%problem 19.4
y2=@(x) 1-x-4*x.^3+2*x.^5
z2=quad(y2,-2,4)
%problem 19.6
clear all
a=@(x,y) x.^2-3*y.^2+x*y.^3
I=dblquad(a,0,4,-2,2)
SOL-
a = function_handle with value:
@(x,y)x.^2-3*y.^2+x*y.^3
I = 21.3333
%problem 19.7
clear all
syms x y z
I=int(int(int(x^3-2*y*z,x,-1,3),0,6),-4,4)
Sol-
I = 960
%problem
clear all
y=@(x) x.*exp(2*x)
I=quad(y,0,3)
Sol-
@(x)x.*exp(2*x)
I = 504.5360
SLIDE 39
%problem 2
clear all
syms x y
I1=int((3-x^2),x,-1,2)
I2=int((x^2-4*x+3),x,1,3)
Sol-
I1,I2
I1 = 6
I2 = -4/3
%problem 3
clear all
syms x y1
y1=(3/2)*cos(x)
I1=int(y1,0,pi)
y2=4*cos(x)
I2=int(y2,0,pi/2)
Sol-
I1 = 0
I2 = 4
%problem 4
clear all
syms x y1
y1=2*sin(2*x)
I1=int(y1,pi/6,pi/3)
y2=@(t) 3.*sin(t)
I2=quad(y2,0,2)
Sol-
y1 = 2*sin(2*x)
I1 = 1
@(t)3.*sin(t)
I2 = 4.2484
%problem 5
clear all
syms x y1
y1=@(x) 5.*cos(3*x)
I1=quad(y1,0,1)
y2=@(t) 3*sec(2*t).^2
I2=quad(y2,0,pi/6)
Sol-
y1 = function_handle with value:
@(x)5.*cos(3*x)
I1 = 0.2352
@(t)3*sec(2*t).^2
I2 = 2.5981
%problem 5
clear all
syms x y1 t
y1=1/(sin(4*t)^2)
I1=int(y1,1,2)
y2=3*sin(2*x)-2*cos(3*x)
I2=int(y2,pi/4,pi/2)
Sol-
y1 = 1/sin(4*t)^2
I1 = Inf
y2 = 3*sin(2*x) - 2*cos(3*x)
I2 = 2^(1/2)/3 + 13/6
%problem 7
y1=@(t) 3*exp(3*t)
I1=quad(y1,0,1)
y2=@(x) (2/3)*exp(-2*x)
I2=quad(y2,-1,2)
y1 = function_handle with value:
@(t)3*exp(3*t)
I1 = 19.0855
@(x)(2/3)*exp(-2*x)
I2 = 2.4569
%problem 1
syms x y
I1=int(int(exp(x+y),0,log(y)),1,log(8));
round(I1,4)
ans = 3.3538
%problem 2
syms x y
I1=int(1/(1+x^2+y^2),0,sqrt(1+x^2))
int(I1,0,1)
clear all
%problem 3
syms x y a
I1=int(int(cos(x+y),0,pi/2),0,pi/2)
%problem 4
I2=int(int(sqrt(a^2-x^2-y^2),0,sqrt(a^2-y^2)),0,a)
Sol-
I1 = 0
I2 = (pi*a^3)/6
%problem 8
syms x y z
I1=int(int(int((x*y*z),0,y*z),1,z),0,2)
Sol-
I1 = 7/2
%problem 9
syms x y z
I1=int(int(int(cos(z/x),0,x*y),x,pi/2),0,pi/2)
SLIDE 48
clear all
%problem 1
ode1 = @(t, x, y, z) -10*x + 10*y;
ode2 = @(t, x, y, z) -y - x*z;
ode3 = @(t, x, y, z) (-8/3)*z + x*y - (224/3);
initial_conditions = [-8, 8, 27];
[t, solution] = ode45(@(t, u) [ode1(t, u(1), u(2), u(3)); ode2(t, u(1), u(2), u(3)); ode3(t, u(1),
u(2), u(3))], [0, 10], initial_conditions);
SLIDE 49
%problem 2
clear all
syms y(x)
cond=y(3)==5
a=diff(y,x)==4*x-2
ysol(x)=dsolve(a,cond)
%problem 3
clear all
syms y(x)
a=diff(y,x)==2*x+3
y(x)=dsolve(a)
Sol-
a(x) = diff(y(x), x) == 2*x + 3
ans = C1 + x*(x + 3)
%problem 4
clear all
syms y(x) a
c=diff(y,x,2)+3*diff(y,x)+a*y==0
y(x)=dsolve(c)
y(x)=simplify(y(x))
%problem 5
clear all
syms y(x) a
a=diff(y,x)==(2*x+3*cos(y))/(3*x*sin(y)-2*y)
y(x)=dsolve(a,'Implicit',true)
Sol-
ans = 3*x*cos(y(x)) + y(x)^2 == - x^2 + C1
SLIDE 50
%problem 6
clear all
syms y(x)
a=diff(y,x,2)+8*diff(y,x)+2*y==cos(x)
Dy=diff(y,x)
cond=[y(0)==0,Dy(0)==1]
y(x)=dsolve(a,cond)
simplify(y(x))
%problem 7
clear all
syms y(x)
a=diff(y,x,2)+y==sin(x)
Dy=diff(y,x)
cond=[y(0)==1,Dy(0)==-2]
y(x)=dsolve(a,cond)
simplify(y(x))
%problem 8 i)
syms y(x)
a=diff(y,x)==7*x
y(x)=dsolve(a)
ans = (7*x^2)/2 + C1
syms y(x)
%problem 8 ii)
b=diff(y,x,2)+5*diff(y,x)+2*y==0
ysol(x)=dsolve(b)
%problem 9
syms y(x)
a=x*diff(y,x)==2-4*x^3
y(x)=dsolve(a)
Sol-
y(x) = C1 + 2*log(x) - (4*x^3)/3
SLIDE 51
%problem 10
syms y(x)
a=diff(y,x)==(x-1)^2
y(x)=dsolve(a)
Sol-
%problem 11
syms y(x)
a=a=diff(y,x)==-(1/exp(3*x))
y(x)=dsolve(a)
%problem 12
syms y1(x) y2(x) y(x)
eqn1=diff(y,x)+2*x*y^2==0
y1(x)=dsolve(eqn1)
%problem 13
eqn2=diff(y,x)==exp(3*x+2*y)
y2(x)=dsolve(eqn2)
SLIDE 52
%problem 1
syms y(x)
a=diff(y,x)==cos(4*x)-2*x
y(x)=dsolve(a)
%problem 2
clear all
syms y(x)
b=2*x*diff(y,x)==3-x^3
y1(x)=dsolve(b)
sol-
%problem 3
syms y(x)
a=diff(y,x)+x==3
y(x)=dsolve(a,'y(1)=2')
%problem 4
clear all
syms y(x)
b=3*diff(y,x)+sin(x)==0
y1(x)=dsolve(b,'y(pi/3)=2/3')
sol-
%problem 5
syms y(x)
a=exp(-x)+2==x-3*diff(y,x)
y(x)=dsolve(a,'y(0)=1')
Sol-
%problem 4
clear all
syms y(x)
b=2*y*(1-x)+x*(1+y)*diff(y,x)==0
y1(x)=dsolve(b,'y(1)=1','Implicit',true)
Sol-
y1(x) =log(y(x)) + y(x) == 2*x - 2*log(x) – 1
%problem 5
clear all
syms y(x)
b=(y/x)*diff(y,x)==(y^2+1)/(x^2+1)
y1(x)=dsolve(b,'Implicit',true)
SOL-
log(y(x)^2 + 1) == C1 + log(x^2 + 1)
(y(x)^2 + 1)/y(x) == 0
%problem 6
clear all
syms y(x)
b=x*y==(1-x^2)*diff(y,x)
y1(x)=dsolve(b,'y(0)=1')
Sol-
y1(x) = 1i/(x^2 - 1)^(1/2)
SLIDE 59
% problem 1
clear all
c=0.517;
tau=3;
[t,y]=ode45(@(t,y) (c-y)/tau,[0 25],0.5)
plot(t,y)
grid on
xlabel('t')
ylabel('x')
% problem 1(b)
clear all
[t,y]=ode45(@(t,y) (y+3*t)/t^2,[1 4],-2)
plot(t,y)
grid on
xlabel('t')
ylabel('x')
% problem 1(c)
clear all
[t,y]=ode45(@(t,y) -y*sin(t),[-1 2],1)
plot(t,y)
grid on
axis('equal')
xlabel('t')
ylabel('y')
% problem 2
clear all
[t,y]=ode45(@(t,y) 2*y+t^(-2),[0.1 6],2.4)
plot(t,y)
grid on
xlabel('t')
ylabel('y')
SLIDE 60
% problem 3
clear all
[t,y]=ode45(@(t,y) (-t*y)/(sqrt(2-y^2)),[0 6],1)
plot(t,y)
grid on
xlabel('t')
ylabel('y')
Problem 4
function solve_ode()
tspan = [0, pi/2];
y0 = [1; -1; 0];
[t, y] = ode45(@odefun, tspan, y0);
plot(t, y)
grid on
xlabel('t')
ylabel('y')
legend('y1', 'y2', 'y3')
end
SLIDE 61
syms y x
y=3*sin(x)+7*cos(5*x);
lhs=diff(y,x,2)+y
rhs=-5*cos(2*x);
if(isequal(lhs=rhs))
disp('Yes, the equation is true')
else
disp('No,the equation is false')
end
SOL-