r****y
发帖数: 1437 | 1
Give you a sample I once played with,
关于j 的循环是高斯消去的算法, 很简单.
do 30 l = 1, lmp
cxl = fxl(l) - xkss
b(jmm1) = cxl
do 20 jw = 4, jm
j = jm - jw +2
vort(l, j) = vort(l, j) -vort(l, j+1)/b(j+1)
20 b(j) = cxl - 1.0/b(j+1)
vort(l, 2) = vort(l, 2) / b(2)
do 30 j = 3, jmm1
30 vort(l, j) = (vort(l, j)-vort(l, j-1))/b(j) |
|
