Discussion :

[Noni23]

Bonjour,

Dans le PDF ci-dessous se trouve un système d'équations aux dérivées partielles. Il m'est demandé de le résoudre numériquement par la méthode des lignes et c'est ce que j'ai fait. Cependant, je ne sais si c'est correct car lorsque je prends tvec=linspace(0,60,61) par exemple, ça a l'air de fonctionner mais quand je prends tvec = (0,1500,31) , ça ne fonctionne plus. Y a-t-il un problème dans mon raisonnement ? dans le code ? Je ne vois honnêtement pas l'erreur.

Le lien de mes scripts sont des fichiers txt copiés de Octave ( mais il me semble que c'est la même syntaxe de Matlab). Le PDF comprend aussi mon raisonnement quant à la manière dont j'ai attaqué le problème.

Merci d'avance !
PDE (4).pdf
pde (3).txt
main (2).txt

[Invité]

Bonjour,
Erreurs dans la fonction PDE : la dérivée première de U au point i pour un schéma centré est D1(i)=(1/(2*dz))*(U(i+1)-U(i-1)) et sa dérivée seconde est D2(i)=(1/(4*dz*dz))*(U(i+2)-2*U(i)+U(i-2))