Scribd
Upload
10 views

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
Available Formats
Download as PDF, TXT or read online on Scribd
10 views

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
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