La condition initiale u(0)=1=U1 est prise en compte simplement en réécrivant le système comme suit :
Soit :
Cette étape est ainsi résumée : F=F-U1*K( :,1)
On supprime ensuite la première ligne de K et celle du second membre car la valeur de U1 est connue. Vu que u(π)=0=U6 , la dernière colonne de K n'a pas lieu d'être. Et comme U6 est également connue, la dernière ligne de K et celle du second membre devront être supprimées elles aussi, d'où :
Le système à résoudre au final se présente ainsi :
La resolution conduit à
La représentation graphique de cette solution figure ci-contre : |
|