epsrhs=1e-13;
%epsrhs=0;
digits=14;
 [normpf, normdf, phi, X, y, iter, primalres] = n20(epsrhs,digits);

