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
