[houkensjtu]

Maybe you know the sample program:
http://www.students.ncl.ac.uk/f.j.he...webs/codes.php

It seems a very good one but when I tried to raise the grid number from 10*10 to 12*12,it leads to a divergence.The main error happened in TDMA:
//Forward substitution
for (int i=1; i < n; i++) {
gam(i) = c(i-1) / bet;
bet = b(i) - a(i) * gam(i);

if (bet == 0.0)
throw("Solver Error 2 in triag: Zero in denominator \n");
u(i) = ( r(i) - a(i) * u(i-1) ) / bet;
}
I can't figure out why it happen and I need your help!