[過去ログ]
高校数学の質問スレ Part437 (1002レス)
高校数学の質問スレ Part437 http://rio2016.5ch.net/test/read.cgi/math/1721071007/
上
下
前
次
1-
新
通常表示
512バイト分割
レス栞
このスレッドは過去ログ倉庫に格納されています。
次スレ検索
歴削→次スレ
栞削→次スレ
過去ログメニュー
653: 132人目の素数さん [sage] 2024/08/04(日) 09:19:05.64 ID:wQ387Xvj デカルトの葉線の重心を計算 https://i.imgur.com/51eMWPw.png Wolframのソルバーで x^3+y^3-3*x*y=0の実数解を出して ループ形成する範囲での差を重みとして原点周りのモーメントを積分すれば 重心の座標が出せるはず、 というアルゴリズムでコーディング x2dy[x_] :=( y0=y /. Solve[x^3+y^3==3*x*y, y,Reals]; re=Simplify@Differences@Select[y0,0<= # <=2^(2/3) &]; If[x==0 || re=={},0,re[[1]]] ) NIntegrate[x*x2dy[x],{x,0,2^(2/3)}]/NIntegrate[x2dy[x],{x,0,2^(2/3)}] エラーが返ってきてバグの原因がわからんので 同じアルゴリズムでRで書くことにした。 Rにはソルバーはないが二次項がないのでカルダーノの式を使う。 計算結果 > S=integrate(x2dy,0,x1,rel.tol = 1e-12)$value > G=integrate(\(x) x*x2dy(x)/S,0,x1,rel.tol = 1e-12)$value > c(G,G) [1] 0.8061331 0.8061331 Rでの算出なので小数表示になる。 東大卒やエリート高校生による検証を希望します。 http://rio2016.5ch.net/test/read.cgi/math/1721071007/653
メモ帳
(0/65535文字)
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 349 レスあります
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.009s