高校数学の質問スレ（医者・東大卒専用） Part438 (979ﾚｽ)
高校数学の質問スレ（医者・東大卒専用） Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
抽出解除
必死ﾁｪｯｶｰ(本家)
(べ)
自ID
ﾚｽ栞
あぼーん
445: １３２人目の素数さん [sage] 2024/12/22(日) 10:04:08.66 ID:rWRciXK9 >>439 では毒薬が２本のときのWolframのコードを書いてみ！ http://rio2016.5ch.net/test/read.cgi/math/1723152147/445
447: １３２人目の素数さん [sage] 2024/12/22(日) 11:05:38.63 ID:rWRciXK9 >>445 Wolframを使えないのが確定。 やはり、Fラン卒なのか？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/447
452: １３２人目の素数さん [sage] 2024/12/22(日) 12:33:00.32 ID:rWRciXK9 >>444 >>444 Wolframで思い通りに作図するスキルはないので R言語で最小値と最大値を求めてそのときの図を作図。 r=0.3のとき > c(min,max) [1] 0.3314719 2.0285281 最小のときの図 https://i.imgur.com/U9IMSgw.png 最大のときの図 https://i.imgur.com/olgz6rn.png 最小最大のときの形状が判明したので あとは一般解を求めればいい。 http://rio2016.5ch.net/test/read.cgi/math/1723152147/452
453: １３２人目の素数さん [sage] 2024/12/22(日) 12:35:08.21 ID:rWRciXK9 >>452 おまけ、R言語でのクズコード 理工系卒ならWolframで書けるはず。 但し、Fラン卒は除く。 r=0.3 source("toolmini.R") p=0.5+0.5i q=-0.5+0.5i r0=-0.5-0.5i s=0.5-0.5i fP=\(t) r*exp(-1i*t)+p fQ=\(t) r*exp(-1i*t)+q fR=\(t) r*exp(-1i*t)+r0 fS=\(t) r*exp(-1i*t)+s f=\(tt,verbose=FALSE,...){ P=fP(tt[1]) Q=fQ(tt[2]) R=fR(tt[3]) S=fS(tt[4]) if(verbose){ Plot(-1,1,zero=FALSE) Polygon(p,q,r0,s,Col=8) Cir(p,r,col=8) Cir(q,r,col=8) Cir(r0,r,col=8) Cir(s,r,col=8) Polygon(P,Q,R,S,...) pta(P);pta(Q);pta(R);pta(S) } ABC2S(P,Q,R)+ABC2S(P,S,R) } opt=optim(runif(4,-pi,pi),f) opt=optim(opt$par,f) tt=opt$par min=f(tt,TRUE,Col=2) opt=optim(runif(4,-pi,pi),f,control=list(fnscale=-1)) opt=optim(opt$par,f,control=list(fnscale=-1)) tt=opt$par max=f(tt,TRUE,Col=4) c(min,max) http://rio2016.5ch.net/test/read.cgi/math/1723152147/453
454: １３２人目の素数さん [sage] 2024/12/22(日) 12:40:56.83 ID:rWRciXK9 >>453 算出のために作った関数の動作確認の図 https://i.imgur.com/0mvHN8w.png http://rio2016.5ch.net/test/read.cgi/math/1723152147/454
455: １３２人目の素数さん [sage] 2024/12/22(日) 12:42:11.71 ID:rWRciXK9 毒薬が２本のときのWolframのコード まだぁ？ http://rio2016.5ch.net/test/read.cgi/math/1723152147/455
463: １３２人目の素数さん [sage] 2024/12/22(日) 18:41:54.83 ID:rWRciXK9 >444に答えてあげればいいのに。 それもできないのでFラン確定！ fd=\(lat1,lon1,lat2,lon2,R){ R*acos( sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos(lon1-lon2) ) } fd(pi/5,0,pi/3,pi/2,1) http://rio2016.5ch.net/test/read.cgi/math/1723152147/463
466: １３２人目の素数さん [sage] 2024/12/22(日) 21:03:18.86 ID:rWRciXK9 spherical_triangle_area <- function(R, a, b, c) { # 直交座標への変換 to_cartesian <- function(lat, lon) { c(cos(lat) * cos(lon), cos(lat) * sin(lon), sin(lat)) } p <- to_cartesian(a[1], a[2]) q <- to_cartesian(b[1], b[2]) r <- to_cartesian(c[1], c[2]) # 辺のベクトルを計算 u_ab <- q - p u_bc <- r - q u_ca <- p - r # 内角を計算 angle_between <- function(u, v) { acos(sum(u * v) / (sqrt(sum(u^2)) * sqrt(sum(v^2)))) } alpha <- angle_between(u_ab, -u_ca) beta <- angle_between(u_bc, -u_ab) gamma <- angle_between(u_ca, -u_bc) # 超過角を計算 E <- alpha + beta + gamma - pi # 球面三角形の面積を計算 R^2 * E } http://rio2016.5ch.net/test/read.cgi/math/1723152147/466
ﾒﾓ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.028s