>>697
説明不足ですいません・・・ソースというのはプログラムのことですよね。
こんな感じで書いてみました。

program
implicit none
integer :: i,j
real(8) :: a(50),b(50)
do i=1,49
a(i)=0.0d0
b(i)=0.0d0
enddo
do j = 1, 500
a(1) = 0.0d0
a(50)=1.0d0
do i = 2, 49
b(i) = 式の右辺
enddo
a(j)=b(j)
enddo
do i=1,50
write(*,*),a(i)
enddo
end program

これでいけるかと思ったのですが、出力するとa(48)が0.5になりあとは
全てゼロになるということになってしまいます。これは一次元熱伝導方程式
についての温度分布の時間変化を表したものらしく、だとするとこのような
値になるはずがないのに、という話です。