méthode de richardson pour les systèmes linéaires