How do I manipulate a function for Euler's method to create a function for an improved Euler's method?

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

Answer 1

[Dave Behera]

The following link describes the Euler method and the improved Euler method:

http://www.math.ubc.ca/~israel/m215/impeuler/impeuler.html

You will need to modify the line inside the for loop for the new function.