您所在位置：主题：帮忙看看这个三对角矩阵的程序有错么？已解决！
LELEMA
[专家分：0] 发布于 2011-12-22 11:48:00
subroutine sdj2(n,a,b,c,f)
doubleprecision a(1:n),b(1:n),c(1:n),f(1:n)
doubleprecision W
f(1)=f(1)/b(1)
w=b(1)
do k=1,n-1
b(k)=c(k)/w
w=b(k+1)-a(k+1)*b(k)
f(k+1)=(f(k+1)-a(k+1)*f(k))/w
end do
do k=n-1,1,-1
f(k)=f(k)-b(K)*f(K+1)
END DO
RETURN
END
矩阵的形式在附件里！
doubleprecision a(1:n),b(1:n),c(1:n),f(1:n)
doubleprecision W
f(1)=f(1)/b(1)
w=b(1)
do k=1,n-1
b(k)=c(k)/w
w=b(k+1)-a(k+1)*b(k)
f(k+1)=(f(k+1)-a(k+1)*f(k))/w
end do
do k=n-1,1,-1
f(k)=f(k)-b(K)*f(K+1)
END DO
RETURN
END
矩阵的形式在附件里！