442
Appendix A. Numerical Miscellany
c Forward elimination sweep
q(l)=-c(l)/b(l)
f(l)= f(l)/b(l)
s(l)=-a(l)/b(l)
do j=2,jmx
p=l.O/(b(j)+a(j)*q(j-l»
q(j)=-c(j)*p
f(j)=(f(j)-a(j)*f(j -l»*p
s(j)= -a(j)*s(j -l)*p
end do
c Backward pass
q(jmx)=O .O
s (jmx)=l. 0
do j=jmx-l,l,-l
s(j)=S(j)+q(j)*s(j+l)
q(j)=f(j)+q(j)*q(j+l)
end do
c Final pass
f(jmx)=( fmx -c(jmx)*q(l)-a(jmx)*q(jmx -l) )/
&
( c(jmx)*s(l)+a(jmx)*s(jmx-l)+b(jmx)
do j = 1 , j mx- 1
f(j)=f(jmx)*s(j)+q(j)
end do
return
end