function [ X, y, Z] = psd_strt( C, a, store, nonz);
[n, n1] = size( C);
X=eye(n);
y = zeros(n,1);
y(n)=max(sum(abs(C)))*1.1;
Z = op_at( y, store, nonz, n) - C;
