function [ X] = initx(A,b);
% Find an initial X >0 
n=sqrt(size(A,2));
X=eye(n);
