disp(' '), disp('question 2)b)')
n=100;
U=zeros(n,n);
for i=1:n, for j=1:n, U(i,j)=sin(i)*sin(j); end; end;
for i=1:n, U(i,i)=n; end;
B=ones(n,1);

alphaoptimal=2/(min(eig(U))+max(eig(U)));
[sol_gpc,nb_iter_gpc]=gpc(U,B,alphaoptimal,zeros(n,1),1000,0.001);
[sol_gpv,nb_iter_gpv]=gpv(U,B,zeros(n,1),1000,0.001);
fprintf(' \nLa solution de UX=B est obtenue\n')
fprintf('par gpc optimal au bout de %d itérations avec une erreur de %.1e\n'...
,nb_iter_gpc,norm(sol_gpc-U\B))
fprintf('et par gpv au bout de %d itérations avec une erreur de %.1e.\n\n'...
,nb_iter_gpv,norm(sol_gpv-U\B))