fortran90 の簡単なコード

質問

初めて質問させていただきます。

fortran を使い15を任意の数で割るプログラムを作成しようとしているのですが
実行時に n を read し結果を書き出すとn+1個の数値が書き出されてしまいます。かつ、n+1個目の数が全く予想しない数値になってしまいます…

コンパイルエラーは発生しません。
初心者なりにひとつひとつ精査したのですが原因がわかりません。
お助けください、よろしくお願いいたします。

画像：実行時n=10としたときの例

PROGRAM work6

IMPLICIT NONE

REAL , ALLOCATABLE :: x(:)
REAL :: dx
INTEGER :: n

PRINT *,' How many numbers do you separate it? '

READ(*,*) n

ALLOCATE (x(n))

dx=15.0/REAL(n)

do n=1,n
x(n)=dx*REAL(n)
enddo

　PRINT *,x(1:n)

　READ (*,*)

ENDPROGRAM work6

回答No.1

asuncion · ベストアンサー

＞do n=1,n

ループの制御変数がnでいいんでしょうか？
nとは別の変数を使って、
例えば

do i = 1, n

とかいう風にしないとまずいのでは？

でもって、このdoループを抜けたときのnの値は
もともとのnよりも1大きいはずなので、

＞PRINT *,x(1:n)

このPRINT文で最後に出力する
x(n)
は、ALLOCATE文で定義してあった
x(n)
の領域の範囲を1個だけ超えているはずです。
訳のわからない値を出力するのは必然かと。

お礼

ありがとうございました、これで詰まっていた課題を次の段階に進めさせられます！

補足

早速のご回答ありがとうございます！やはり別の変数にしなければいけないのですね、いろいろな可能性を考えて「同じ変数にしないといけないのか？」というふとした疑問でそのように変えていました。

”doループを抜けたときのnの値はもともとのnよりも1大きい”
どうやらここに原因がありそうですね。
doループというのは do i=1.n と書けば　
”iを1からnまで一つずつ大きくしていく”
という意味ではないのですか？どうしてもとのnより1大きくなるのでしょうか。

回答No.3

Hayashi_Trek

＞どうしてもとのnより1大きくなるのでしょうか。
doループを続けるか、抜けるかを判断する箇所が、ループの先頭にあるからです。
つまり「 do i=1.n」の場合では
iに、1を加算して、nを超えたから、ループを抜けるのです。
したがって、ループを抜けた時点では、iはn+1になっています。

お礼

なるほど！！
DO文の本質がやっとわかりました！！
ご丁寧にありがとうございました

回答No.2

asuncion

というわけで、

INTEGER :: n, i

do i = 1, n
x(i) = dx * REAL(i)
enddo

と書けばうまくいくのではないか、という、
Fortran77までしか知らない素人のたわごとでした。