fortranプログラムについての質問です。

質問

fortranプログラムについての次の問題について質問させてください。

ｉ＝１、２、・・・、５０に対し（i,a500(i)）(i,a2000(i))を次の演算に従い計算せよという問題です。わかりにくいですが
５００と２０００、さらに次の式のaのあとのnやn+1,n-1などは添字です。

an+1(i)=0.5{an(i+1)+an(i-1)-2.0d0an(i)}+an(i)　　(i=2,3,・・・49),　(n=1,2,・・・,500,・・・,2000)
ただしa1(i)=0.0(i=1,2,・・・,49), a1(50)=1.0, an(1)=0.0, an(50)=1.0 (n=1,2,3,・・・,500,・・・,2000)　である。

というものです。自分は添字のｎなどについてはプログラム上では無視し、一次元配列ａとｂを宣言して、「ＤＯループも用いて
ｉ＝１～４９を回してまずa2（これがb）を求め、そのa2（すなわち）を元の漸化式にの右辺に結果を入れ、それを５００回と２０００回
別々に回す」という操作を考えたのですがうまくいきません。i-1やi+1をどう扱うかや、そもそもこの考え方は合っているのか、
どうするほうが良いのかなどについてアドバイスをいただきたいです。よろしくお願いします。

回答No.2

Tacosan · ベストアンサー

「do j=1,500→enddoをiをはさんだりしてやったみました」というのが
わからないんですが, どうやったんでしょうか? 単純に
do j = 1, 500
a(1) = 0
a(50)=1
do i = 2, 49
b(i) = 式の右辺
enddo
a=b
enddo
で動かない?

補足

説明不足でした。すみません。色々とdo j=1,500・・・enddoの配置を変えたりしながらやってみたということです。
Tacosanさんの書いてくださったプログラムも書いてみたのですがどうもうまくいきませんでした・・・。

回答No.3

Tacosan

だからさぁ,
・どんなプログラムを書いたのか
・どんな結果が得られたのか
・どのような結果が期待されるのか
くらいは書こうよ. 「do～enddo の配置を変えたりしながらやってみた」っていってもどのようなプログラムを使ったのか全くわからんし, 「どうもうまくいきませんでした」といっても何がどう「どうもうまくいかないのか」がわからんの.
もちろん「回答をもらうつもりはなくただ単に愚痴を零しているだけ」ならいいけど, それはこんなところでしないでもらいたい.

補足

do j = 1, 500
a(1) = 0
a(50)=1
do i = 2, 49
b(i) = 式の右辺
enddo
a=b
enddo
とか
a(1) = 0
a(50)=1
do i = 2, 49
b(i) = 式の右辺
enddo
do j = 1, 500
a=b
enddo

とか思いつく限りの組み合わせでＤＯループを回したということです。なんでそこまで怒るのかよく分かりませんが、お返事ありがとうございました。
自力でうまく回せました。

回答No.1

Tacosan

考え方はそれであってると思う.
「うまくいかない」とは, どのようなプログラムでどううまくいかないと言っているんでしょうか?

補足

ＤＯループの配置（？）などがよくわかりません。自分は、初期値をプログラムで書いた後、

do i=2,49
b(i)=上記の演算の右辺
enddo
としたあと、
a(i)=b(i)　としてまた同じ演算を５００回（またそれと別に２０００回）行うことを考えているのですが、どう書けばいいのか思い
つきません・・・。do j=1,500→enddoをiをはさんだりしてやったみましたがそれらしき値が得られなかったということです。