jacmodhs52.m

function jac = modhs52_jac(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];