矩形波のフーリエ展開の式である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でやっています。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2021/02/24 02:35
2021/02/24 11:57