function f=fdjur(t,y)
   f=[2*y(1)-0.01*y(1)*y(2)   -y(2)+0.01*y(1)*y(2)]';


