
A=2*rand(3,3)-ones(3,3)+10*diag([1,1,1])
b=2*rand(3,1)-ones(3,1)

ti=timer()

solEvo = algoEvo(A,b,0.01,4,0.7)

disp(timer()-ti)

solExact = linsolve(A,-b)


S=algoEvo2(A,b,0.01,4,0.7)
E=[]

n=size(S)

for i=1:n(2)
   E=[E,critere(A,b,S(:,i))] 
end

plot2d(E)

ti = timer()

SolJacobi=jacobi(A,b,0.1)

disp(timer()-ti)
