L'équation de la chaleur tombe à zéro
Bonjour !
J'essaye depuis quelques jours d'approcher l'équation de la chaleur en dimension 1 par différences finies (schéma d'Euler implicite) et la température de la barre de cuivre tombe inéxorablement à 0 partout... sauf un peu au bord... Dans un pdf en pièce jointe je présente mon travail, l'autre pièce jointe est mon programme en fortran 90 qui inverse la matrice à chaque pas de temps.
Je ne sais pas quoi faire, même en mettant la température initiale à 10 elle tombe à 0 toute seule !
-----