MATLAB: I can’t plot 2nd order differential equation

differential equationsodeplotplotting

I can't plot 2nd order differential equation.

syms x(t); Dx = diff(x); Dx2 = diff(x,2); x(t) = dsolve(Dx2+(5*Dx)+(6*x)==100,x(0)==0,Dx(0)==0)

Best Answer

I can!

Try this:

syms x(t);
Dx = diff(x);
Dx2 = diff(x,2);
x(t) = dsolve(Dx2+(5*Dx)+(6*x)==100,x(0)==0,Dx(0)==0)
figure(1)
fplot(x, [0 10])
axis([xlim    0  17])
grid

Related Solutions

MATLAB: How to resolve this

Following code solves the problem. The command t(0)==0 does not make any sense since t is independent variable. Also, you used dsolve command incorrectly. See the following code:

syms x(t) m g l e a
Dx=diff(x,t);
x(t)=dsolve(m*diff(x,2) == -m*g*(x-l)/e,[x(0) == a, Dx(0) == 0])

Then, you can verify the solution according to the initial conditions. It is verified.

MATLAB: Using dsolve to solve a system of differential equations analytically

The differential equation system is nonlinear, so it likely does not have an analytic solution.

Try this instead:

m                     = 3.125*0.001;                   % Mass of the body in slug
Cd                   = 4.71*10^-7;                   
g = 32.2;    
syms x(t) y(t) T Y                                      
Dx = diff(x,1);                                    
D2x = diff(x,2);                                    
Dy = diff(y,1);                                  
D2y = diff(y,2);                                  
[VF,Sbs] = odeToVectorField([-Cd*Dx*sqrt(Dx^2+Dy^2)==m*D2x,-Cd*Dy*sqrt(Dx^2+Dy^2)-m*g==m*D2y],...
                  [x(0)==0  ,Dx(0)== 186*8.8/6*cosd(11.2), y(0) == 0, Dy(0)==186*8.8/6*sind(11.2)]);  
Fcn = matlabFunction(VF, 'Vars',{T,Y});
[T,Y] = ode45(Fcn, [0 10], [0 0 0 1]); 
figure
plot(T, Y)
grid
legend(string(Sbs), 'Location','SW')

Best Answer

Related Solutions

MATLAB: How to resolve this

MATLAB: Using dsolve to solve a system of differential equations analytically

Related Question