clear();

exec("TP1.sci");
n= 10;
k=1E1 ;
iter_max=1.3E5;
w=1.0;

[H] = initialiserMatriceH(n)
[D] = initialiseMatD(n,k)
[A,b,xi,x1]=ortho(H,D)
m_min(A,n);
//omega_opti(A,n);
[i,Xk]=Jacobi(A,b,n,iter_max);
//[i,Xk]=Jacobi_relax(A,b,n,iter_max,w);
