Error Message: Improved Euler Step method code

Question

Catherine Zdunek 2021-2-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/755534-error-message-improved-euler-step-method-code

评论： Catherine Zdunek 2021-2-25

Hi everyone! I am working on a code that solves an equation using the improved euler method. I have all my variables in place, but whenever I try to run the code I get multiple error messages which don't tell me explicity what the error is. I'm not quite sure what I could be doing wrong, but I can put my code here if anyone has any suggestions. The error message is occurs at the dV2=... line in the rhs function, the y_low=... line in the improved_euler_step function, and in the first line after the while loop.

The first function is the function I am trying to run. I uses function handles to call the other two functions. The second is the improved euler step function. The third is a function that solves the rhs of the equation provided to be used for the improved euler method. I am trying to run with the given parameters defined in the param statement and over a timespan of 0-100 (which is input as tspan). The initial conditions for y are y(0)=0, y'(0)=0. Thank you all so much for your help!!

[t_out, y_out] = ode(@rhs, @improved_euler_step, [0 100], [0;0], 0.5);
param=[0.5 0.1];
y=[V;dV];
yp=[dV;dV2];
function [t_out, y_out] = ode(rhs, method, tspan, y_init, h)
     t=tspan(1)
     t_out=[t];
     y=y_init;
     y_out=[y_init(1), y_init(2)];
     while t<tspan(2)
         [t_2,y_low,y_high]= method(rhs,t,y,h);
         y=y_high;
         t_out=[t_out;t_2];
         y_out=[y_out; y_high];
      end
  end
function [t_2, y_low, y_high] = improved_euler_step(rhs, t_2, y, h)
% calculate the next time by increasing the current time, t_in by one time step h
t_2 = t_2 + h;
% calculate an approximation of y at t_out using Euler
y_low = y + h*rhs(t_2,y);
% calculate an approximation of y at t_out using improved Euler.
y_high = y + h*(rhs(t_2,y) + rhs(t_2, y_low))/2;
end
function yp = rhs(t, y, param)
t=t;
V=y(1);
dV=y(2);
dV2= -dV-param(1)*V +cos(param(2)*t);
yp= [dV;dV2]
end