function jac = jacmodhs38(p)
  m=4;

  jac(1, 1:m)=[-20.0*p(1), 10.0, 0.0, 0.0];
  jac(2, 1:m)=[-1.0, 0.0, 0.0, 0.0];
  jac(3, 1:m)=[0.0, 0.0, -2.0*sqrt(90.0)*p(3), sqrt(90.0)];
  jac(4, 1:m)=[0.0, 0.0, -1.0, 0.0];
  jac(5, 1:m)=[0.0, sqrt(10.1), 0.0, 0.0];
  jac(6, 1:m)=[0.0, 0.0, 0.0, sqrt(10.1)];
