北航计算流体力学第11课.docVIP

用于求解Laplace方程的网格系统 , , 式中 求解三对角方程的追赶法简介：（以方程组为例） （1） （2） （3） （4） （5） 由式（1）得： （即） 代入式（2）得： （即） 继续上面的工作，可得： 最后求得： 然后回代，依次求出从而完成了求解工作。 以下是一段求解任意三对角方程组的Fortran程序 Subroutine TRIDAG SUBROUTINE TRIDAG (AA,BB,CC,FF,N1,N) C C******************************************************************************************* C SOLUTION OF A TRIDIAGONAL SYSTEM OF N-N1+1 EQUATIONS OF THE FORM C C AA(K)*X(K-1) + BB(K) * X(K) + CC(K) * X(K+1) = FF(K) K=N1, … ,N C C K RANGING FROM N1 TO N C