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

頼むよ

ワトソン君

スレタイ
大人の人に乗ってもらう悩み相談スレ
内容
お願いします

誰か立ててください

方程式系は以下の通り
0=-(Γ/h) * 4πr^2 * ∂p/∂m - M/r^2 - 4πp
∂M/∂m = 4πr^2 *　∂r/∂m ρ(1 + e)
P=P(ρ)

変数はMとρ？
他の文字の扱いをどうすればいいのか分からないので
実際に解いてる人の意見を聞きたい

原子核系の研究室の人なら解いたことあるはず

おい、こういう過疎のときだからこそ、何かを作り上げようぜ

現状の問題は以下の三つ

・独立変数以外の文字の扱い
・初期条件をどうするか
・状態方程式をどうするか

一週間くらいで作れたらいいなと思ってる

そんなんよりエロゲーでも作れよ

絵かけないもん

中性子星がキーワードのSFエロゲ

限界を超えるとブラックオｈールになるの？

>>11　そうだよ　でもまだ限界質量はわかってないんだ

エクセルでできねえの

fortranで作りたいと思ってる　まあエクセルでも出来るかもしれんけど

ちゅうせいしせいマンコ
あまりにきょうりょくなじゅうりょくなためチンコにもうれつなちょうせきりょくがはたらき
チンコはスパゲッティじょうにのばされてしまう
ないぶはとてつもないあつりょくがあるのでいっしゅんにしてイッてしまう
そうにゅうしたらさいごピストンうんどうしつづけるしかない

おっ　人集まってきたな

本当は状態方程式を用意しないといけないんだけど
データを持ってないので
P=ρ^γとでも置いて解こうかと思ってる

おい、今こそ既卒の力をみせるときだぞ

誰かー；ω；

そんなものよりエロゲ作れよ

じゃあ、真冬に捨てられてた猫を拾ってあげたら
次の日、謎の美少女が訪ねてきて洗濯とか掃除してくれるエロゲ作るか？

（∂Д∂）

よっしゃ作ってくるから4時間ほど起きて待っててくれ

ヒロインは黒髪のロリ系で頼むな！！

数学科のやつ恋や！！

マジで作ろうよー

じゃあ>>15のエロゲ作るわ

せめて大学生のいる板でやりなさい・・・

お前らｋｓだな

方程式だけみると簡単なんだけどなぁ

あああ

ニュートンが万有引力を発表した時に
質量が十分大きければ光でも抜け出せない、と仮設がたてられた

日本では江戸時代かな？
そんな時にブラックホールを数学的に予測していたなんて凄い

>>32　おっぱいみせて

あれ、引力と重力って同じものって認識でいいんだっけ

おなじ

私の胸は可愛くないし

可愛いかどうかはオレが決めるんだよ！！！！！！！！！！！！！！！！！！！！

何その他は可愛いみたいなニュアンス

>>1
まずは中性子星の限界質量を計算してる、あるいは記述してる論文を探しては如何？

>>39
理論は分かってるんだよ
初期条件とか状態方程式の細かいところをどうするか悩んでるだけ

>>40
論文は見つけてるわけだ、ならその著者にメールで質問してみるのは？

趣味でやってるのにそこまでしたくねえよｗｗだいたいどうやって解くのか自分で考えるから面白いんだろ

方程式はこれなんだってば
0=-(Γ/h) * 4πr^2 * ∂p/∂m - M/r^2 - 4πp
∂M/∂m = 4πr^2 *　∂r/∂m ρ(1 + e)
P=P(ρ)

連立微分方程式

>>42
趣味だって良いと思うけどな。
研究者って偏屈だと思われがちだけど、自分の専門分野の事を尋ねられると大抵の人は喜ぶよ。

大生のみんなで作りたいねん！！

信用してないくせに。。。

F欄文系の俺にできることある？

0=-(Γ/h) * ♀д♀ * ∂p/∂m - M/r^2 - 4πp
∂M/∂m = ♀д♀ *　∂r/∂m ρ(1 + e)
P=P(ρ)
これにξ＝*^-^)2K^2を代入すれば重力の限界質量とシュバルツシルト半径を求められるよ

一から十まで完璧にアルゴリズムを説明してくれたら作ってもいいよ

まず何の方法使って方程式解くかが謎やねん

dP(r)/dr = -G/r^2( ρ(r) + P(r)/c^2 ) (M(r) + 4πr^2P(r)/c^2 )( 1 - 2GM(r)/c^2r )
dM(r)/dr = 4πr^2ρ(r)
解きやすい形に変形してみた
独立変数はM(r),ρ(r)

