fortranでペットボトルロケットの軌道計算をしています。
コンパイルは問題ないのですが、正しい数値が出てきません。
なぜなのか教えてください。
e=1000
S=0.0006
m=0.055
P0=100000
L=0.002
Ps=600000
Vs=0.0008
P=600000.0
V=0.0008
a=0.0
d=0.0
h=0.0
dt=0.01
t=0.0
write(6,600)P,V,a,d,h,t
600 format(4f10.4)
do i=1,10000
if(V>0)then
w=(2*(P-P0)/e)**(1/2)
else
w=0
end if
Mm=m+e*v
V=V-w*S*dt
a=a+e*S*w*w/Mm
d=d+a*dt
h=h+d*dt
P=Ps*(L-Vs)/(L-V)
t=t+dt
write(6,600)P,V,a,d,h,t
enddo
end
No.2ベストアンサー
- 回答日時:
各文字が何を表しているのか
どういう微分方程式を解こうとしているのか
「正しい数値が出ない」とは具体的にどういう結果が出力されたという意味か
みたいな事を貴方の書いたプログラムだけから推測するクイズを解く気にはなれないので、クイズを解かずにわかる範囲の事しか書けませんが、
言語によっては整数/整数の結果が常に整数(商を返す)になります。つまり1/2が0.5ではなく0になります。fortranの仕様は覚えてませんが、この事を意識せずに1/2と書いたのならこれがいくらになっているのか確かめた方がいいでしょう。
>a=a+e*S*w*w/Mm
>d=d+a*dt
>h=h+d*dt
この部分の計算の順番はこれで正しいですか?
例えば、hの式にあるdはその直面で計算したdの値(多分、次の時刻での値)が使われます。
No.3
- 回答日時:
Fortran では 1/2 は 0 になります>#2. あと, 暗黙の型宣言が有効になっていると L は整数型.
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Excel(エクセル) 日付で矢印マクロ 4 2023/07/25 16:47
- その他(プログラミング・Web制作) 物理の斜方投射で目盛りに数値を入れたい 2 2023/05/27 06:32
- その他(プログラミング・Web制作) Pythonでのかんたんな物理シミュレーションについての書籍 5 2023/06/02 07:37
- 物理学 相対論的運動方程式 1 2022/07/04 06:20
- その他(プログラミング・Web制作) 物理の斜方投射のシミュレーションにおける位置や速度の単位について 4 2023/05/31 09:50
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
- PHP PHP MySql 画像を取得 1 2022/06/04 14:05
- C言語・C++・C# 大量のデータを読み込んで表示する速度を改善したい 8 2023/05/07 13:29
- 数学 x=r・cosθの2回微分 θ=ωtとすると? 5 2022/05/10 23:53
- C言語・C++・C# クラスのメンバ変数を基準に並べ替えをしたい 5 2022/12/25 17:40
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
電流の時間微分、電圧の時間微分
-
質量流量の記号「・ の読み方を...
-
d^2r/dt^2の意味
-
最後のdv/dtは何でしょうか。
-
微分積分のdの意味
-
誘導起電力について 誘導起電力...
-
単振り子の運動方程式
-
機械力学の問題です!!!
-
微分記号“d”について
-
d/dx=dt/dx * d/dt =d/dt * dt/...
-
高校物理、位置、速度、加速度...
-
ブックオフは潰れないよね? キ...
-
物理:ひもでつながれた2物体の...
-
交流電圧の時間変化が正弦波の...
-
EXCEL上の数字を自動で振り分け...
-
運動方程式の成分表示
-
この計算がよくわかりません dx...
-
dx/dt=√(1-x^2)の一般解の求め...
-
シュレーディンガー方程式の変...
-
力学的エネルギー保存則の導出...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
電流の時間微分、電圧の時間微分
-
d^2r/dt^2の意味
-
ラプラス方程式
-
Debug.Printで表示される内容を...
-
極座標の運動方程式の計算の間...
-
質量流量の記号「・ の読み方を...
-
微分積分のdの意味
-
最後のdv/dtは何でしょうか。
-
加速度 a=dv/dt = (d^2 x) /dt^2
-
雨滴の運動質量が変化する落体...
-
力学について質問です。 1.棒の...
-
運動方程式を求めてください
-
機械力学の問題です!!!
-
微分記号“d”について
-
v^2-v0^2=2ax 今日この式を習っ...
-
運動方程式の微分積分の計算
-
運動方程式の成分表示
-
ブックオフは潰れないよね? キ...
-
EXCEL上の数字を自動で振り分け...
-
難しくてわからないので教えて...
おすすめ情報
w={2×(P-P0)/e}^1/2
だと思って入れました。
fortranでペットボトルロケットの軌道計算しなさい。というのと、
1/2×e×w^2=P-P0
だけ与えられて自分で組んだのでこれだけです。
何か足りないですか?