Solving a system of second order differential equations
Show older comments
I get the following error message when I run the code below (although this error does not always occur; sometimes the plot is generated). I do not see why because matlabFunction (for example, in the line "M = matlabFunction(V,'vars', {'t','Y'})" (source: https://www.mathworks.com/help/symbolic/solve-differential-equation-numerically-1.html)) is a standard name.
Error message
Error using sym/matlabFunction>checkVars (line 223)
Variable names must be valid MATLAB variable names.
Error in sym/matlabFunction (line 158)
vars = checkVars(funvars,opts);
Code
syms y(t) z(t )
Dy = diff(y,t );
Dz = diff(z,t );
eqns = [diff(y,t,2) == (29.4*cos(z)+4*Dz*Dz)*sin(z)/(1+3*sin(z)*sin(z )),...
diff(z,t,2) == -3*(19.62+2*Dz*Dz*cos(z))*sin(z)/(2*(1+3*sin(z)*sin(z )))];
fun = matlabFunction(odeToVectorField(eqns),'Vars',{'t','Y '});
tspan = [0 10]; % Time interval for integration
y0 = [0 1 0 1]; % initial conditions
[t,y] = ode45(fun,tspan,y0 );
plot(t,y)
Accepted Answer
More Answers (0)
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!