MATLAB Answers

2nd order differential errors

1 view (last 30 days)
John Masters
John Masters on 15 Apr 2021
Edited: John Masters on 16 Apr 2021
Hi,
I am having issues with ODE45,
My function
function dydt=solver(y,t)
Caller
tspan = [0:0.01:3.12];
thrust=6;
g=9.81;
R=1.3;
m12=4.3;

Answers (1)

William Rose
William Rose on 15 Apr 2021
Are you getting an error message? If so, what?
I notice that your dydt is mising a "*" after C.
  3 Comments
William Rose
William Rose on 15 Apr 2021
@John Masters The code works when you fix the (t,y) ordering. By the way, you can specify just the start and end time in [tspan]. YOu don;t have to tell it the time step . That results in a more efficient execution because in that case, ode45() adapts the step size as it goes, as need to achieve a certain level of accuracy. See attached image. When I used you code as is, there were 313 evenly space time steps. When I only specified the start and end
tspan = [0,3.12];
there were 105 unevenly spaced times, with a concentration around the "glitch" (staging, perhaps?). See figure below. Code attached.

Sign in to comment.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!