function [x1]=SP_lp_re(K,phi,y,p)
[xfinal,That]=SPnew(K, phi, y);
[x1] = lp_re1(phi,y,p,xfinal);