これをどうやって解けばいいかが分からない

>>47
設計図どおりにコード作るくらいなら出来るんじゃないでしょうか

底辺かける高さ÷2

>>53　おっぱいみせて

そもそも初期条件なんかいらないのか
勘違いしてた

谷間？

今自分で考えてる案は
P(r)i - P(r)j / �决 = -G/r^2( ρ(r) + P(r)/c^2 ) (M(r) + 4πr^2P(r)/c^2 )( 1 - 2GM(r)/c^2r )
M(r)i - M(r)j / �决 = 4πr^2ρ(r)

と差分化してニュートンラプソン法で解く方法

�决は100cmにしときましょうか

>>56　谷間でもいいよ＾＾

僕もおっぱいみたい！

みんなで一緒にみよう＾＾

コテなのに頭いいの？

×コテなのに　○コテだから

おっぱい！！おっぱい！！

おっぱいまだー？？ﾁﾝﾁﾝ

今ちんぽ出してます！早くしないとこの子死んじゃうんです！！
助けてください！！！！＞＜＞＜＋＋＋っっっｌ；；；

もう目ニコには期待しない

ビップでdeSlove使えって言われたけど意味分かりません

フォートランでルンゲクッタのプログラムっての見つけたぞ

dP/dr = f(P,r) , P(r=r0) = 0

P(r + dr) = P(r) + hP( f1 + f2 + f3 + f4)

これがルンゲクッタっていうらしい

有限要素法なら

手抜き紙芝居でいいならフラッシュエロゲ作るぜ
絵は俺が書くから話とキャラ設定くれよ

登場人物

ドラドリ・・・一人暮らししている大学生
小夜子・・・近所に住んでる小学生の女の子。暗くて友達がいない。

話

ころんで怪我をしてる小夜子を手当てして仲良くなる。
友達がいない同士でいっしょに遊ぶ。
いちゃいちゃする。
小夜子の保護者が出てきて強引に別れる。
そして・・・

セックス！！！！！！！！！！！！！！！！！！！！

大生板にまともな理系がいるとでも思ってるのか>>1は・・・
文系のおれはググっても一ミリたりともわからんけどな

わからん方程式はマトラボにぶっこめばいいってじっちゃんがいってた

ルンゲクッタが分からない

専門用語も小難しい式もよくわからんから何したいのかよくわからんかったけど
多項式で近似してプログラムに落として計算しようって話？

そんなに難しい式でもないだろ
余計な項とか落とすと
dP(r)/dr = -Mρ/r^2
dM(r)/dr = 4πr^2 ρ
P(r) = ρ^n

これを解きたいの

数値解析学でルンゲクッタ法すら習わない大学生て…

そんな名前の講義ねーし

見慣れなさ過ぎて式が見えてこないｗ
Ｐ（）とかＭ（）って関数なの？右辺のＭは何でつか

じゃあ物理の式じゃなくて単純にxとかyで書く

dx(t)/dt = f(x,y,t)
dy(t)/dt = g(x,y,t)
x(t) = z(t)^n

こんな感じ

うんちゲーリー

文系の俺に何をどうしたいのか初めから教えてくれ
そうしたら、機械的にプログラムを書くよ

居室で近くの席の人が超新星爆発とかニュートリノとかブラックホールの
数値計算してるけど，自分は全然違う分野だからさっぱりわからんね

これを解いたところで何の得になるの

♀д♀ってやつスレ見てる感じ似非理系そう
本当は何もわかってなさそう

理系の先端科学と哲学の境目が分からん

ドラドリってホントに理系なのか気になる

>>87
それ言ってたら何の進歩もない

>>87
まさにこれ
とりあえずこれが解けてなにが嬉しいの？
それによって誰が得するの？

ドラドリの肩持つのも癪だが、それが趣味ってもんじゃないの

オレが理系じゃなかったら誰が理系なんだよ

新参増えたなぁ

>>86
超新星爆発すると、中心部に中性子星が出来るんだけど
爆発に失敗すると周りの物質が外に飛び散らないで、中性子星に降り積もっていくんだよ
中性子星には限界質量があるから
降り積もった物質が限界質量を超えるとブラックホールになる

ちなみに、そのときに陽子が電子と反応してニュートリノが生成される

>>88　とりあえずドラドリについて勉強しろ

ひも理論とか関係あるの？

関係ないよ

100

学問で言うとなんて分野の話？
入門書的な教科書とか参考書ある？

中性子星の話？
宇宙物理の本にふつうに載ってる。
方程式の名前はトルマン・オッペンハイマー・ボルコフ方程式
wikiにも載ってるがな

