H1 = heun(@f3,0,2,1/10,10)
H2 = heun(@f3,0,2,1/10,20)

 

E1 = euler(@f3,0,2,1/10,10)

E2 = euler(@f3,0,2,1/10,20)

 R1 = rk4(@f3,0,2,1/10,10)
R2 =rk4(@f3,0,2,1/10,20)

t = (linspace(0,2,11))'

F = .1*exp(-2.*t)+t.*exp(-2.*t)

ErrorEst1 = [(E1(:,2)-F) (H1(:,2)-F) (R1(:,2)-F)]

t = (linspace(0,2,21))'

F = .1*exp(-2.*t)+t.*exp(-2.*t)

 

ErrorEst2 = [(E2(:,2)-F) (H2(:,2)-F) (R2(:,2)-F)]