clear;

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

for k = 1 : 20
    tmpA = randn(100, m);
    tmpy = tmpA*x + randn(100,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 : 20;
plot(k,err,'r',k,diff,'b');