function [dx]=stiffnl(x,t)
  dx(1)=-0.013*x(1)-1000*x(1)*x(3);
  dx(2)=-2500*x(2)*x(3);
  dx(3)=-0.013*x(1)-1000*x(1)*x(3)-2500*x(2)*x(3);
endfunction  
