clear
N=128; 
K=20;
x=[rand(1,K) zeros(1,N-K)]'; 
x=x(randperm(128));    
m=floor(3.5*K);
Phi=randn(m,N)/sqrt(m);aPhi=Phi';          
y=Phi*x;
figure; plot(x);
nx=cspocs_l1(Phi,aPhi,norm(x,1),y,10000);      
% Stop at n=2157, score=9.987234e-11 ...
% Final score = 9.987234e-11 ...
hold on; plot(nx,'ro');