% file name: gradh.m % problem 39: h Jacobian (Miele et al.) n=4,m=3 function M = gradh(x) M=zeros(4,2); M(1,1)=-3*x(1)^2; M(2,1)=1; M(3,1)=-2*x(3); M(4,1)=0; M(1,2)=2*x(1); M(2,2)=-1; M(3,2)=0; M(4,2)=-2*x(4);