Définir les CL de flux dans le programme principal :
On procède comme pour les conditions de Dirichlet en décrivant les conditions de flux par un tableau clne, de dimension (nclne,2) :
la 1ère colonne de clne contient les numéros des noeuds d’extrémités où le flux est imposé,
la 2ème colonne de clne contient les valeurs imposées,
le nombre de ligne est le nombre des extrémités où le flux est imposé
Ici :
une seule extrémité est concernée :
et
est la valeur imposée à l’extrémité de numéro
sous scilab : qo =1 ; clne = [ nx , qo ] ;
remarque : ne pas oublier d’adapter les CL de Dirichlet.