echo on
t=clock;
[T2,W2ab2]=ab_2step(Fun,t0,tlast,w0,N2);
tab2=etime(clock,t);t=clock;
[T2,W2ab3]=ab_3step(Fun,t0,tlast,w0,N2);
tab3=etime(clock,t);t=clock;
[T2,W2ab4]=ab_4step(Fun,t0,tlast,w0,N2);
tab4=etime(clock,t);t=clock;
[T2,W2ab5]=ab_5step(Fun,t0,tlast,w0,N2);
tab5=etime(clock,t);t=clock;
[T2,W2pc54]=precor54(Fun,t0,tlast,w0,N2);
tpc54=etime(clock,t);t=clock;
[T2,W2rkmid]=rkmid_2(Fun,t0,tlast,w0,N2);
tmid=etime(clock,t);t=clock;
[T2,W2mode]=modeul_2(Fun,t0,tlast,w0,N2);
tmode=etime(clock,t);t=clock;
[T2,W2heun]=heun_2(Fun,t0,tlast,w0,N2);
theun=etime(clock,t);t=clock;
[T2,W2rk4]=rk4(Fun,t0,tlast,w0,N2);
trk4=etime(clock,t);

hmax=.1;hmin=.001;TOL=.001;
t=clock;
[Tv1 Wrkf Hrkf]=rkf53_2(Fun,t0,tlast,w0,hmax,hmin,TOL);
trkf=etime(clock,t);
t=clock;
[Tv2 Wpc55 Hpc55]=precor55(Fun,t0,tlast,w0,hmax,hmin,TOL);
tpc55=etime(clock,t);