理学部宇宙物理学科みたいな感じ？
ここ卒業する奴らはどんな職につくんだい？

梅干しの種を除く仕事

ドラドリって超おもしれーな

だろ？？

将来ノーベル賞でも取るの？

まーねぃ

ドラドリさん宇宙レベルのギャグってありますか？

おちんちん歯磨き！！！！！！！！！！！！！！！！！！！！！！！これが今年のイチオシギャグやね

僕のおちんちんも磨いてください（ポロリータ

ドラドリのおてぃんてぃんうｐ

数学板で解き方ならってきた
dP(r)/dr = f(P,M,ρ,r)
dM(r)/dr = g(P,M,ρ,r)

P = Kρ^n
log P = logK + nlogρ
1/P * f - logK = n * 1/ρ dρ/dt

dρ/dt =ρ/n (1/P * f - logK)

こんなエラーが出てコンパイル通らない。誰か助けて＞＜

fortcom: Error: 230124.f, line 28: A specification statement cannot appear in the executable section.
real(8) x(0:pmax,0:mmax,0:rhomax,0:rmax),
------^
fortcom: Error: 230124.f, line 56: This name has not been declared as an array or a function. [X]
x(P,M,rho,r) = - (g / r**2) * (rho + (P/cvc2))
------^
fortcom: Error: 230124.f, line 59: This name has not been declared as an array or a function. [Y]
y(P,M,rho,r) = pi4 * rho * r**2
------^
fortcom: Error: 230124.f, line 60: This name has not been declared as an array or a function. [Z]
z(P,M,rho,r) = rho / gam * (x(P,M,rho,r)/P - dlog(k))
------^

A specification statement cannot appear in the executable section. ってなんだよ

mathematicaあるけど使う？

いやもうプログラム書いたからいい。それより上のエラーメッセージを教えてくれ

すげーアホな質問だけど
配列でx(実数、実数)みたいなのって許されんの？

実数を引数にする配列ってことね

無理だよ

じゃあコレおかしくない？？

c 2階常微分方程式のときは、以下のようにして、連立の
c 1階微分方程式にして解く。
c d^2x/dt^2 = - x --->> dx/dt = v
c dv/dt = -x
c 一般的に dx/dt = f(x,v,t)
c dv/dt = g(x,v,t) の形を想定。

parameter(pi=3.14159265358979) ! pi の値をここで宣言しておく。
f(x,v,t) = v
g(x,v,t) = -x

c 初期条件
t = 0.0
x = 0.0
v = 1.0

c t = 0 ~ 2 pi までを n 分割して計算。
n = 8
c きざみ幅
h = 2.0 * pi / n

できないんじゃない？
それより俺にはこのスレで何をやってるのかすら分からない

FORTRANか？
line28のエラーは宣言式だと思うけど、多分構文が間違ってる
他の３つはrhoの型の不一致ってとこだろ、ポインタ関連見直せ

f(x,v,t) = v
g(x,v,t) = -x
の左辺って配列じゃないの？なんで右辺は実数なの？

ルンゲクッタが分からん！！
だれかソースコードくれ！！

実数の配列って意味で言ってるなら出来るけど
double x[3]={1.0, 2.2, 3.3} みたいな具合に、
double x[1.5]　とかは無理
C系列の言語しか分からないから悪いけど

>>125
>>121はとあるサイトに載ってたルンゲクッタのサンプルコードなんだけど
おかしくない？
なんでfは配列なのに左辺は実数なの？

>>126
配列なの？

そもそも何語で書いてるか分からないから定義文なのか代入文なのか分からない
関数型言語っぽいけどコメントの構文にc使うとか普通の言語じゃねぇ。。。。

配列じゃないの？だってf(x,v,t)って書いてあるし

コメントにcはフォートランだよ

邪魔なコメント取るとこんな感じか
//==============================
parameter(pi=3.14159265358979)
f(x,v,t) = v
g(x,v,t) = -x

t = 0.0
x = 0.0
v = 1.0

n = 8
h = 2.0 * pi / n
//==============================
最初の２行はf()とg()の定義文っぽいな
それ以降は変数の定義だろう
肝心の計算部分のコードがないが

>>129
Cで言うと関数っぽいのを想定するんじゃないの　f　は？
コメントアウトに書いてある微分方程式のプログラム部は無いの？

このサイトでつ
http://windom.phys.hirosaki-u.ac.jp/~kasai/wiki.cgi/su04?page=%BE%EF%C8%F9%CA%AC%CA%FD%C4%F8%BC%B0%A4%CE%BF%F4%C3%CD%B2%F2%CB%A1

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

以後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でやれよ