TP1 - Partie 1 - Aides : 4. 1

TP1 - Partie 1 - Aides : Etape 4 : résolution : 4. 1 - Résolution

4. 1 - Résolution

Sous-programme de résolution :

Ecriture d’un sous-programme devant résoudre le problème pour tout

$nx$ : nombre de points $x_ j$ du domaine,

$x$ , tableau des abscisses $x_ j$ , $j=1,nx$ ,

$k_ c$ , tableau des coefficients $kc$ donnés par élément,

$cldi$ , tableau décrivant les conditions de Dirichlet,

fournis par la liste d’appel du sous-programme,
la sortie étant le vecteur $T$ , solution de $KT=F$ .

structure :

- initialisation de $F$ à zero,

- construction de $K$ ,

- modification de $K$ pour tenir compte des CL de dirichlet,

- modification de $F$ pour tenir compte des CL de dirichlet,

- résolution du système $KT=F$ ,

sous scilab : écrire, en complétant les zones grises, dans le fichier
chaleur1DP1_plein.sci (pour un stockage plein) .

$\includegraphics[width=7cm]{./prog-images/chaleur1DP1dipl.png}$

Besoin d’aide pour remplir les zones grises ? pour vous aider, voir le QCM suivant :

Tester ces sous programmes en ajoutant au programme principal précédent

...

exec(’chaleur1D1P_plein.sci’,-1) ;

...

T=chaleur1D1P_di_plein(nx,x,kc,cldi) ; T

On doit obtenir $T=\left(\begin{array}{c}0\\ 2.5\\ 5.\\ 7.5\\ 10\end{array} \right)$ .

si cela n’est pas le cas, chercher l’erreur ...