function [ y,Z] = inityz(A,L);
% Find an initial Z >0 
n2=size(A,2);
n=sqrt(n2);
m=size(A,1);
y=ones(m,1);
Z=A'*y-reshape(L,n2,1);
Z=reshape(Z,n,n);
Z=Z+abs(min(eig(Z)))*eye(n);
