function jac = jacmodhs52(p)
  m=5;

  jac(1, 1:m)=[4.0, -1.0, 0.0, 0.0, 0.0];
  jac(2, 1:m)=[0.0, 1.0, 1.0, 0.0, 0.0];
  jac(3, 1:m)=[0.0, 0.0, 0.0, 1.0, 0.0];
  jac(4, 1:m)=[0.0, 0.0, 0.0, 0.0, 1.0];