質問: sum + sum0 K1= (pi/2.0d0)*sum enddo write(fo,*) K1 close(fo) end program first_eliptic_integral ### 試したこと nをカウントしたいので、do i=1,n文の外に更にnをカウントするdo文を置いてみたり、if文を使ってn-1を取り出そうとしましたが、出力結果は全部同じ数値になってしまいます。 ### 補足情報(言語/FW/ツール等のバージョンなど) 現在配列は勉強中なので、do文, if文(exit, cycle)のみで教えていただけたらうれしいです。(GOTO文もなしでお願いします.) 参考書籍は数値計算のためのFortran90/95プログラミング入門(牛島)です。 環境はgfortranでコンパイルしています。 ※ここを利用するのは初めてなので、よろしくお願いします。
回答: 今プログラム中の pi/2 * sum0 が第N項にあたるわけだから、それをループの内側に書けばいいのではないかと思う。 write(*, *) i, pi/2.0d0 * sum0