5-求解三实对角线方程组-p66.txt

来自「地球物理大师Claebout的几个有用的程序」· 文本 代码 · 共 21 行

#公式见P67
# real tridiagonal equation solver

subroutine rtris(n,endl,a,b,c,endr,d,q,e,f)
real q(n),d(n),f(n),e(n),a,b,c,den,endl,endr
e(1)=-a/endl; f(1)=d(1)/endl

do i=2,n-1
{
	den=b+c*e(i-1)
	e(i)=-a/den
	f(i)=(d(i)-c*f(i-1))/den
}
q(n)=(d(n)-c*f(n-1))/(endr+c*e(n-1))

do i=n-1,1,-1
	q(i)=e(i)*q(i+1)+f(i)

return
end