6.2.1 Méthode de résolution basique

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.