function dx=hires(x,t)
 dx=zeros(x);
      dx (1) = -1.71D0*x(1) + 0.43D0*x(2) + 8.32D0*x(3) + 0.0007D0
      dx (2) = 1.71D0*x(1) - 8.75D0*x(2)
      dx (3) = -10.03D0*x(3) + 0.43D0*x(4) + 0.035D0*x(5)
      dx (4) = 8.32D0*x(2) + 1.71D0*x(3) - 1.12D0*x(4)
      dx (5) = -1.745D0*x(5) + 0.43D0*x(6) + 0.43D0*x(7)
      dx (6) = -280.0D0*x(6)*x(8) + 0.69D0*x(4) + 1.71D0*x(5) - 0.43D0*x(6) + 0.69D0*x(7)
      dx (7) = 280.0D0*x(6)*x(8) - 1.81D0*x(7)
      dx (8) = -dx (7)
 endfunction     
