% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
% empezar con h=0.1, e.g.
[t,y]=rk4('ejer10',[0,0.96],1,h);
subplot(2,1,1); plot(t,y); xlabel('eje t'); ylabel('eje y'); title('solucion numerica,Rk4')
dt=diff(t); tdt=t(2:length(t));
subplot(2,1,2); plot(tdt,dt,'o');xlabel('t');ylabel('tamaño paso');
pause
%
[t,y]=ode45('ejer10',[0,0.96],1);
subplot(2,1,1); plot(t,y); xlabel('eje t'); ylabel('eje y'); title('solucion numerica,ode45')
dt=diff(t); tdt=t(2:length(t));
subplot(2,1,2); plot(tdt,dt,'o');xlabel('t');ylabel('tamaño paso');
