% ex1eul,  y'=sin(ty), y(0)=6 med Eulers metod
  clear, clf
  tslut=2; 
  disp('     h        yslut ')                 
  t=0; y=6;  T=t; Y=y; h=0.2;         % stega med h=0.2 
  while t<tslut-h/2
    f=sin(t*y); 
    y=y+h*f;  t=t+h;  T=[T; t]; Y=[Y; y]; 
  end
  disp([h y])
  plot(T,Y,'r', T,Y,'o')
  axis([0 tslut  5.8  6.4]), hold on
  for n=2:4                           % räkna om allt med halverat steg
    t=0; y=6;  T=t; Y=y; h=h/2; 
    while t<tslut-h/2
      f=sin(t*y); 
      y=y+h*f;  t=t+h;  T=[T; t]; Y=[Y; y]; 
    end
    disp([h y])
    if n==2, plot(T,Y,'g--', T,Y,'x')
    elseif n==3, plot(T,Y,'c-.'), else plot(T,Y,':'), end
  end
  text(0.5,6.36,'Eulers metod på  dy/dt=sin(ty),  y(0)=6')



