FORTRAN、フーリエ展開、出力

Question

矩形波のフーリエ展開の式である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でやっています。

Answer 1

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

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

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

Comments

[Naruto_Netu]

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

[curehoney]

うまくいって良かったです！