function res = performance_vs_epsilon

	figure; hold on;

	n = 10;
	x0 = rand(n,1);
	A = generate_connected(10);

%	epsilons=[0.02 0.05 0.1 0.2 0.3]

	epsilons = [0.2 0.1 0.05 0.025 0.0125 0.0064]
	
	colors = {'r','g','b','k','m'};

	for i=1:numel(epsilons)
		epsilon=epsilons(i);

	 	r = cons00(A,x0, epsilon, 5, 2000);

		p = plot(log_or_small(r.mse, -30), '-');

		set(p, 'color', colors{mod(i, numel(colors))+1});
%		set(p, 'markersize', 2);
		
		pause(0.1)
		
		res{i} = r;
	end


	for i=1:numel(epsilons)
		legends{i} = sprintf('\\epsilon = %f',epsilons(i));
	end
	legend(legends);


function r = log_or_small(k, min)
	for i=1:numel(k)
		r(i) = log(max(k(i),exp(min)));
	end