Salut
Le code pour le schéma implicite a l'air correct, il marche bien? Désolé je l'ai pas testé, j'ai pas scilab.
S'il marche pas je peux faire une ou deux adaptations et le tester avec autre chose, dans ce cas dis-moi.
Pour le schéma implicite, au lieu de discrétiser la dérivée spatiale au temps n, tu discrétises la dérivée spatiale au temps n+1.
Au lieu d'avoir U(n+1)=BU(n) ou B=I+M comme tu avais pour le schéma explicite, si je me trompe pas tu auras B'U(n+1)=U(n) ou B'=I-M. J'ai trouvé un lien sur l'équation de la chaleur où ça parle de schéma explicite/implicite si tu veux jeter un oeil:
http://www-cast3m.cea.fr/cast3m/html/Dabbene/R03035/node19.html
Contrairement au schéma explicite qui est stable si 2adt<dx², le schéma implicite est inconditionnellement stable.
Le code pour la méthode implicite devrait être un peu plus compliqué puisque tu as à inverser une matrice. Mais c'est pas dramatique non plus. Par contre le code sera plus lent.
A plus
--
Some folks are born made to wave the flag, ooo, they're red, white and blue.
And when the band plays "Hail to the Chief", ooo, they point the cannon at you, y'all!