%Lösning av y'=sin(ty),y(0)=6 med ode45
tslut=2;y0=6;
[T45,Y45]=ode45(@fsin,[0,tslut],y0);
[T23,Y23]=ode23(@fsin,[0,tslut],y0);
plot(T45,Y45,'b',T23,Y23,'r')
title('Lösning med ode45 (blå)och ode23 (röd)')
pause
tol=odeset('RelTol',1e-8);
[T45tol,Y45tol]=ode45(@fsin,[0,tslut],y0,tol);
plot(T45,Y45,'o',T45tol,Y45tol,'.')
title('ode45 standardtolerans (o), tol=1e-8 (.)')
pause
[T23tol,Y23tol]=ode23(@fsin,[0,tslut],y0,tol);
plot(T23,Y23,'ro',T45tol,Y45tol,'r.')
title('ode23 standardtolerans (o), tol=1e-8 (.)')