%Example to compare the error of Euler's Method to RK4.
%  Uses euler.m and solve RK4.m

f=@(t,y)(t*y+t.^3);
interval=[0,1];
starty=1;
ExactAnswer=3*(exp(1/2)-1);

h=[0.2 0.1 0.05 0.025 0.0125 0.00625 0.003125 0.0015625];
N=length(h);
for j=1:N
  [t1,y1]=eulerBasic(f,interval,starty,h(j));
  [t2,y2]=solveRK4(f,interval,starty,h(j));
  Error1(j)=abs(y1(end)-ExactAnswer);
  Error2(j)=abs(y2(end)-ExactAnswer);
end
[Error1' Error2']