r****y
发帖数: 1437 | 1
图书馆啊, 网上以前有盗板的CD, 上面该有源程序.
就是没程序也很简单, 自己手写该不难, 大概思想是消元, 从
第一个方程消起, 到最后一个, 就剩下一个未知数, 解出, 然后
步步回溯, 搞定! | r****y
发帖数: 1437 | 2
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) |
|
