Неявная разностная схема, метод прогонки

for(stept = 0 ; stept <= tend ; stept+=dt)
                 {          
 
                     alfa[0] = 0;                
                     if(stept==0)
                     beta[0] =  Tg;
                     else
                     beta[0] = T[0];                     
                    
                     for(i = 1,steph=dh;i<=N;i++,steph+=dh)
                     {
 
                         A=alf*stept/pow(steph,2);
                         B=A;
                         C=1+2*A;
                         if(stept==0)
                         F=Tg-exp(-steph);
                         else
                         F=-T[i];
 
                         gamma = C-B*alfa[i-1];
                         alfa[i]=A/gamma;
                         beta[i]=(B*beta[i-1]-F)/gamma;
 
                     }
                     
                     if(stept==0)
                         T[N]=exp(-h)-Tg;
                     else
                         T[N]=T[N-1]-alf*dh/lambda*(T[N-1]-Tg);
                     
                     for(i=N-1;i>=0;i--)
                     {
                         T[i]=alfa[i]*T[i+1]+beta[i];
 
                     }
                 }