![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
以下の問題について
『原点に質量Mの天体があり、点(x.y)上にある質量mの天体が万有引力で引かれているとする。M>>mとして、質量Mの天体は動かないと仮定する。そのとき、質量mの天体の運動方程式は以下のようになる。
mdVx/dt=-GMmx/r^3
mdVy/dt=-GMmy/r^3
ここでGは万有引力定数でr=√x^2+y^2である。
問: 簡単のためにGM=1として、初期条件x(t=0)=0:5、y(t=0)=0:0、vx(t=0)=0、vy(t=0)=1:63のもとで、上の微分方程式をオイラー法によって解け。適当な力学の教科書を見れば、この初期条件では、答えは閉じた楕円軌道になることが、この系の持つ力学的エネルギーの値から分かる。また、長軸、短軸、周期なども、力学的エネルギーや角運動量の値から計算できる。この答えを実現するためには時間刻み
Δtの値はどの程度でないといけないかを調べよ。』
使用するc++のソースコードに間違いがあるのですが、よくわかりません。知識のある方に教えていただきたいです。↓ソースコード
//n.C
#include <iostream>
#include <cmath>
using namespace std;
main()
{
int n;
float vx=0.0;
float vy=1.63;
float x=0.5;
float y=0.0;
float t=0.0;
float dt=0.05;
float r, r3;
float x1, y1;
cout << x << " "<< y <<"\n";
for (n=1; n<=100; ++n){
x=x+vx*dt;
y=y+vy*dt;
r=sqrt(x*x+y*y);
r3=r*r*r;
vx=vx-(x/r3)*dt;
vy=vy-(y/r3)*dt;
cout << x << " "<< y <<"\n";
}
}
よろしくお願い致します!!
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- その他(プログラミング・Web制作) 物理の斜方投射のシミュレーションにおける位置や速度の単位について 4 2023/05/31 09:50
- その他(プログラミング・Web制作) Pythonでのかんたんな物理シミュレーションについての書籍 5 2023/06/02 07:37
- その他(プログラミング・Web制作) 物理の斜方投射で目盛りに数値を入れたい 2 2023/05/27 06:32
- その他(プログラミング・Web制作) 物理の斜方投射の目盛り線とx軸、y軸の追加について 3 2023/05/26 21:11
- その他(プログラミング・Web制作) Pythonにおける物理のシミュレーションでの単位変換について 2 2023/06/02 17:11
- 物理学 なめらかな水平面の床の上に、質量 200 g の物体がある。床の面を xy 面とし、鉛直方向に z 1 2022/07/23 11:28
- その他(プログラミング・Web制作) ボールの動きがスムーズに動いてかつ目盛り線描画を維持するためには 4 2023/05/31 10:01
- 物理学 物体に一定の大きさfの力をx軸の正の向きに加える。またこの物体には抵抗係数がγの速度に比例する抵抗力 2 2023/07/06 04:01
- その他(プログラミング・Web制作) Pythonによる物理の斜方投射の位置座標表示について 2 2023/06/05 12:46
- 物理学 鉛直下向きにx軸をとる。真上に投げた物体の運動方程式をかけ。ただし物体には抵抗係数αの空気抵抗が働く 5 2023/06/20 02:04
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
C言語 exitの使い方
-
数字の位ごとの値を表示するプ...
-
Excel-vba 文字列と変数を...
-
世界のナベアツ
-
VBAの定数の使い方で、計算値を...
-
Access2003 オートナンバーの現...
-
VB6.0-整数と余りを求める
-
フォームを開くときに、コンボ...
-
VBAで配列のNULL判定
-
Exit Subのような・・・。
-
Nullってどういう意味ですか?
-
c言語 偶数個
-
ファイルの特定行を取り出したい.
-
jsp~jspにhiddenを使って変数...
-
エラーの意味は? Lvalue req...
-
VBAでダブルコーテーション入り...
-
VB6で配列を文字列に変換する方...
-
[VBS]変数を定数に変換する方法...
-
vbsでテキストファイル内の文字...
-
エクセルへのデータの貼付時に...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
C言語 exitの使い方
-
フォームを開くときに、コンボ...
-
Excel-vba 文字列と変数を...
-
エクセルのマクロについて教え...
-
VB6.0-整数と余りを求める
-
足して100になるような乱数のア...
-
数字の位ごとの値を表示するプ...
-
VBAで配列のNULL判定
-
世界のナベアツ
-
ラジオボタンの値の取得につい...
-
相関係数p値の出し方
-
1つ前の値を変数に保存する方法
-
VBAの定数の使い方で、計算値を...
-
DWORDって
-
4択問題のプログラムでランダム...
-
Access2003 オートナンバーの現...
-
C#で動的にコントロールを取得...
-
スピンボタンで小数点
-
UWSCのcallについて
-
フリーランタイマーの時間差分...
おすすめ情報
回答ありがとうございます。これは課題として出されたものでして、前提として「このソースコードには間違いがあります」とされているのですが質問文に記載するのを忘れていたようです。申し訳ありません。