L'équation de la chaleur tombe à zéro

[invite589b6513]

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 !


-----

[Fanch5629]

Bonjour.

L'expression de beta semble incorrecte dans le code ( dx doit être au carré au dénominateur).

@+

[invite589b6513]

Bon sang mais c'est bien sûr !

Comment j'ai pu laisser passer ça !

Tout marche impec' ! Merci mille fois !

JP