以下のようなプログラム(速度ベレル法による時間発展)を書いたのですが,エラーがでてしまいます.
どこが悪いのかを教えていただけたらと思います.
PROGRAM verlet
implicit none
integer::i,n,p
real(8)::dt,k,m,ratio,ti
real(8),allocatable::x(:),v(:),f(:)
allocate(x(0:n))
allocate(v(0:n))
allocate(f(0:n))
dt=0.01d0
n=1000
k=1.0d0
m=1.0d0
ratio=k/m
p=4
ti=0.0d0
x(0)=0.0d0
v(0)=1.0d0
f(0)=-k*(x(0))**(p-1)
do i=0,45
x(i+1)=x(i)+dt*v(i)+f(i)/(2.0d0*m)*dt**2.0d0
f(i+1)=-k*(x(i+1))**(p-1)
v(i+1)=v(i)+dt*(f(i)+f(i+1))/(2.0d0*m)
end do
do i=0,n
ti=i*dt
write(6,*) ti,x(i),v(i)
end do
deallocate(x)
deallocate(v)
deallocate(f)
END PROGRAM verlet
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- IT・エンジニアリング 三菱製PLC Qシリーズで技術的なことをご教示いただければ幸いです。 下記の回路で0.1秒の値をD0 1 2023/02/27 18:44
- 哲学 《自性清浄心》と《聖なる甘え》――ブッダとは目が覚めたという意味である―― 0 2023/05/23 04:15
- 物理学 電界による電子線の偏向についてなのですが 極板AとBの間隔をd0、極板Aの電位を-V0、極板Bの電位 1 2022/10/10 14:51
- その他(プログラミング・Web制作) 物理の斜方投射のシミュレーションにおける位置や速度の単位について 4 2023/05/31 09:50
- その他(プログラミング・Web制作) 物理の斜方投射の目盛り線とx軸、y軸の追加について 3 2023/05/26 21:11
- 数学 x=r・cosθの2回微分 θ=ωtとすると? 5 2022/05/10 23:53
- その他(プログラミング・Web制作) 物理の斜方投射で目盛りに数値を入れたい 2 2023/05/27 06:32
- 物理学 誘導起電力について 誘導起電力Vはファラデーの法則より、φを回路を貫く磁束として、 V=-(dφ)/ 1 2023/03/01 05:13
- C言語・C++・C# 大量のデータを読み込んで表示する速度を改善したい 8 2023/05/07 13:29
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
関連するカテゴリからQ&Aを探す
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
UserForm1.Showでエラーになり...
-
VBAでfunctionを利用しようとし...
-
【VBA】ワークブックを開く時に...
-
String""から型'Double'への変...
-
実行時エラー'-2147467259(8000...
-
文字列内で括弧を使うには
-
マクロで"#N/A"のエラー行を削...
-
Filter関数を用いた結果、何も...
-
実行時エラー 438 の解決策をお...
-
8876017c とエラー表示がされる
-
Excel実行時エラー[80004005]に...
-
ExecuteNonQueryメソッドの戻り値
-
VB6とVB.NETでNullの扱いが違う?
-
ACCESSで値を代入できないとは?
-
On ErrorでエラーNoが0
-
Excel2013で実行時エラー9...
-
オブジェクト型の変数にフォー...
-
Oracle Case文でのエラー(デー...
-
データベース操作時エラーについて
-
実行したSQLの結果を返す方法
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
UserForm1.Showでエラーになり...
-
String""から型'Double'への変...
-
VBAでfunctionを利用しようとし...
-
お助けください!VBAのファイル...
-
文字列内で括弧を使うには
-
実行時エラー 438 の解決策をお...
-
マクロで"#N/A"のエラー行を削...
-
On ErrorでエラーNoが0
-
エクセルVBA 「On Error GoTo...
-
VBA データ(特定値)のある最...
-
実行時エラー'-2147467259(8000...
-
【VBA】ワークブックを開く時に...
-
VBでSQL文のUPDATE構文を使った...
-
ACCESSで値を代入できないとは?
-
マクロの「SaveAs」でエラーが...
-
Excel vbaについての質問
-
インポート時のエラー「データ...
-
Filter関数を用いた結果、何も...
-
「実行時エラー '3167' レコー...
-
バッチファイルで、あるスクリ...
おすすめ情報