[chberbelali]

I have just searched the Nusselt number delicate and compared it with those of other literatures,
I found that I have to do the test of stop inside the loop "While",
here is what I have added in the Matlab script:

%************************************************* *
initialize Theta;
Tnew = Theta;
Told = Tnew;

while (norm(Tnew-Told,1)/norm(Tnew),1) <= 1e-6


Told = Tnew;
.....
..... The body of the program.
.....
Tnew = Theta (Calculated);


end (while)

%*************************************************

(it doesn't work!!! it just passes the 1st iteration.)

i'll be grateful for any help !