Next: 6.2.2 Remarque sur cette
Up: 6.2 Résolution d'un système
Previous: 6.2 Résolution d'un système
  Table des matières
  Index
Pour résoudre un système linéaire dont la matrice est carrée, Scilab
utilise une factorisation LU avec pivot partiel suivie de la résolution des
deux systèmes triangulaires. Cependant ceci est rendu transparent pour
l'utilisateur par l'intermédiaire de l'opérateur }, essayer:
-->R=rand (5,5);
-->y=rand(5,1);
-->x=R\y // resolution de Rx=y
x =
! - 0.9448615 !
! 1.8150216 !
! 1.0487447 !
! 0.7257294 !
! - 0.4538296 !
Pour se souvenir de cette instruction, il faut avoir en tête le système
initial Ax=y puis faire comme si on multipliait à gauche par l'inverse de A, (ce que
l'on symbolise par une division à gauche par A) d'òu la synthaxe utilisée
par scilab. En général, il est important de noter qu'il peut y avoir des
erreurs d'arrondi dues à l'arithmétique flottante.
root
2000-02-07