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で質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
-
許せない心理テスト
私は「あなたの目の前にケーキがあります。ろうそくは何本刺さっていますか」と言われ「12本」と答えたら…
-
フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
あなたが普段思っている「これまだ誰も言ってなかったけど共感されるだろうな」というあるあるを教えてください
-
映画のエンドロール観る派?観ない派?
映画が終わった後、すぐに席を立って帰る方もちらほら見かけます。皆さんはエンドロールの最後まで観ていきますか?
-
海外旅行から帰ってきたら、まず何を食べる?
帰国して1番食べたくなるもの、食べたくなるだろうなと思うもの、皆さんはありますか?
-
天使と悪魔選手権
悪魔がこんなささやきをしていたら、天使のあなたはなんと言って止めますか?
-
C言語でケプラーの法則
物理学
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
ダークマター論進展
-
ぶんぶんごまがブンブンまわる条件
-
放物線と双曲線の力学上の違い...
-
C言語でケプラーの公式
-
ダークマターが太陽系にないのって
-
太陽系がほぼ同一平面状にある理由
-
惑星の運行は永久機関ではない?
-
太陽は惑星を伴い銀河を63度の...
-
万有引力の法則の発見
-
近日点で惑星の速さが最大にな...
-
運動量保存則がちょっとわから...
-
新物理入門持っている人教えて...
-
2球体間の静電容量を求める
-
楕円の二つ焦点の両方に太陽が...
-
ケプラー時間について
-
ガリレオが楕円運動を否定して...
-
C言語でケプラーの法則
-
衛星がSwing byで運動エネルギ...
-
面積速度一定の法則の証明について
-
地球の公転を証明する観測について
おすすめ情報