if true
clear all
tic;
h=0.1;
format short g
[t,ya]=euler(@funObj,0,1,h,1);
[t,yc]=heun(@funObj,0,1,h,1);
[t,yd]=runge4(@funObj,0,1,h,1);
[t,zf]=ode45(@funObj,tspan,1);
dim = length(t);
yx=yd;
disp('------METODO EULER, EULER MOD., HEUN , RUNGE-KUTTA, ODE23 , ODE45 , ODE113 , SOL. EXATA------');
disp('t(i) ye(i) yem(i) yh(i) yk(i) y23(i) y45(i) y113(i) sol. exacta');
disp('------------------------------------------------------------------------------------------------');
disp([t(1:dim), ya(1:dim) , yb(1:dim) , yc(1:dim) , yd(1:dim) , y23(1:dim) , y45(1:dim) , y113(1:dim) , yx(1:dim)])
disp('------------------------------------------------------------------------------------------------');
end