🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

Q&A

解決済

1回答

1146閲覧

FORTRAN、フーリエ展開、出力

Naruto_Netu

総合スコア8

FORTRAN

FORTRAN(フォートラン)は科学時術計算に向いた手続き型プログラミング言語です。 並列計算の最適化が行いやすい特性上、数値予報および気候モデルなどの大規模な計算を行う分野のスーパーコンピュータで使われています。

0グッド

0クリップ

投稿2021/02/23 06:41

矩形波のフーリエ展開の式であるf(x)=4/π*(sinx+sin3x/3+sin5x/5+...+sin13x/13)(0≦x≦2π)
のグラフの値をFORTANで画面上に出力したいのですが、自分の書いているコードだと思った通りに出力することができません。

コードは下のとおりです。
'
INTEGER I,J
REAL*4 X(632),Y1(632),SUM(632,9),PI,

PI=4.0*ATAN(1.0) DO 80 I=1,630 X(I)=0.+FLOAT(I-1)/100 DO 90 J=1,7 SUM(I,J)=SUM(I,J-1)+(4/PI)*(SIN((2*J-1)*X(I))/(2*J-1))

90 CONTINUE
WRITE(,) X(I),SUM(I,2),SUM(I,3),SUM(I,4),SUM(I,5),SUM(I,6),
$ SUM(I,7)
80 CONTINUE```
といった感じです。
このコードをどう改良したらよいのか教えていただきたいです。もしくはこのやりたいことを書く時のヒントやこのプログラムとは違うソースコードを明示していただきたいです。
ちなみにFORTRAN77でやっています。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

fortran

1DO 90 J=1,7 2 SUM(I,J)=SUM(I,J-1)+云々 3・・・・

のところで J=1 の時、J-1 が 0 となって、配列 SUM の宣言長からはみ出しているのが問題の気がします。
ただどういう問題が生じているのか書いていないので、気がするだけです。

>自分の書いているコードだと思った通りに出力することができません。
この部分、どう思っていて、どうなったのかを質問に書かないと、誰も答えようがないです。

投稿2021/02/24 02:03

編集2021/02/24 02:04
curehoney

総合スコア249

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

Naruto_Netu

2021/02/24 02:35

情報が足りなくてすみません。正直配列に関してよく分かっておらず、後のものに前のものを代入するという考えのもと、しました。そして、curehoneyさんの回答したことをもとに改良したら、思った通りの値が出ました。ありがとうございます。
curehoney

2021/02/24 11:57

うまくいって良かったです!
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問