FORTRANのDO WHILE問題です。
次の漸化式で与えられる数列{Xk}が、
|Xk - Xk-1|< 0.00001
となる時のXkの値を求めよ。
Xk+1 = Xk -
(Xk)^3 - 6(Xk)^2 + 11Xk - 6
/ 3(Xk)^2 - 12Xk + 11
(Xk)^3はXkの3乗です。
私の回答は、
implicit none
integer :: k
real,dimension(100000000) :: x
k = 1
x(0) = 0
do while ( abs(x(k)-x(k-1)) >= 0.00001 )
x(k+1)=x(k)-(x(k)**3-6*(x(k)**2)+(11*x(k))-6)&
/(3*(x(k)**2)-12*x(k)+11)
x(k) = x(k+1)
k = k +1
end do
write(*,*) x(k)
stop
というものなんですが、作動しません。
どこが違うのか教えてください!!お願いします!!
No.3ベストアンサー
- 回答日時:
うーん、新しいFORTRANは読みにくい
簡単に指摘すると
・FORTRANでは、配列の添え字は何も書かないと1からです(0:1000とかで指定すればいいのかな)。なので、最初のx(0)への代入はバウンダリエラー
・初期化されていないx(1)を、最初に比較としてdo whileで使っているので、FORTRANの仕様上、動作は不定。実際には、大概のプログラムローダーでは0クリアされるが、その場合、0と0の比較でループは回らない
・偶々一回目のループが回ったとして、その場合、x(1)は不定なので、x(2)は不定の値を元に計算することになり、計算する度に結果が変わる。
なお、配列で処理するのもあまり感心できません、100000000までに必ず収束するかは場合によって変わりますからね。
というわけで、FORTRAN77で書きますが
integer k
real*8 x0, x1
k = 1
x0 = 0
x1 = xfunc(x0)
do while ( abs(x1-x0) >= 0.00001 )
x0=x1
x1=xfunc(x0)
k = k +1
end do
write(*,*) k,x1
stop
real*8 function xfunc(x)
xfunc=x-x**3-6*(x**2)+11*x-6/(3*(x**2)-12*x+11
return
最近、FORTRANさわってないので、シンタックスは危ないと思う
No.2
- 回答日時:
かなり昔の経験なので自信はありませんが・・・
初期値はx(1)であり、x(0)はおかしいと思います。
判別の不等号が逆では?
計算式のミスなどにより永久ループとなるのを防止するために、do whileの判別式には直接の計算結果を使わない方が良いでしょう。
k = 1
x(1) = 0
do while ( k > 10000000 )
x(k+1)=x(k)-(x(k)**3-6*(x(k)**2)+(11*x(k))-6)&
/(3*(x(k)**2)-12*x(k)+11)
if ( abs(x(k)-x(k-1)) <= 0.00001 ) k=10000001
x(k) = x(k+1)
k = k +1
end do
間違っていたら、ごめんなさいね
そうなんですかっ!!
アドバイスを参考にもう一度考え直してみます☆
補足にも書きましたが、
問題文が十分でなかったことをお詫びします(><)
ありがとうございました!!!
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 リーマン和 1 2022/12/01 13:32
- 工学 Pythonの3Dグラフ表示に関する質問です。 1 2022/12/06 15:03
- 数学 外測度について 2 2023/01/12 07:28
- BTOパソコン BTOPCのレベルアップ 8 2023/02/15 20:06
- 数学 データの分析と標準偏差 5 2022/03/25 12:55
- Visual Basic(VBA) Excelで下記のようにマクロを作ったところ、一回目は実行できたのですが、二回目以降「実行時エラー1 1 2022/03/25 08:08
- Visual Basic(VBA) ExcelVBA No.を自動連番で設定をしながらデータ入力をしたい 2 2022/08/03 18:19
- UNIX・Linux 次の要件を満たすにはどのように修正したらよろしいでしょうか 1 2022/11/24 20:57
- Visual Basic(VBA) Vbaで数式をポーランド記法に変換するコードを作って実行しようとするとフリーズします。 1 2022/05/24 17:53
- 数学 「FFTの基本は、DFTはサンプル数Nが偶数なら 2つのDFTに分解できるということ。 分解するとD 3 2022/03/31 21:01
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
フローチャートで。
-
フローチャートの菱形が狭い。。。
-
fortran errorについて
-
フローチャート以外の設計方法
-
65536は2の何乗なのでしょうか?
-
正しい五十音順について
-
C言語の課題で、1年の秒数を計...
-
Notepad++の関数リスト表示でC...
-
VBAで仕様書は書きますか?
-
Excelに埋め込んだVBAのプログ...
-
UWSCで指定のフォルダを開きたい。
-
グループを均等に分けるには?...
-
シャープのポケコンについて
-
めんどくさがり屋はプログラマ...
-
読み込み中にアクセス違反が発...
-
排他的論理和 BCC(水平パリテ...
-
画像処理
-
ドロップダウンリストの文字を...
-
EXCELなどで「返す」という表現
-
科学技術計算の仕事について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
フローチャートで 変数に代入す...
-
フローチャートの菱形が狭い。。。
-
フローチャートで。
-
fortran errorについて
-
フローチャート以外の設計方法
-
二分法(FORTRAN)
-
fortran go to 文
-
フローチャート(本当に困って...
-
TeXでフローチャート
-
フローチャートが書けません
-
【fortran】フーリエ級数について
-
フローチャートのループ
-
フローチャートを書きたい
-
配列 x に入っているデータの最...
-
連立1次方程式を計算するFortra...
-
フロートチャートをVBAで表すには
-
カシオ fx-4800pにて
-
for文のフローチャート
-
FORTRAN subroutineと配列と繰...
-
正しい五十音順について
おすすめ情報