% Piirtää differentiaaliyhtälön x'(t)=f(t,x(t)) ratkaisuja
% hiirellä annetusta alkuarvosta välillä 0<t<1.
% Keski- tai oikeanpuoleinen hiirennäppäin lopettaa.

button=1; hold on
while button<2 ,
   [x1,x2,button]=ginput(1);
   [t,X]=ode45(f,[0,1],[x1;x2]);
   plot(X(:,1),X(:,2),'r'); end
hold off
