function [dg,dh]= grestr(x)

x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
x5=x(5);

dg=[];

dh(1,:)=[(exp(x1)) 0 (exp(x3)) (exp(x4)) 0 ]

dh(2,:)=[(4*exp(x1)) (2*exp(x2)) 0 0 (2*exp(x2)) ]

dh(3,:)= [0 (exp(x2)) (exp(x3)) (2*exp(x4)) 0 ]