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);

exec gpc.sci; exec gpv.sci;
alphaoptimal=2/(min(spec(U))+max(spec(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);
printf(' \nLa solution de UX=B est obtenue\n..
par gpc optimal au bout de %d itérations avec une erreur de %.1e\n..
et par gpv au bout de %d itérations avec une erreur de %.1e.\n'..
,nb_iter_gpc,norm(sol_gpc-U\B),nb_iter_gpv,norm(sol_gpv-U\B))