% odesin, ode23 resp ode45 på y'=sin(ty), y(0)=8
  clear, clf
  tslut=4.8;  y0=8;
  [T,Y]=ode23('fsin',[0 tslut],y0);  n=length(T), yslut=Y(n)
  subplot(1,2,1), plot(T,Y, tslut,yslut,'*'), axis([0 tslut 7.3 8.3])
  title('ode23'), text(0.5,7.5,['yslut=' num2str(yslut)])

  [T,Y]=ode45('fsin',[0 tslut],y0);  n=length(T), yslut=Y(n)
  subplot(1,2,2), plot(T,Y, tslut, yslut,'*'), axis([0 tslut 7.3 8.3])
  title('ode45'), text(0.5,7.5,['yslut=' num2str(yslut)]), hold on

  tol=odeset('RelTol',1e-8);
  [T,Y]=ode45('fsin',[0 tslut],y0,tol);  n=length(T), yslut=Y(n)
  plot(T,Y,'r--', tslut,yslut,'rd')
  text(0.5,7.4,['Reltol=1e-8. yslut=' num2str(yslut)])



