clc ;clear ;
N=input('donner la taille de la matrice ' ) ;
A=input('donner la matrice a:') ;
b=input('donner la vecteur b:' ) ; 

for k=1:N -1 
    for m=k+1:N
        f=A(m,k)/A(k,k) ;
        for n=k+1:N
            A(m,n) =A(m,n) -f*A(k,n) ;
        end ;
        b(m) = b(m) -f*b(k) ;
    end;
    end ;
    x=zeros(1,N) ;
    for m=N:-1:1
        somme =0 ;
        for n=m+1:N
            somme = somme + A(m,n)*x(n) ;
        end;
        x(m)=(b(m) - somme)/A(m,m) ;
    end;
    x