Programmation à l'aide du logiciel Matlab

Le logiciel Matlab est un outil très performant pour la résolution de grands systèmes linéaires. L'utilisation de Matlab permet de faire des expériences numériques sans avoir à programmer en C ou Fortran.

Matlab contient une bibliothèque numérique très complète pour la résolution de grands systèmes linéaires creux (tapez help sparfun pour avoir des précisions). Par exemple, vous pouvez effectuer l'élimination de Gauss d'une matrice tridiagonale en procédant de la manière suivante.

N = 10;
a = 2;
c = 1;
A = sparse(2:N,1:N-1,-c,N,N) + sparse(1:N,1:N,a,N,N) + sparse(1:N-1,2:N,-c,N,N);
full(A)
b = ones([N,1]);
x=A\b

Le résultat est :

x =

    5.0000
    9.0000
   12.0000
   14.0000
   15.0000
   15.0000
   14.0000
   12.0000
    9.0000
    5.0000