![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?c9bd177)
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は太陽系のすべての惑星に共通で一定のあたいになる。
を数値的に示せ。
というプログラムが作れません。
教えてください。
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_03.png?c9bd177)
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を見た人はこんなQ&Aも見ています
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
惑星の運行は永久機関ではない?
-
面積速度一定の法則の証明について
-
物理学基礎の単位換算のコツを...
-
潮の満ち干きと人の命の関係に...
-
自転しない星ってありますか? ...
-
1ヶ月って何週間ですか?
-
月の下部分が欠けて見えるんで...
-
次の満月はいつですか??教え...
-
世界一月が大きく見える場所の画像
-
水金地火木~ 「冥海」から「...
-
天動説の矛盾
-
地球の自転方向にそって飛行機...
-
太陽から半径100光年の範囲に恒...
-
月の名前
-
中国のルナカレンダー(旧暦?...
-
彗星とダークマター
-
宇宙人は地球へ来れるのでしょ...
-
2036年小惑星・アポフィス...
-
太陽に生物は存在するか?
-
ウミガメや珊瑚などの産卵が満...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ぶんぶんごまがブンブンまわる条件
-
放物線と双曲線の力学上の違い...
-
第三宇宙速度はどうやって求め...
-
C言語でケプラーの公式
-
ダークマターの回収
-
物理学基礎の単位換算のコツを...
-
太陽は惑星を伴い銀河を63度の...
-
C言語でケプラーの法則
-
ケプラー時間について
-
太陽系がほぼ同一平面状にある理由
-
月の公転の速さ
-
近日点で惑星の速さが最大にな...
-
新物理入門持っている人教えて...
-
物理について 物理学 宇宙 地球...
-
物理の惑星の問題なのですが、...
-
ケプラーの第三法則
-
太陽の末期には巨大化し地球を...
-
面積速度一定の法則の証明について
-
ケプラーの法則(3)の証明方法
-
運動量保存則がちょっとわから...
おすすめ情報