主题:[讨论]fortran 有限差分法,热扩撒方程调试出错,无结果
PROGRAM MAIN
IMPLICIT NONE
REAL,DIMENSION(0:10,0:6)::T
INTEGER::I,N
REAL::B,TT,XX,S
READ*,B,TT,XX
S=B*TT/(XX**2)
DO I=0,6
T(10,I)=100
END DO !x=1恒为100度
DO I=0,9
T(I,0)=0
END DO !t=0时其他点出为0度
DO I=0,6
T(0,I)=0
END DO !x=0恒为0度
DO N=0,6
DO I=9,1,-1
T(I,N+1)=T(I,N)-S*(T(I+1,N)-2*T(I,N)+T(I-1,N))
END DO
END DO
WRITE(*,*)T
END PROGRAM
问题是:一根1米长的棒子,一端温度衡为100度,另一端一直未0度。求T(X,T)的值?
TT是时间步长,xx是位置步长。结果编译没有问题,运行出错。求大家分析问题。。谢谢