function x = modhs52(p)
  n=4;

  x(1)=4.0*p(1)-p(2);
  x(2)=p(2)+p(3)-2.0;
  x(3)=p(4)-1.0;
  x(4)=p(5)-1.0;