function [h] = hess_my(ip,iq,ry,rz,ryz,rzy) 
      
for j=1:ip 
      	for k=1:ip 
		h(j,k)=2.*ry(abs(j-k)+1);
	end
end
for j=ip+1:ip+iq 
	for k=ip+1:ip+iq 
		h(j,k)=2.*rz(abs(j-k)+1);
	end
end
for j=1:ip 
	for k=ip+1:ip+iq 
		if j >= (k-ip)
			h(j,k)=-2.*ryz(j-(k-ip)+1);
		end
		if j < (k-ip)
			h(j,k)=-2.*rzy((k-ip)-j+1);
		end
	end
end
for j=ip+1:ip+iq 
	for k=1:ip 
		h(j,k)=h(k,j);
	end
end