clear;

m = 200;
k = 50;
noisePower = 0.001;
x = zeros(200,1);
x(1:k) = 1;
A = [];
err = [];
diff= [];
y   = [];
oldSol = x*0;

for k = 1 : 50
    tmpA = randn(20, m);
    tmpy = tmpA*x + randn(20,1)*sqrt(0.01);
    A = [A; tmpA];
    y = [y; tmpy];
    sol = SolveBP(A, y, m, 100, 0.1, 1e-3);
    diff = [diff, norm(sol-oldSol)];
    err = [err, norm(x-sol)];
    oldSol = sol;
end

k = 1 : 50;
plot(k,err,'r',k,diff,'b');