趣味で中性子星の限界質量を作るプログラム作りたいから手伝って

ああ、今理解した、マクロか

以後f(,,)はvにg()は-ｘとして置き換えるってだけのマクロ文だわ

マクロス？

ただの置き換えだよ、プログラムの序盤で置換するだけの命令

あとルンゲクッタの方がいいんだけど
はじめは、オイラー法とかで試してみては？

>>137　じゃあただの実数ですか？

ただの変数ですvと-xが実数なら実数です

>>138　うん

>>140　なんて宣言すればいいの？

置き換えたい文字列 = 置き換え後

例)
func(x) = x+20　と記述すると
以後
func(x)*2　は
（x+20）*2　と解釈される

>>114のエラーの原因
line28は配列の宣言構文が間違ってる
ｘという名前の配列を作りたいならば
real x(8)　と記述すべき

他の３つは多分rhoの宣言を忘れてる

よって寝る

みんなありがとう！！コンパイルできたよ！！
でも新たな問題が・・・
計算結果がNaNになってしまうのです！！
うｐするから誰かどこが間違ってるか教えて！！

http://loda.jp/tomfile/?id=86

メインの計算部分は
c loop----------
って書いてあるところの下。
初期条件はちゃんと入ってるみたいです

誰か勝手にデバックしてくれたら神と呼びます

あああ　なんでうまくいかん

この先なんねんもパソコンとにらめっこしながらデバックし続ける人生だとおもうと
悲しくなる

でもパフュームのみんなも売れるまで何年も頑張ったんや　おいらも頑張るじょ！！！！！！！！！！！！！！！！！！！！

よくわからないけど、このプログラムで計算するとだいたいどれくらいになるかはもうわかってるの？

http://www.phys.u-ryukyu.ac.jp/~apweb_user/diss/s2008.html
このサイトの解系列ってグラフになるはず

最大質量は0.7solorくらいになるはずです

そうなんだ
まあ頑張って

っておーーーーーーーーーーーーーーーーーい！！！！！！！！！！１それだけかーーーーーーーーーーーーーーい！！！！！！！！！！！！！！！！！！！！！！(爆死

NaNってことはどっかでゼロ除算してるとかだろ

なるへそ！！！！！！！！！！！！！

>>156　デバック頼むわ＾ｐ＾

何がなるへそだよｗ

ざっと目を通したけどx()y()z()のマクロ宣言が見当たらない

一番下に書いてあるよ

と思ったら関数にしてたのか。。。。

じゃぁｘ（）ｙ（）ｚ（）関数内にバグがないかそれぞれの関数を一度だけ呼び出すプログラム作って確認
場所特定はデバッグの基本

大学のサーバ使ってるからデバックしてるときは２ちゃんに書き込めないんだよ
代わりにやってくれよ

いやだよ実行環境ないし、フォートランのコード今日始めてみたし

うんちょこちょこちょこぴーーーーーーーーーーーーーーーーー！！！！！！！！！！！！！！！！！！！！！！１

ああ、分かった、関数ｘで最初の方に
real(8) P,M,r,rho　って書いてるだろ、この３つは引数で渡されてるから宣言の必要はないと思う
他の関数も同様

ｃｖｃ２も初期化できてない、引数で渡してくるかなんかしないと

バグっていうか関数の使い方が出来てないわ
関数内で宣言してる変数の殆どが意味成してないけど

もうツんだ

なんだ諦めたのか、アルゴリズム書いといたら情報系の誰かがコード書くかもよ

ほんと、もう諦めたの？

結局諦めたの？

お前らｋｓだな

今週は忙しいので金曜日くらいにデバック再開する。

今週忙しすぎてやばす

プログラムなら手伝うから
アルゴリズムを簡単でいいから教えて

アルゴリズムがよくわからないから困ってるわけだけど。

単位はcgs単位系

初期条件の読み込み
↓
変数、定数のノーマライズ
↓
方程式の読み込み
dP / dr = - (g / r**2) * (rho + (P/cvc2))
& * (M + (pi4 * P / cvc2))
& / (1.0d0 - (2.0d0 * goc2 * M / r))
c
dM / dr = pi4 * r**2 * rho
c
drho / dr = - rho * (g / r**2) * (rho + (P/cvc2))
& * (M + (pi4 * P / cvc2))
& / (1.0d0 - (2.0d0 * goc2 * M / r))
& / (gam * P)
↓
ルンゲクッタ法ループ開始
↓
出力

cvcは光速の二乗、goc2はG/c^2

初期条件をどうするか問題
r=0にしちゃうと式が発散してしまう

わかった
時間あるときにやってみるよ

あ

Simulimkでやれよ