function udot=lorenz(t,u,flag,a,b,c)
 % odefile for lorenz system
   udot(1)=a*(u(2)-u(1));
   udot(2)=u(1)*(b-u(3))-u(2);
   udot(3)=u(1)*u(2)-c*u(3);
   udot=udot';
 return