Fortran Coder

标题: TDMA算法一直改不对不知道是哪里出了问题请求指点 [打印本页]

作者: 纯小白 时间: 2023-3-8 18:01
标题: TDMA算法一直改不对不知道是哪里出了问题请求指点
TDMA
c(1)=c(1)/b(1)
x(1)=x(1)/b(1)
  ！！max=3
do i=2,max
      middle=1.0d0/b(i)-c(i-1)*a(i)
      c(i)=c(i)*middle
      x(i)=(x(i)-x(i-1)*a(i))*middle
end do

do i=max-1, 1,-1
      x(i)=x(i)-c(i)*x(i+1)
end do

   open(unit=10,file='result.txt')
      do i=1,max
         write(10,*) x(i)
      end do

作者: 纯小白 时间: 2023-3-8 18:12
发现少了个括号但是想请问一下这个算法应该怎么验证

作者: 纯小白 时间: 2023-3-8 18:16
加了个括号还是不对唉

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)