// Résolution de l'équation AX=B

A=[1 4 ; 1 0 ; 1 2 ; 1 1] ; B=[1 ; 2 ; 3 ; 4] // par exemple

[n,p]=size(A) ; r=rank(A)

if (n==p)&(r==n) then
    X=inv(A)*B
    disp("Une solution unique :") ; disp(X)
else
    [X,kerA]=linsolve(A,-B)
    if X==[] then
        disp("Pas de solution")
    elseif kerA==[] then
        disp("Une solution unique :") ; disp(X)
    else
        disp("Solution particulière :") ; disp(X)
        disp("À une combinaison linéaire près de :") ; disp(kerA)
    end
end