采用Thomas方法计算不可压缩Couette流的fortran程序,Calculated using Thomas incompressible Couette flow fortran program

poisson.f90 是矩形区域上泊松方程Dirichlet边值问题差分求解的fortran程序，可选择Jacobi或者Gauss-Seidel迭代子程序，其中Jacobi迭代次数为925，Gauss-Seidel迭代次数为478。结果输出至solution.dat文件 poisson.ncl 是NCL脚本，运行该脚本得到的图为result.pdf thomas.f90 是追赶法解代数方程组的fortran程序。控制分辨率n的不同取值，输出以下文件： thomas(n=3).