Scribd
Upload
26 views11 pages

MATLAB Assessment: 3: Dave Meet Hareshbhai 19BEC0041 L53 + L54

The document appears to be MATLAB code for solving differential equations both analytically and numerically. The code takes in coefficients and initial conditions for a differential equation, then uses symbolic solving or Laplace transforms to find the general or complete solution which can be plotted.

Uploaded by

Meet Dave
Copyright
© © All Rights Reserved
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
26 views11 pages

MATLAB Assessment: 3: Dave Meet Hareshbhai 19BEC0041 L53 + L54

The document appears to be MATLAB code for solving differential equations both analytically and numerically. The code takes in coefficients and initial conditions for a differential equation, then uses symbolic solving or Laplace transforms to find the general or complete solution which can be plotted.

Uploaded by

Meet Dave
Copyright
© © All Rights Reserved
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/ 11

MATLAB Assessment : 3

Dave Meet Hareshbhai


19BEC0041
L53 + L54

CODE:
clear all
close all
clc
syms c1 c2 x m
F=input('Enter the coefficients [a,b,c]: ');
f=input('Enter the RHS function f(x): ');
a=F(1);b=F(2);c=F(3);
AE=a*m^2+b*m+c;
m=solve(AE);
m1=m(1); m2=m(2);
D=b^2-4*a*c;
if(D>0)
y1=exp(m1*x);y2=exp(m2*x);
elseif (D==0)
y1=exp(m1*x);y2=x*exp(m1*x);
else
alfa=real(m1);beta=imag(m1);
y1=exp(alfa*x)*cos(beta*x);
y2=exp(alfa*x)*sin(beta*x);
end
yc=c1*y1+c2*y2;
fx=f/a;
W=y1*diff(y2,x)-y2*diff(y1,x);
u=int(-y2*fx/W,x);
v=int(y1*fx/W,x);
yp=y1*u+y2*v;
y_gen=yc+yp;
check=input('If the problem has initial conditions then
enter 1 else enter 2: ');
if(check==1)
cn=input('Enter the initial conditions [x0, y(x0),
Dy(x0)]:');
dy_gen=diff(y_gen);
eq1=(subs(y_gen,x,cn(1))-cn(2));
eq2=(subs(dy_gen,x,cn(1))-cn(3));
[c1 c2]=solve(eq1,eq2);
y=simplify(subs(y_gen));
disp('The complete solution is');
disp(y);
ezplot(y, [cn(1),cn(1)+2]);
else
y=simplify(y_gen);
disp('The General Solution is ');
disp(y);
end

Input & Output:


1.

2.
3.

4.

5.
6.
7.
CODE:
clc;
clear all;
syms t s Y
F=input('Enter the coefficients of the differential
equation [a b c]=');
a=F(1);b=F(2);c=F(3);
f=input('Enter the function on RHS=');
fs=laplace(f,t,s);
IC= input('Enter the initial conditions [y0,dy0]');
y0=IC(1);dy0=IC(2);
dy1=s*Y-y0;
dy2=s^2*Y-s*y0-dy0;
eqn=a*dy2+b*dy1+c*Y-fs;
y=solve(eqn,Y);
yt=simplify(ilaplace(y,s,t));
disp(yt)
ezplot(yt,[y0,y0+2])

Input & Output:


8.
9.
10. a.
10. b.

You might also like