How do I manipulate a function for Euler's method to create a function for an improved Euler's method?
5 views (last 30 days)
Here is the code I have for my Euler's method. It works pretty well. I'm unclear as to what modifications I would make to this so that I may create a function for the improved Euler's method. function [ t,y ] = euler( f,tspan,y0,N ) m=length(y0); t0=tspan(1); tf=tspan(2); h=(tf-t0)./N; t=linspace(t0,tf,N+1); y=zeros(m,N+1); y(:,1) = y0'; for n=1:N y(:,n+1) = y(:,n) + h*f(t(n),y(:,n)); end t=t';y=y'; end