1.惑星は太陽を一つの焦点とする楕円軌道を描いて太陽のまわりをまわるというプログラムは
#include <stdio.h>
#include <math.h>
#define GM 4*M_PI*M_PI
#define M 1.0
#define dT 1.0/256
#define QX0 1.0
#define QY0 0.0
#define PX0 0.0
#define PY0 sqrt(GM)*M
int main(void)
{
FILE*fout;
double qx, qy;
double px, py;
double T;
double qxk1, qyk1, pxk1, pyk1;
double qxk2, qyk2, pxk2, pyk2;
double qxk3, qyk3, pxk3, pyk3;
double qxk4, qyk4, pxk4, pyk4;
double tqx, tqy, tpx, tpy;
double q, qi3;
fout = fopen("data16.dat","w");
qx=QX0;qy=QY0;px=PX0;py=PY0;
for(T=0.0;T<10.00;T+=dT)
{
tqx = qx; tqy = qy; tpx = px; tpy = py;
q = hypot( tqx, tqy ); qi3 = 1.0/(q*q*q);
qxk1 = tpx / M * dT;
qyk1 = tpy / M * dT;
pxk1 = -GM * M * tqx * qi3 * dT;
pyk1 = -GM * M * tqy * qi3 * dT;
tqx = qx+0.5*qxk1;
tqy = qy+0.5*qyk1;
tpx = px+0.5*pxk1;
tpy = py+0.5*pyk1;
q = hypot( tqx, tqy ); qi3 = 1.0/(q*q*q);
qxk2 = tpx / M * dT;
qyk2 = tpy / M * dT;
pxk2 = -GM * M * tqx * qi3 * dT;
pyk2 = -GM * M * tqy * qi3 * dT;
tqx = qx+0.5*qxk2;
tqy = qy+0.5*qyk2;
tpx = px+0.5*pxk2;
tpy = py+0.5*pyk2;
q = hypot( tqx, tqy );
qi3 = 1.0/(q*q*q);
qxk3 = tpx / M * dT;
qyk3 = tpy / M * dT;
pxk3 = -GM * M * tqx * qi3 * dT;
pyk3 = -GM * M * tqy * qi3 * dT;
tqx = qx+qxk3;
tqy = qy+qyk3;
tpx = px+pxk3;
tpy = py+pyk3;
q = hypot( tqx, tqy );
qi3 = 1.0/(q*q*q);
qxk4 = tpx / M * dT;
qyk4 = tpy / M * dT;
pxk4 = -GM * M * tqx * qi3 * dT;
pyk4 = -GM * M * tqy * qi3 * dT;
qx += (qxk1 + 2*qxk2 + 2*qxk3 + qxk4)*(1.0/6);
qy += (qyk1 + 2*qyk2 + 2*qyk3 + qyk4)*(1.0/6);
px += (pxk1 + 2*pxk2 + 2*pxk3 + pxk4)*(1.0/6);
py += (pyk1 + 2*pyk2 + 2*pyk3 + pyk4)*(1.0/6);
printf("qx=%f qy=%f px=%f py=%f \n",qx,qy,px,py);
fprintf(fout,"%f, %f\n",qx,qy);
}
fclose(fout);
return 0;
}
と作れたのですが、
2.面積速度は一定である。
3.惑星の公転周期Tの2乗と軌道長半径aの3乗の比T^2/a^3は太陽系のすべての惑星に共通で一定のあたいになる。
を数値的に示せ。
というプログラムが作れません。
教えてください。
No.1ベストアンサー
- 回答日時:
>2.面積速度は一定である。
3.惑星の公転周期Tの2乗と軌道長半径aの3乗の比T^2/a^3は太陽系のすべての惑星に共通で一定のあたいになる。
を数値的に示せ。
自力でどこまで計算出来ていますか?
http://cms.phys.s.u-tokyo.ac.jp/~naoki/CIPINTRO/ …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 経済学 異質財のクールノー競争について教えてください。 3 2022/11/24 10:35
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
- 経済学 経済学の問題について教えてください。 4 2022/11/23 19:16
- 物理学 運動量演算子の交換関係について 5 2022/06/14 17:17
- JavaScript ソースコードのいじる場所が分かりません。 1 2022/12/23 02:06
- 経済学 ミクロ経済の問題 準線形効用関数 4 2022/08/14 00:58
- 防災 【自衛隊のPX品のPXってどういう意味ですか?】 3 2022/07/02 13:12
- 楽器・演奏 電子キーボードのCASIO Privia PX-S1100 とPrivia PX-S1000のどちら 1 2023/07/31 20:56
- プリンタ・スキャナー エプソンの PX 105 プリンター ブラックのみカートリッジ使用する設定 教えてください 1 2023/03/20 17:02
- その他(教育・科学・学問) 画用紙 26cm×37.5cmは 何px でしょうか? 何×何ですか? 1 2023/01/09 09:53
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
物理について 物理学 宇宙 地球...
-
物理学 地学 地球 物理 宇宙 地...
-
物理学 宇宙には多くの惑星系が...
-
ぶんぶんごまがブンブンまわる条件
-
潮の満ち干きと人の命の関係に...
-
なぜ地球型惑星と木星型惑星に...
-
【エクセル】CTRL+FでA...
-
どうして木星や土星の自転速度...
-
天動説の矛盾
-
「月」を使ったきれいな言葉、...
-
ある日、突然、公転軌道から外...
-
中三理科。 月の動きは西から東...
-
天の川銀河の端っこの天体
-
次の満月はいつですか??教え...
-
20億光年
-
世界一月が大きく見える場所の画像
-
地球は丸い~遠くの船の見え方
-
太陽の自転の向き 地球の自転の...
-
地球が丸いことの証明について ...
-
昼間の月
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ぶんぶんごまがブンブンまわる条件
-
放物線と双曲線の力学上の違い...
-
新物理入門持っている人教えて...
-
月の公転の速さ
-
ダークマターの回収
-
惑星の運行は永久機関ではない?
-
質問です。実測値がケプラーの...
-
ガリレオが楕円運動を否定して...
-
太陽は惑星を伴い銀河を63度の...
-
赤道付近の時計は、回転運動分...
-
物理学基礎の単位換算のコツを...
-
物理について 物理学 宇宙 地球...
-
物理について 地球 宇宙 物理学...
-
万有引力の法則の発見
-
C言語でケプラーの公式
-
ケプラーの第三法則
-
面積速度一定の法則の証明について
-
楕円の二つ焦点の両方に太陽が...
-
潮の満ち干きと人の命の関係に...
-
1ヶ月って何週間ですか?
おすすめ情報