作者t70349 (JAY)
看板Fortran
标题[请益] 关於TDMA的程式码
时间Fri Aug 26 19:27:11 2011
以下是我的程式码 COMPILE无错误
executing却有问题
请各位高手帮忙
我是新手请大家多指教
program aaa
parameter n=5
dimension a(n),b(n),c(n),q(n),bb(n),x(n)
do 10 i=1,n
a(i)=1.
c(i)=1.
10 continue
b(1)=1
b(2)=2
b(3)=3
b(4)=4
b(5)=5
q(1)=3
q(2)=8
q(3)=15
q(4)=24
q(5)=29
call trdiaq(5,a,b,c,x,q)
write(*,*)x
stop
end
subroutine trdiaq(n,a,b,c,x,q)
dimension a(n),b(n),c(n),q(n),bb(n),x(n)
do 1 i=1,n
bb(i)=b(i)
1 continue
do 2 i=2,n
t=a(i)/bb(i-1)
bb(i)=bb(i)-c(i-1)*t
q(i)=q(i)-q(i-1)*t
2 continue
x(n)=q(n)/bb(n)
do 3 i=1,n-1
j=n-i
x(j)=(q(j)-c(j)*x(j+1))/bb(j)
3 continue
return
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.124.33.174
1F：推 terryys:有什麽问题？起码说明一下吧 08/26 19:57
2F：→ t70349:一个error,但完全不知道错在哪^^ 08/26 20:14
3F：推 terryys:我编译(ifort test.f90)执行没有问题，输出为1 2 3 4 5 08/26 20:36
4F：→ terryys:不知道你编译有没有加上别的档案？ 08/26 20:36
5F：→ t70349:我用compaq怎跑不出来 08/26 20:39
6F：→ t70349:完全没加其他档案 08/26 20:42
7F：→ t70349:terryys谢谢你...确实加上其他档案了 08/26 20:49