23/07/27 20:25:40.02 uVmATNWsa.net
尿瓶ジジイ、素人からも素人臭いと言われるw
323:132人目の素数さん:[sage]:2023/07/27(木) 20:23:10.61 ID:n4ZdqQfJ
>>327
残念なのは、都合の悪いことは全部スルーするお前のオツムだろ
臨床ネタとかいう表現も何か素人臭いし
同業者ってのもお前の自演なんじゃねーの
680:卵の名無しさん (アウアウウー Sabd-TCAR [106.128.68.47])
23/07/27 20:34:17.64 uVmATNWsa.net
神のお告げとかいいながら速攻で間違いを指摘されてダンマリ決め込む尿瓶最高に傑作だわw
681:卵の名無しさん (JP 0Hd3-yMPU [217.138.252.98 [上級国民]])
23/07/28 07:30:51.85 r9+Sd36JH.net
ans=NULL
for(a in 100001:109999){
if((a^2-a) %%10000==0) ans=c(ans,a)
}
ans
682:卵の名無しさん (JP 0Hd3-yMPU [217.138.212.98 [上級国民]])
23/07/28 14:15:09.98 UlP+1rs6H.net
三次方程式の数値解をRで算出する。
cubic_equation =\(a,b,c,d){ #ax³+bx²+cx+d=0
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('i','1i',x)
f7=\(x) gsub('^.=','',x)
str="a!=0, x = (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)/(3 2^(1/3) a) - (2^(1/3) (3 a c - b^2))/(3 a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7() -> x1
x1[2] |> str2lang() |> eval() -> x1
str="a!=0, x = -((1 - i sqrt(3)) (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 + i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7() -> x2
x2[2] |> str2lang() |> eval() -> x2
str="a!=0, x = -((1 + i sqrt(3)) (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 - i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7() -> x3
x3[2] |> str2lang() |> eval() -> x3
c(x1,x2,x3)
}
683:卵の名無しさん (JP 0Hd3-yMPU [217.138.212.98 [上級国民]])
23/07/28 14:28:43.20 UlP+1rs6H.net
Wolframの出力をRで使えるように変換する。
Wolfram2R =\(str,n){ # n : n番目の変数
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('i','1i',x)
f7=\(x) gsub('^.=','',x)
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7() -> x
x[n] |> str2lang()
}
str="a!=0, x = (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)/(3 2^(1/3) a) - (2^(1/3) (3 a c - b^2))/(3 a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
Wolfram2R(str,2)
str="a!=0, x = -((1 - i sqrt(3)) (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 + i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
Wolfram2R(str,2)
str="a!=0, x = -((1 + i sqrt(3)) (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 - i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
Wolfram2R(str,2)
684:卵の名無しさん (JP 0Hd3-yMPU [217.138.212.98 [上級国民]])
23/07/28 17:48:39.94 UlP+1rs6H.net
OutCircle3d <- function(A,B,C){ # 三次元で三角形ABCの外心を返す
dot=function(x,y) sum(x*y)
sss=\(x) sqrt(sum(x^2))
a=sss(B-C)
b=sss(C-A)
c=sss(A-B)
cosA=dot(C-A,B-A)/(b*c)
cosB=dot(A-B,C-B)/(c*a)
cosC=dot(B-C,A-C)/(c*a)
center=(a*cosA*A+b*cosB*B+c*cosC*C)/(a*cosA+b*cosB+c*cosC)
radius=sss(center-A)
c(center=center,radius=radius)
}
A=c(1,2,3)
B=c(4,5,6)
C=c(3,2,1)
OutCircle3d(A,B,C)
685:卵の名無しさん (JP 0H52-RvK4 [217.138.212.98 [上級国民]])
23/07/29 02:21:01.72 qmAKGzOBH.net
OutCircle3d <- function(A,B,C){ # 三次元で三角形ABCの外心を返す
dot=function(x,y) sum(x*y)
sss=\(x) sqrt(sum(x^2))
a=sss(B-C)
b=sss(C-A)
c=sss(A-B)
cosA=dot(C-A,B-A)/(b*c)
cosB=dot(A-B,C-B)/(c*a)
cosC=dot(B-C,A-C)/(c*a)
center=(a*cosA*A+b*cosB*B+c*cosC*C)/(a*cosA+b*cosB+c*cosC)
radius=sss(center-A)
c(x=center[1],y=center[2],z=center[3],radius=radius)
}
incircle3d = \(P,Q,R){ # 三次元で三角形ABCの内心を返す
sss=\(x) sqrt(sum(x^2))
p=sss(Q-R)
q=sss(R-P)
r=sss(P-Q)
s=(p+q+r)/2
S=sqrt(s*(s-p)*(s-q)*(s-r))
radius=S/s
center=(p*P+q*Q+r*R)/(p+q+r)
c(x=center[1],y=center[2],z=center[3],radius=radius)
}
Orthocenter3d=\(A,B,C){ # 三次元で三角形ABCの垂心を返す
bac3d=\(b,a,c){
ab=b-a
ac=c
686:-a dot=sum(ab*ac) # |ab|*|ac|*cos(θ) 内積 bac=acos(dot/sqrt(sum(ab^2))/sqrt(sum(ac^2))) } Aa=bac3d(B,A,C) Ba=bac3d(C,B,A) Ca=bac3d(A,C,B) center=(tan(Aa)*A+tan(Ba)*B+tan(Ca)*C)/(tan(Aa)+tan(Ba)+tan(Ca)) c(x=center[1],y=center[2],z=center[3]) }
687:卵の名無しさん (JP 0H52-RvK4 [217.138.212.98 [上級国民]])
23/07/29 02:35:44.76 qmAKGzOBH.net
excenter3d <- function(A,B,C){ # 三次元で三角形の頂点の座標から傍心の座標と傍心円の半径を返す
sss=\(x) sqrt(sum(x^2))
a=sss(B-C)
b=sss(C-A)
c=sss(A-B)
s=(a+b+c)/2
S=sqrt(s*(s-a)*(s-b)*(s-c))
excenters=c((-a*A+b*B+c*C)/(-a+b+c),(a*A-b*B+c*C)/(a-b+c),(a*A+b*B-c*C)/(a+b-c))
radii=c(2*S/(-a+b+c),2*S/(a-b+c),2*S/(a+b-c))
ans=cbind(matrix(excenters,ncol=3,byrow=TRUE),radii)
colnames(ans)=c('x','y','z','radius')
ans
}
688:卵の名無しさん (JP 0H52-RvK4 [217.138.212.98 [上級国民]])
23/07/29 07:29:38.02 qmAKGzOBH.net
u=pi/180
Vega = c( (18+36/60+56.33635/3600)/24*360*u,(38+47/60+ 1.2802/3600)*u,25.04)
Altair = c( (19+50/60+46.99855/3600)/24*360*u,( 8+52/60+ 5.9563/3600)*u,16.73)
Deneb = c( (20+41/60+25.91514/3600)/24*360*u,(45+67/60+49.2197/3600)*u,1400 )
Stars=rbind(Vega,Altair,Deneb)
dat=cbind(Stars[,3],pi/2- Stars[,2],Stars[,3])
colnames(dat) = c('距離','天頂角','方位角')
dat
colnames(dat)=NULL
P2C=\(x){
r=x[1]
θ=x[2]
φ=x[3]
c(x=r*sin(θ)*cos(θ),y=r*sin(θ)*sin(φ),z=r*cos(θ))
}
data.P=t(apply(dat,1,P2C))
data.P
Oxyz=OutCircle3d(data.P[1,],data.P[2,],data.P[3]) ; Oxyz # 外心の直交座標、外接円半径
Ixyz=incircle3d(data.P[1,],data.P[2,],data.P[3]) ; Ixyz # 内心の直交座標、内接円半径
C2P=\(xyz){
x=xyz[1]
y=xyz[2]
z=xyz[3]
r=sqrt(x^2+y^2+z^2)
θ=acos(z/r)
φ=sign(y)*acos(x/sqrt(x^2+y^2))
c(距離=r,赤経=φ/u,赤緯=(pi/2-θ)/u)
}
C2P(Oxyz) # 外心の極座標
C2P(Ixyz) # 内心の極座標
689:卵の名無しさん (JP 0H52-RvK4 [217.138.212.98 [上級国民]])
23/07/29 07:33:54.33 qmAKGzOBH.net
夏の大三角を形成する星の諸元を以下の通りとする(Wikipediaより)
こと座のベガ(織姫星)
赤経 18h 36m 56.33635s
赤緯 +38° 47′ 01.2802″
距離 25.04 ± 0.07 光年
わし座のアルタイル(彦星)
赤経 19h 50m 46.99855
赤緯 +08° 52′ 05.9563″
距離 16.73 ± 0.05 光年
はくちょう座のデネブ
赤経 (RA, α) 20h 41m 25.91514s
赤緯 (Dec, δ) +45° 16′ 49.2197″
距離 約 1400 光年
問題
(1) 夏の大三角の外接円、内接円の半径を求めよ。
(2) 外心、内心の赤経、赤緯、距離を求めよ。
690:卵の名無しさん (アウアウウー Sac5-53/D [106.128.68.193])
23/07/29 08:39:38.20 HYZ99Dvwa.net
アンタが偽医者でほんとによかったよだってさw
691:卵の名無しさん (スッププ Sda2-HzX1 [49.105.91.210])
23/07/29 11:30:06.37 VAJyuymVd.net
>>659-660
全然、医者板に相応しくない書き込みだな
その書き込みこんなとこで書いて何したいんだ?誰も答えるわけないと思うし
数学板でやってもらっていいっすか?
692:卵の名無しさん (ワッチョイ ee12-53/D [153.137.214.121])
23/07/29 11:57:58.73 ANMrxTku0.net
数学板でも相手にされてないからねw
693:卵の名無しさん (アウアウキー Sa09-53/D [182.251.35.46])
23/07/29 16:10:34.39 JbK1n4r2a.net
389:132人目の素数さん:2023/07/29(土) 16:00:18.62 ID:mTiQ7rLF
>>392
バカすぎるわw
天文学的にみてもまったく無意味な問題のための問題。しかも有効数字が揃ってないという欠陥問題。
アホすぎて話にならん。
せめて球面三角の問題にするくらいの知恵すらないバカ。
尿瓶ジジイ、数学板でも素人丸出しの模様
694:卵の名無しさん (アウアウキー Sa09-53/D [182.251.35.46])
23/07/29 21:26:30.89 JbK1n4r2a.net
>>660
替え歌荒らしにすらセンスないってバカにされてて草
695:卵の名無しさん (ブーイモ MMe6-4qtN [133.159.149.113])
23/07/30 08:16:00.03 4FPQnO7UM.net
# 移動i 経度k を 直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
696:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 11:45:48.69 Jn2avB5O0.net
# 移動i 経度k を直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R=1) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
US=c(47.15,-122.48)
RU=c(55.83,37.62)
CN=c(39.93,116.28)
# ミサイルの走行距離の和
missile=\(i,k){
ik2arc(i,k,US[1],US[2]) + ik2arc(i,k,CN[1],CN[2]) + ik2arc(i,k,CN[1],CN[2])
}
ARC=rbind(US,RU,CN)
missile(colMeans(ARC)[1],colMeans(ARC)[2])
optim(colMeans(ARC),\(ik) missile(ik[1],ik[2]))
697:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 11:56:41.42 Jn2avB5O0.net
# 移動i 経度k を直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
US=c(47.15,-122.48)
RU=c(55.83,37.62)
CN=c(39.93,116.28)
R=6378137
# ミサイルの走行距離の和
missile=\(i,k){
ik2arc(i,k,US[1],US[2],R) + ik2arc(i,k,CN[1],CN[2],R,) + ik2arc(i,k,CN[1],CN[2],R)
}
ARC=rbind(US,RU,CN)
# 平均値での値
missile(colMeans(ARC)[1],colMeans(ARC)[2],R)
# 平均値を初期値にNelder-Mead 法で探索
optim(colMeans(ARC),\(ik) missile(ik[1],ik[2]))
698:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 11:56:55.75 Jn2avB5O0.net
# 移動i 経度k を直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
US=c(47.15,-122.48)
RU=c(55.83,37.62)
CN=c(39.93,116.28)
R=6378137
# ミサイルの走行距離の和
missile=\(i,k){
ik2arc(i,k,US[1],US[2],R) + ik2arc(i,k,CN[1],CN[2],R,) + ik2arc(i,k,CN[1],CN[2],R)
}
ARC=rbind(US,RU,CN)
# 平均値での値
missile(colMeans(ARC)[1],colMeans(ARC)[2],R)
# 平均値を初期値にNelder-Mead 法で探索
optim(colMeans(ARC),\(ik) missile(ik[1],ik[2]))
699:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 11:58:35.94 Jn2avB5O0.net
# 緯度i 経度k を直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
US=c(47.15,-122.48)
RU=c(55.83,37.62)
CN=c(39.93,116.28)
R=6378137
# ミサイルの走行距離の和
missile=\(i,k){
ik2arc(i,k,US[1],US[2],R) + ik2arc(i,k,CN[1],CN[2],R,) + ik2arc(i,k,CN[1],CN[2],R)
}
ARC=rbind(US,RU,CN)
# 平均値での値
missile(colMeans(ARC)[1],colMeans(ARC)[2],R)
# 平均値を初期値にNelder-Mead 法で探索
optim(colMeans(ARC),\(ik) missile(ik[1],ik[2]))
700:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 11:58:45.13 Jn2avB5O0.net
地球を完全な球と仮定する。
仮想敵国の首都と緯度経度は以下の通りである。
ワシントン北緯47.15° 西経122.48°
モスクワ 北緯55.83° 東経37.62°
北京 北緯39.93° 東経116.28°
3つの首都を爆撃するためにミサイル基地を建設する。
(1) 3首都へのミサイルの飛行距離の総和が最小になる基地の緯度経度を求めよ。
(2) (1)で求めた基地はどこの国の領土かを答えよ。
ミサイルの飛行する高度は考慮しなくてよい。
701:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 12:04:02.34 Jn2avB5O0.net
# 緯度i 経度k を直交座標化
ik2xyz=\(i,k,r=1) r*c(cos(i)*cos(k),cos(i)*sin(k),sin(i))
# 緯度ai 経度ak と緯度bi 経度bkの間の球面上距離
ik2arc=\(ai,ak,bi,bk,R) R*acos(sum(ik2xyz(ai,ak,1)*ik2xyz(bi,bk,1)))
US=c(47.15,-122.48)
RU=c(55.83,37.62)
CN=c(39.93,116.28)
R=6378137
# ミサイルの走行距離の和
missile=\(i,k){
ik2arc(i,k,US[1],US[2],R) + ik2arc(i,k,CN[1],CN[2],R) + ik2arc(i,k,CN[1],CN[2],R)
}
ARC=rbind(US,RU,CN)
# 平均値を初期値にNelder-Mead 法で探索
optim(colMeans(ARC),\(ik) missile(ik[1],ik[2]))
702:卵の名無しさん (ワッチョイ 61c7-4qtN [124.47.76.122])
23/07/30 12:06:21.59 Jn2avB5O0.net
par
[1] 52.49637 15.74904
$value
[1] 5975135
$counts
function gradient
141 NA
$convergence
[1] 0
$message
NULL
703:卵の名無しさん (アウアウキー Sa09-53/D [182.251.35.46])
23/07/30 13:41:31.74 aqgQ3jdVa.net
尿瓶まだやってたのか
どうせ誰にも相手にされないのにw
704:岩田純一は犯罪者
23/07/30 17:18:23.78 UXRCimkHk
■■■■■■■■■■■■■■■■■■■■■■■■■■■■
岩田純一は、A医師に暴力をふるった。それを見ていた証人は、多数いる。
磯辺医師、池之上医師、小室医師、関医師、半井英夫医師など。
岩田純一に殴られたA医師が、本富士警察署に訴えてやると言ったところ、
半井英夫が、そんなことをしたら、A医師をクビにすると強迫した。
すべて実名なので、警視庁サイバー犯罪対策課に 通報してください。
■■■■■■■■■■■■■■■■■■■■■■■■■■■■
705:卵の名無しさん (ブーイモ MMe6-4qtN [133.159.149.218])
23/07/30 21:40:16.20 f5OP0fsKM.net
(ai,ak) -> (pi/2,0) -> (0,0,1)
(bi,bk) -> (pi/2-(ai-bi),bk-ak) -> sin(ai-bi)cos(bk-ak),sin(ai-bi)sin(bk-ak),cos(ai-bi)
(ci,ck) -> (pi/2-(ai-ci),ck-ak) -> sin(ai-ci)cos(ck-ak),sin(ai-ci)sin(ck-ak),cos(ai-ci)
(sin(ai-bi)cos(bk-ak)sin(ai-ci)cos(ck-ak)+sin(ai-bi)sin(bk-ak)sin(ai-ci)sin(ck-ak))/(sss*sss)
bac.sph = \(ai,ak,bi,bk,ci,ck) {
((sin(ai-bi)*sin(ai-ci)*( cos(bk-ak)*cos(ck-ak) + sin(bk-ak)*sin(ck-ak) ))
/( abs(sin(ai-bi)) *abs(sin(ai-ci)) )) |> acos()
}
706:卵の名無しさん (ワッチョイ ee12-53/D [153.137.214.121])
23/07/31 10:50:23.43 WsDDKQKG0.net
尿瓶、早朝発狂なしw
707:卵の名無しさん (アウアウキー Sa09-53/D [182.251.35.46])
23/07/31 11:04:11.56 1ZcxZ36Da.net
911:卵の名無しさん:[sage]:2023/07/29(土) 21:38:30.79 ID:JbK1n4r2
アンポンタン、アンポンタン♪
数学できない偽医者♪
アンポンタン、アンポンタン♪
708:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 14:36:54.05 xgtz5I6M0.net
地球を完全な球とする。
E,W,S,Nを下記の点とする
日本の東西南北端点の経度緯度
区分 場所 世界測地系
経度 緯度
E 最東端 東京都 南鳥島 153°59′12″ 24°16′59″
W 最西端 沖縄県 与那国島 122°55′57″ 24°27′05″
S 最南端 東京都 沖ノ鳥島 136°04′11″ 20°25′31″
N 最北端 北海道 択捉島 148°45′08″ 45°33′26″
地表に沿った弧でE-W,S-Nを結ぶとき2つの弧の交点の経度、緯度を求めよ。
709:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 14:37:34.37 xgtz5I6M0.net
E=c(24+16/60+59/3600,153+59/60+12/3600) ; E
W=c(24+27/60+5/3600,122+55/60+57/3600) ; W
S=c(20+25/60+31/3600,136+4/11) ; S
N=c(45+33/60+26/3600,148+45/60+8/3600) ; N
# uE+(1-u)W == vS+(1-v)N
f = \(u,v) ((u*E+(1-u)*W) - (v*S+(1-v)*N))^2 |> sum()
opt=optim(runif(2),\(x) f(x[1],x[2]),method='L') ; opt
u0=opt$par[1]
v0=opt$par[2]
u0*E+(1-u0)*W
v0*S+(1-v0)*N
710:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 15:09:34.00 xgtz5I6M0.net
# (緯度・経度) 弧A(a1,a2)-B(b1,b2) と 弧C(c1,c2)-D(d1,d2)の交点の緯度経度を算出
intsect.sph = \(a1,a2,b1,b2,c1,c2,d1,d2){
x=(b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))
y=(a1*(b2-d2)+a2*(d1-b1)+b1*d2-b2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))
c(x*a1+(1-x)*b1, x*a2+(1-x)*b2)
}
intsect.sph(E[1],E[2],W[1],W[2],S[1],S[2],N[1],N[2])
711:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 15:10:16.61 xgtz5I6M0.net
> intsect.sph(E[1],E[2],W[1],W[2],S[1],S[2],N[1],N[2])
[1] 24.36805 138.30719
Nelder-Meadで検算
> # uE+(1-u)W == vS+(1-v)N
> f = \(u,v) ((u*E+(1-u)*W) - (v*S+(1-v)*N))^2 |> sum()
> opt=optim(runif(2),\(x) f(x[1],x[2]),method='L-BFGS-B') ; opt
$par
[1] 0.4950927 0.8431172
$value
[1] 8.582807e-28
$counts
function gradient
6 6
$convergence
[1] 0
$message
[1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
> u0=opt$par[1]
> v0=opt$par[2]
> u0*E+(1-u0)*W
[1] 24.36805 138.30719
712:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 18:01:12.58 xgtz5I6M0.net
"
地球を完全な球とする。
E,W,S,Nを下記の点とする
日本の東西南北端点の経度緯度
緯度 経度
E 最東端 東京都 南鳥島 24.28306 153.9867
W 最西端 沖縄県 与那国島 24.45139 122.9325
S 最南端 東京都 沖ノ鳥島 20.42528 136.3636
N 最北端 北海道 択捉島 45.55722 148.7522
日本の領土領海の球面四角形ESWNの中に点Pをとり
4つの球面三角形PES,PSW,PWN,PNEの面積が等しくなるように4分割して
米英中露が分割統治する。
点Pの緯度・経度を求めよ。
"
E
S
W
triangle.sph(E[1],E[2],S[1],S[2],W[1],W[2])
M=colMeans(rbind(E,S,W,N))
i=M[1]
k=M[2]
triangle.sph(i,k,E[1],E[2],S[1],S[2])
ABC2S.sph=\(A,B,C,r=1){
r^2*(sum(triangle.sph(A[1],A[2],B[1],B[2],C[1],C[2]))-pi)
}
ABC2S.sph(M,E,S)
fn=\(i,k){
P=c(i,k)
PES=ABC2S.sph(P,E,S)
PSW=ABC2S.sph(P,S,W)
PWN=ABC2S.sph(P,W,N)
PNE=ABC2S.sph(P,N,E)
diff(c(PES,PSW,PWN,PNE))^2 |> sum()
}
opt=optim(M,\(x) fn(x[1],x[2])) ; opt
713:卵の名無しさん (ワッチョイ 3224-RvK4 [149.50.210.17 [上級国民]])
23/07/31 18:37:44.96 xgtz5I6M0.net
>>679
答が海上だと今ひとつ興味がわかないので数値を変えた練習問題。
四則演算で答が出せる。
A:嘉手納飛行場 Kadena Air Base 北緯26度21分06秒 東経127度46分10秒
B:三沢飛行場(三沢空港) Misawa Air Base (Misawa Airport)北緯40度42分19秒 東経141度22分19秒
C:朝鮮の東倉里という場所にあるミサイル基地 経緯度: 北緯39度39分48秒 東経124度42分8秒
D:内閣総理大臣官邸 北緯35度40分23.15秒 東経139度44分32.3秒
(1)地球を完全な球として大円の弧A-Bと大円の弧C-Dの交わる点の緯度経度を計算せよ。
(2)その地点のある市町村名を答えよ。
714:卵の名無しさん (アウアウキー Sa09-53/D [182.251.21.101])
23/07/31 22:44:42.41 xVkK+rZea.net
188:卵の名無しさん:[sage]:2023/07/31(月) 22:35:21.47 ID:6IPqu6HL
ネタはまだストックがありますか?脳内妄想医さん
715:卵の名無しさん (ブーイモ MMe6-RvK4 [133.159.149.109 [上級国民]])
23/08/01 08:28:04.29 oRgPBjXSM.net
a=c(replicate(65,sum(sample(9)*10^(0:8))),sum(sample(8)*10^(0:7)),replicate(34,sum(sample(9)*10^(0:8))))
paste(a,collapse='+') |> cat()
716:卵の名無しさん (ブーイモ MMb9-PQn6 [220.156.12.242])
23/08/01 14:43:21.17 RSC+6NqWM.net
a=4
b=5
c=6
A=abc2ABC(a,b,c)[1]
B=0i
C=a+0i
f=\(xy){
P=xy[1]+1i*xy[2]
if(!in3(P,A,B,C)) return(NA)
sum(diff(c(ABC2S(A,B,P), ABC2S(B,C,P),ABC2S(C,A,P)))^2)
}
optim(c(Re(A+B+C)/3,Im(A/3)),f)
717:卵の名無しさん (JP 0H52-RvK4 [217.138.252.98 [上級国民]])
23/08/02 06:43:27.14 3NNjYFcbH.net
f= \(a1,a2,b1,b2,c1,c2,d1,d2){
lat=((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1)))*a1+(1-((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))))*b1
lon=((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1)))*a2+(1-((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))))*b2
c(lat,lon)
}
f(E[1],E[2],W[1],W[2],S[1],S[2],N[1],N[2])
718:卵の名無しさん (JP 0H52-RvK4 [217.138.252.98 [上級国民]])
23/08/02 06:57:25.81 3NNjYFcbH.net
f= \(a1,a2,b1,b2,c1,c2,d1,d2){
lat=((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1)))*a1+(1-((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))))*b1
lon=((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1)))*a2+(1-((b1*(d2-c2)+b2*(c1-d1)-c1*d2+c2*d1)/((a1-b1)*(c2-d2)+a2*(d1-c1)+b2*(c1-d1))))*b2
c(lat,lon)
}
E=c(24+16/60+59/3600,153+59/60+12/3600) ; E
W=c(24+27/60+5/3600,122+55/60+57/3600) ; W
S=c(20+25/60+31/3600,136+4/11) ; S
N=c(45+33/60+26/3600,148+45/60+8/3600) ; N
f(E[1],E[2],W[1],W[2],S[1],S[2],N[1],N[2])
# optimでの検証
# uE+(1-u)W == vS+(1-v)N
f = \(u,v) ((u*E+(1-u)*W) - (v*S+(1-v)*N))^2 |> sum()
opt=optim(runif(2),\(x) f(x[1],x[2]),method='L-BFGS-B') ; opt
u0=opt$par[1]
v0=opt$par[2]
u0*E+(1-u0)*W
v0*S+(1-v0)*N
719:卵の名無しさん (ワッチョイ 0958-PQn6 [14.13.16.0])
23/08/02 14:05:14.09 LakLPy4e0.net
妻がとうもろこを買ってきてくれたので
バターコーンの低温調理をやってみる。
80℃30分に設定。
バターレモンチキンが美味だったのでバターレモンコーンにアレンジ。
720:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/02 15:09:29.63 3XtkjKmWH.net
library(RcppAlgos)
f=\(n) (apply(permuteGeneral(6,n,repetition = TRUE),1,\(x) sum(x)%%7==0)) |> mean()
y=sapply(1:9,f)
plot(y)
721:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/03 07:39:11.94 5oQ9dj8SH.net
on.sphere =\(n){ # 球面上から一様分布でn個選んだ地点の緯度・経度を返す
ans=NULL
z=runif(n,-1,1)
t=runif(n,-pi,pi)
x=sqrt(1-z^2)*cos(t)
y=sqrt(1-z^2)*sin(t)
i=asin(z)
k=sign(y)*acos(x/sqrt(x^2+y^2))
cbind(i=i,k=k)
}
722:卵の名無しさん (ワッチョイ 0958-PQn6 [14.13.16.0])
23/08/03 15:52:50.76 t4ZfvwI30.net
P = \(n){
Grid <- function(v,n){ # Grid(1:6,3) == expand.grid(1:6,1:6,1:6)
paste0("expand.grid(",paste0(rep(deparse(v),n),collapse=","),')') |> str2lang() |> eval()
}
723:d=Grid(1:6,n) apply(d,1,\(x) sum(x)%%7==0) |> mean() } sapply(1:7,P) |> MASS::fractions()
724:卵の名無しさん (スッププ Sda2-Xleb [49.105.85.56])
23/08/03 17:32:32.22 69uB5F4Yd.net
尿瓶ジジイも自分の分をわきまえたのか
間違った脳内医療をここに晒すことはなくなったなwww偉いじゃん
尿瓶ジジイは内視鏡スレには現れるのに、消化器専門医スレには現れない
何故なんでしょうねwww
ヒント:内視鏡はコメディカルでもある程度理解できます!
725:卵の名無しさん (ワッチョイ 0958-PQn6 [14.13.16.0])
23/08/03 17:39:07.97 t4ZfvwI30.net
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、
恥ずかしくね?
726:卵の名無しさん (ワッチョイ 0958-PQn6 [14.13.16.0])
23/08/03 17:40:53.60 t4ZfvwI30.net
クローンやUCとかやりたくないね。
書類の記載も面倒くさいし。
727:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 06:17:13.22 LH6ITclIH.net
r=C/2
R=(r^2+h^2)/(2*h)
theta=2*asin(r/R)
R^2*theta/2 - r*(R-h)
bow = \(C,h){
R=((C/2)^2+h^2)/(2*h)
S=(-4*C^3*h+(C^2+4*h^2)^2*asin((4*C*h)/(C^2+4*h^2))+16*C*h^3)/(64*h^2)
# S=(((C/2) ^2+h^2)/(2*h))^2*asin((C/2)/(((C/2) ^2+h^2)/(2*h)))-(C/2)*((((C/2) ^2+h^2)/(2*h))-h)
c(S=S,R=R)
}
728:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 06:47:37.99 LH6ITclIH.net
(((r^2+h^2)/(2*h))^2*(atan((C/2)/((((C/2)^2+h^2)/(2*h))-h)))/2-((r^2+h^2)/(2*h))*r/2)*2+C*h/2
729:卵の名無しさん (ブーイモ MMe6-NUl3 [133.159.148.83])
23/08/04 07:59:46.23 kNnoYEiTM.net
R: 大元の球半径
r:生地の底面半径
r0:生地の切離側面の円の半径
r0=sqrt(R^2-x^2)
h0=r0-(R-h) # 弓形側面の高さ
C0=2*sqrt(r^2-x^2) # 弦の長さ
730:卵の名無しさん (ブーイモ MMe6-NUl3 [133.159.148.83])
23/08/04 11:53:43.14 kNnoYEiTM.net
bow = \(C,h){
R=((C/2)^2+h^2)/(2*h)
S=(-4*C^3*h+(C^2+4*h^2)^2*asin((4*C*h)/(C^2+4*h^2))+16*C*h^3)/(64*h^2)
# S=(((C/2) ^2+h^2)/(2*h))^2*asin((C/2)/(((C/2) ^2+h^2)/(2*h)))-(C/2)*((((C/2) ^2+h^2)/(2*h))-h)
c(S=S,R=R)
}
calc = \(x,C=10,h=3){
r=C/2
R=(r^2+h^2)/(2*h)
f=\(x){
r0=sqrt(R^2-x^2)
h0=r0-(R-h)
C0=2*sqrt(r^2-x^2)
bow(C0,h0)['S']
}
integrate(f,0,x)$value
}
calc(10)
731:すみちゃん
23/08/04 15:54:58.30 gVuTFyJfV
美容液(シートマスク)というものは、栄養成分の有効性が高いものほど価格も高額になるのが普通です。予算に添って長らく愛用することになっても金銭的に大変にならない程度のものをおすすめします。
URLリンク(ikanakama.ink)
732:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 13:15:17.66 LH6ITclIH.net
bow = \(C,h){ # 弓形の弦の長さC、高さhから半径と面積を返す
R=((C/2)^2+h^2)/(2*h)
S=(-4*C^3*h+(C^2+4*h^2)^2*asin((4*C*h)/(C^2+4*h^2))+16*C*h^3)/(64*h^2)
# S=((h^2+r^2)^2*atan((4*C*h)/(C^2-4*h^2))+2*C*h^3-2*h*r*(h^2+r^2))/(4*h^2)
c(S=S,R=R)
}
733:卵の名無しさん (JP 0H96-r+G9 [61.199.225.248])
23/08/04 13:17:55.05 Kgx2eOG2H.net
数学板でも数学できない偽医者とか言われてるw
734:卵の名無しさん (ワッチョイ 0958-O5nL [14.13.16.0 [上級国民]])
23/08/04 16:17:02.78 CPhIRibW0.net
パン作りをしながら思いついた問題、自分への宿題
1次発酵が終わったパン生地を球を平面で切った形状とする。
球冠spherical capと呼ぶらしい。中華饅頭の形のイメージ。
生地の底面の円の半径は10cm、生地の高さは5cmであった。
生地を平行な2つの平面で切って3分割したい。
鉛直方向に切断するとして端から何cmのところで切ればよいか?
735:卵の名無しさん (ワッチョイ 0958-O5nL [14.13.16.0 [上級国民]])
23/08/04 16:17:32.89 CPhIRibW0.net
>>703
東大卒が正解を出せない問題をRで数値解を出しているからね。
736:卵の名無しさん (ワッチョイ 0958-O5nL [14.13.16.0 [上級国民]])
23/08/04 16:18:30.71 CPhIRibW0.net
ようやく計算が終わった。
R: 大元の球半径 (10^2+5^2)/2*5=12.5
r:生地の底面半径 10
h: 生地の高さ 5
r0:生地の切離側面の円の半径
x: [0,r]
r0=sqrt(R^2-x^2) # 切離面円の半径
h0=r0-(R-h) # 弓形側面の高さ
C0=2*sqrt(r^2-x^2) # 弦の長さ
bow = \(C,h){ # 弓形の弦の長さC、高さhから半径と面積を返す
if(C*h==0) return(0)
R=((C/2)^2+h^2)/(2*h)
S=(-4*C^3*h+(C^2+4*h^2)^2*asin((4*C*h)/(C^2+4*h^2))+16*C*h^3)/(64*h^2)
c(S=S,R=R)
}
# xでの球冠の切離面の面積
area.bow=\(x,r=10,h=5){
R=bow(2*r,h)['R']
r0=sqrt(R^2-x^2) # 切離面円の半径
h0=r0-(R-h) # 弓形側面の高さ
C0=2*sqrt(r^2-x^2) # 弦の長さ
bow(C0,h0)[1]
}
area.bow=Vectorize(area.bow)
# [0,x]までの球冠の体積
v02x = \(x) integrate(area.bow,0,x)$value
v02x = Vectorize(v02x)
V=v02x(10)*2 # 球冠の体積
curve(v02x,0,10)
abline(h=V/6,lty=3)
x.ans=uniroot(\(x) v02x(x)-V/6,c(0,10))$root
ans=10-x.ans
ans
737:卵の名無しさん (ワッチョイ 0958-O5nL [14.13.16.0 [上級国民]])
23/08/04 16:25:04.85 CPhIRibW0.net
高校数学スレにはこれが解ける人はいないみたい。
俺はRで数値解が出せたけど。
長さが2,3,4,5,6,7の6本のチョコ棒を組み合わせて最も鋭利な頂点をもつ三角錐を作る。
底面を三角形ABC、頂点をDとし頂点Dが最も鋭利とする。
(1)∠ADB+∠BDC+∠CDAは何度か?
(2)その三角錐の高さを求めよ
(3)その三角錐の体積を求めよ。
(4)その三角錐を図示せよ。
いずれも数値は有効数字3桁でよい。
738:卵の名無しさん (スッップ Sda2-r+G9 [49.98.212.202])
23/08/04 17:15:24.66 qDvgyvpNd.net
素人相手にも脳内東大合格者の偽医者とか言われてて草
739:卵の名無しさん (スッップ Sda2-r+G9 [49.98.212.202])
23/08/04 17:19:17.79 qDvgyvpNd.net
>>707
シカトされてるだけだろアホが
740:卵の名無しさん (ワッチョイ 0958-5v8d [14.13.16.0])
23/08/04 17:19:38.96 CPhIRibW0.net
>>708
理Iくらい普通は通るだろ。
741:卵の名無しさん (ワッチョイ 0958-5v8d [14.13.16.0])
23/08/04 17:21:14.90 CPhIRibW0.net
>>709
難問に取り組んで解答をする人もいるんだよ。
433 132人目の素数さん 2023/07/30(日) 02:02:36.59 ID:Fyeh8XR9
>>8
(1) 3首都へのミサイルの飛行距離の総和が最小になる基地の緯度経度は、北緯48.77°、東経88.69°です。この場合、ミサイルの飛行距離の総和は約1.57×10^7 mとなります。この値は、地球の半径を6.378×10^6 mとして、大円距離の公式を用いて求めました。
(2) (1)で求めた基地は、モンゴル国の領土内にあります。具体的には、ホブド県の南西部に位置します。
742:卵の名無しさん (ワッチョイ 0958-5v8d [14.13.16.0])
23/08/04 17:22:41.28 CPhIRibW0.net
>>707
これは手書き計算では難しいと思う。
俺はRを使って数値解が出せた。
743:卵の名無しさん (ワッチョイ 0958-5v8d [14.13.16.0])
23/08/04 17:25:13.97 CPhIRibW0.net
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
744:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 17:44:48.93 LH6ITclIH.net
プログラムできる能力と根気がないと解が出せないからね。それなりに時間もかかる。
概略を考えながら寝て朝に実装。朝飯前の問題シリーズ。
745:卵の名無しさん (アウアウウー Sac5-MkdW [106.128.69.32])
23/08/04 22:27:25.74 KsZ3Q2aga.net
早朝から発狂かよ
クソコテ以外には荒らしにもゴミ扱いみたいだね
そのクソコテにも脳内だと思われてるあり様
746:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 22:46:58.93 LH6ITclIH.net
fn=\(p){
P=p+1i*p^2
B=(p^2+1/2)*1i # y-p^2 =-1/(2*p)(x-p)
C=0i
H=ABC2H(P,B,C)$H
c(Re(H),Im(H))
}
fn=Vectorize(fn)
x=seq(-5,5,le=10001)
H=t(fn(x))
plot(H,type='l',xlab='x',ylab='y',asp=1,col='red')
curve(x^2,add=TRUE,lty=3)
747:卵の名無しさん (JP 0H4d-O5nL [146.70.31.52 [上級国民]])
23/08/04 22:53:36.97
748:LH6ITclIH.net
749:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 06:26:51.63 T/2uNk4DH.net
Wolfram2R = \(str){ # Wolframの出力をRで使えるように変換
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('^.=','',x)
f7=\(x) gsub('i','1i',x)
f8=\(x) gsub('s1in','sin',x)
f9=\(x) gsub('s1in^(-1)','asin',x)
f10=\(x) gsub('tan^(-1)','atan',x)
f11=\(x) gsub('cos^(-1)','acos',x)
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7()
}
750:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 06:39:30.79 T/2uNk4DH.net
朝のお勤めw
Wolframのプレーンテキスト出力をRで使える形式に変換するスクリプトを作成使用していたが、
Wolframはatanをtan^(-1)と表示するのでcotと紛らわしい。
これをatanに変換するようにデバッグ。
sin^(-1)のようにiが含まれるとiが虚数単位なのか否かを判断しなくちゃならんことにも気づいた。デバックした。
751:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 07:07:53.65 T/2uNk4DH.net
Wolfram2R = \(str){ # Wolframの出力をRで使えるように変換
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('^.=','',x)
f7=\(x) gsub('sin\\^\\(-1\\)','asin',x)
f8=\(x) gsub('tan\\^\\(-1\\)','atan',x)
f9=\(x) gsub('cos\\^\\(-1\\)','acos',x)
f10=\(x) gsub('i','1i',x)
f11=\(x) gsub('s1in','sin',x)
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |> f7() |> f8() |> f9() |> f10() |> f11()
}
752:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 07:09:02.21 T/2uNk4DH.net
動作確認
x="-((1 - i sqrt(3)) (sqrt((-27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 + i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
Wolfram2R(x)
e=exp(1)
a=b=c=d=1
Wolfram2R(x) |> str2lang() |> eval()
> Wolfram2R(x) |> str2lang() |> eval()
[1] -0.1192194+0.8967725i
753:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 07:12:50.01 T/2uNk4DH.net
Wolframの数値と照合
{x == -1/3 + (2/(-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2]))^(1/3)/3 + (5 (2/(-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2]))^(2/3))/3 - (9 ArcTan[1/(Sqrt[1 - E^(-2)] E)])/(2 2^(1/3) (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2])^(2/3)) - Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2]/(6 2^(1/3) (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2])^(2/3)) + I ((2/(-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2]))^(1/3)/Sqrt[3] - (5 (2/(-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2]))^(2/3))/Sqrt[3] + (9 Sqrt[3] ArcTan[1/(Sqrt[1 - E^(-2)] E)])/(2 2^(1/3) (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2])^(2/3)) + Sqrt[(32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2)/3]/(2 2^(1/3) (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)] + Sqrt[32 + (-20 + 27 ArcTan[1/(Sqrt[1 - E^(-2)] E)])^2])^(2/3)))}
x≈-0.11922 + 0.89677 i
今のところ、きちんと動作している。
まあ、sin以外にiを含む単語がでてくると誤変換するだろうけど。
754:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 07:15:36.40 T/2uNk4DH.net
Wolfram2R = \(str){ # Wolframの出力をRで使えるように変換
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('^.=','',x)
f7=\(x) gsub('sin\\^\\(-1\\)','asin',x)
f8=\(x) gsub('tan\\^\\(-1\\)','atan',x)
f9=\(x) gsub('cos\\^\\(-1\\)','acos',x)
f10=\(x) gsub('i','1i',x)
f11=\(x) gsub('s1in','sin',x)
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |>
f7() |> f8() |> f9() |> f10() |> f11()
}
755:卵の名無しさん (JP 0Hd5-swBy [146.70.31.52 [上級国民]])
23/08/05 09:21:55.73 T/2uNk4DH.net
Wolfram2R = \(str){ # Wolframの出力をRで使えるように変換
f1=\(x) gsub('^ ','',x)
f2=\(x) gsub(' = ','=',x)
f3=\(x) gsub(' \\+ ','\\+',x)
f4=\(x) gsub(' \\- ','\\-',x)
f5=\(x) gsub(' ','*',x)
f6=\(x) gsub('^.=','',x)
f7=\(x) gsub('sin\\^\\(-1\\)','asin',x)
f8=\(x) gsub('tan\\^\\(-1\\)','atan',x)
f9=\(x) gsub('cos\\^\\(-1\\)','acos',x)
f10=\(x) gsub('i','1i',x)
f11=\(x) gsub('s1in','sin',x)
f12=\(x) gsub('sgn','sign',x)
strsplit(str,',') |> unlist() -> str
str |> f1() |> f2() |> f3() |> f4() |> f5() |>f6() |>
f7() |> f8() |> f9() |> f10() |> f11() |> f12() -> ans
cat(ans,'\n')
invisible(ans)
}
x="-((1 - i sqrt(3)) (sqrt((-27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^(1/3))/(6 2^(1/3) a) + ((1 + i sqrt(3)) (3 a c - b^2))/(3 2^(2/3) a (sqrt((-27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^2 + 4 (3 a c - b^2)^3) - 27 a^2 d + 27 a^2 sin^(-1)(1/e) + 9 a b c - 2 b^3)^(1/3)) - b/(3 a)"
Wolfram2R(x)
e=exp(1)
a=b=c=d=1
Wolfram2R(x) |> str2lang() |> eval()
756:卵の名無しさん (スプープ Sd22-lp8H [1.73.148.73])
23/08/05 19:32:38.88 iO3zDLfFd.net
尿瓶ジジイダンマリかよ
757:卵の名無しさん (JP 0Hf8-swBy [217.138.212.122 [上級国民]])
23/08/05 20:12:53.90 xeUk+Em+H.net
A=c(0,0,0)
B=c(0,0,1)
f=\(x){
c1=x[1]
c2=x[2]
d1=x[3]
d2=x[4]
if(c1< -pi/2 | pi/2<c1 | d1< -pi/2 | pi/2<d1) return(NA)
if(c2< -pi | pi<c2 | d2< -pi | pi <d2) return(NA)
C=c(cos(c1)*cos(c2),cos(c1)*sin(c2),sin(c1))
D=c(cos(d1)*cos(d2),cos(d1)*sin(d2),sin(d1))
P=B/2+C/2
Q=C/2+D/2
R=D/2+B/2
((A-P) %.% (B-C))^2+((A-Q) %.% (C-D))^2+((A-R) %.% (D-B))^2
}
opt=optim(runif(4,-pi/2,pi/2),f) ; opt
C=c(cos(c1)*cos(c2),cos(c1)*sin(c2),sin(c1)) ; C
D=c(cos(d1)*cos(d2),cos(d1)*sin(d2),sin(d1)) ; D
(B-C)^2 |> sum() |> sqrt()
(C-D)^2 |> sum() |> sqrt()
(D-B)^2 |> sum() |> sqrt()
ABCD=rbind(A,B,C,D)
library(r
758:gl) plot3d(ABCD,type='s',col='blue',size=1,xlab='x',ylab='y',zlab='z') aspect3d('iso') segments3d(rbind(A,C),col='blue',lwd=2) segments3d(rbind(A,B),col='blue',lwd=2) segments3d(rbind(A,D),col='blue',lwd=2) segments3d(rbind(B,C),col='blue',lwd=2) segments3d(rbind(C,D),col='blue',lwd=2) segments3d(rbind(D,B),col='blue',lwd=2)
759:卵の名無しさん (ワッチョイ 22c9-+YZw [59.139.198.126])
23/08/06 05:40:16.73 8o4WLPlZ0.net
チンパンジーが計算ドリルを電卓叩いてキーキー喚いてるのは滑稽極まりない
760:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/06 05:41:32.76 8CLnscyb0.net
abs = \(z){
if(is.complex(z)){
return(sqrt(Re(z)^2+Im(z)^2))
}else{
return(sqrt(sum(z^2)))
}
}
abs(cos(1)+1i*sin(1))
abs(c(cos(1)*cos(2),cos(1)*sin(2),sin(1)))
761:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/06 05:42:46.20 8CLnscyb0.net
高校数学の質問スレで
>じゃあネットで答えを聞いてもいいのか?
>それのどこが高校数学だよ
ですって
これほどの馬鹿がいるだろうか?
762:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 05:54:02.06 sYHzVRqAa.net
尿瓶ジジイの数学板でのチンパンごっこ()は数学でも算数でもないからトンチンカンなレスしかつかないどころかみんなに馬鹿にされるんだろw
763:卵の名無しさん (ワッチョイ 22c9-+YZw [59.139.198.126])
23/08/06 05:55:15.25 8o4WLPlZ0.net
911:卵の名無しさん:[sage]:2023/07/29(土) 21:38:30.79 ID:JbK1n4r2
アンポンタン、アンポンタン♪
数学できない偽医者♪
アンポンタン、アンポンタン♪
764:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:06:53.06 cNr8zxB+M.net
>>727
R って電卓じゃないんだけど。
臨床医の嗜みのRを使えないの?
765:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:11:08.43 cNr8zxB+M.net
こういうのを作って置くと便利。数値に応じたグラフも出力。
# ================== gamma by minute ====
gamma2hourlyVol=function(mg,mL,kg,gamma){ # mcg/kg/min -> mL/h
# mg/1000mcg * kg/body * 60min/h * mL/mg * gammma
coef=60/1000*mL/mg*kg
velocity=coef*gamma
mat=cbind(gamma,velocity)
colnames(mat)=c('mcg/kg/min','mL/h')
L=ifelse(length(gamma)==1,0,min(gamma))
U=ifelse(length(gamma)==1,5*gamma,max(gamma))
main=paste('concenteration =',mg,'mg /',mL,'mL\nweight =',kg,'kg')
curve(coef*x,L,U,bty='l',xlab='mcg/kg/min',ylab='mL/h',main=main)
n=21
X=seq(L,U,length=n)
Y=coef*X
for(i in 1:n){
for(j in 1:n){
segments(X[i],Y[1],X[i],Y[n],col=c(8,2)[i%%2+1],lty=3)
segments(X[1],Y[j],X[n],Y[j],col=c(8,2)[j%%2+1],lty=3)
}
}
for( i in 1:nrow(mat)){
text(mat[i,1],mat[i,2], as.character(mat[i,2]))
}
mat
}
766:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:11:17.86 cNr8zxB+M.net
# CV用高濃度
#ノルアドレナリン1mgを全量20mLに希釈して体重50kgで0.05~0.30μg/kg/minで投与するときの速度とグラフを出力
# Noradrenaline 1mg/20mL at 0.05-0.3mcg/kg/min
gamma2hourlyVol(mg=1,mL=20,kg=50,gamma=c(0.05,0.30)) ; legend('top',bty='n',legend='for Central Vein',cex=1.25)
# 末梢用低濃度
#ノルアドレナリン1mgを全量50mLに希釈して体重50kgで0.05~0.30μg/kg/minで投与するときの速度とグラフを出力
# Noradrenaline 1mg/20mL at 0.05-0.3mcg/kg/min
gamma2hourlyVol(mg=1,mL=50,kg=50,gamma=c(0.05,0.30)) ; legend('top',bty='n',legend='for Peripheral Vein',cex=1.25)
767:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 08:14:27.90 sYHzVRqAa.net
>>732
話がまるで通じないチンパンジーにはRも電卓も関係ないだろ
768:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:16:13.40 cNr8zxB+M.net
>>720
パイプなしだと()対応で混乱しそう。
ver4 以後、パイプ表現できるので多用している。
自作できなくはな�
769:「けど '%|%' <- function(x,FUN) FUN(x) # パイプ
770:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:17:43.42 cNr8zxB+M.net
'%.%' <- function(x,y){ # ベクトルの内積 加減計算より優先
if(is.complex(x)|is.complex(y)){
a=c(Re(x),Im(x))
b=c(Re(y),Im(y))
}else{
a=x
b=y
}
sum(a*b)
}
771:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:19:46.07 cNr8zxB+M.net
>>735
んで、R使えないの?
結果を動画にするのはPythonよりも便利。
計算機を自分の手足のように使える人のための練習問題
四面体ABCDは以下の条件を満たす。
(ア)AB=AC=AD=BC
(イ)BC,CD,DBの中点をそれぞれP,Q,Rとすると、AP⊥BCかつAQ⊥CDかつAR⊥DB
この四面体は正四面体といえるか。
Rによる実験
条件を満たす四面体を作図
URLリンク(i.imgur.com)
772:卵の名無しさん (ブーイモ MMcd-lnPD [202.214.231.203])
23/08/06 08:27:25.23 cNr8zxB+M.net
どこの国立を落ちたの?
に答えらないで発狂しているのが尿瓶チンパポンコツタコフェチである。
Q.E.D.
さてはシリツだな?
773:卵の名無しさん (スップー Sd22-fPrT [1.73.22.163])
23/08/06 08:36:59.52 OXNG2/+Fd.net
Rとか糞ソフト使ってる日本人の医者周りにはいないけどなwww
フリーソフトなのにそのcodeが本当に正しいって誰が責任持つんだよ
そんなソフトで正式な論文書けないだろ
RでAcceptされてるのはImpact Factor低いクソ論文だけだよ
尿瓶ジジイは学会とか行ったことないんだろうな
774:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 09:56:34.57 sYHzVRqAa.net
尿瓶ジジイは脳内だから現実の医者が全く興味ないことでマウントとりたくて必死なんだよなww
775:卵の名無しさん (ワッチョイ 22c9-+YZw [59.139.198.126])
23/08/06 09:59:34.97 8o4WLPlZ0.net
>>738
胆汁ドレナージといいRといいとことん脳内医療丸出しだなマジで
776:ミカエル!!世界の終わり!!
23/08/06 16:54:08.29 xi31lNV4o
終末予言 「数年後 魚が絶滅して、漁業が終了します」!!
近未来「岸田総理」が「放射能汚染水」を海に流した!
「神」が怒り「人類の大量絶滅期」(たいりょうぜつめつき)がやってきた!
「放射能汚染水」が「ヨーロッパ」に流れついた!
「ヨーロッパの全国民」が怒り、「日本政府」に
「賠償金」(ばいしょうきん)「9千兆円」を請求(せいきゅう)した!
「日本人」は「世界中の人間から」「憎まれた」!!
「放射能汚染水」が原因で「海の水」が「血のように赤くなった」!
「魚」が全滅した!!魚が滅(ほろ)びた!全部死んだ!!
「「日本の漁業が数年後に終了した」」!!「ヨーロッパの漁業も終了した」!
「岸田総理」は「全人類」から「恨(うら)まれ、憎まれた」!!
ミカエル 「魚が絶滅し「食糧危機」がおこり「人類の3分の2」が死んだ」!
「「人間が人間を食べる時代」がやってくるだろう・・」
「この世界はもうすぐ終わります・・・・」
「神」と「天使」が人類をさばいた!!「七つの大災害」がおこった!!」
777:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 11:00:27.48 BOo6XOy90.net
>>742
どこの国立を落ちたの?
に答えらないで発狂しているのが尿瓶チンパポンコツタコフェチである。
Q.E.D.
さてはシリツだな?
778:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 11:02:22.02 BOo6XOy90.net
>>740
FDAはRの解析でも認可する。
有
779:料ソフトでもバグがあるぞ。 エクセルの統計は酷いので有名だった。
780:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 11:03:42.84 BOo6XOy90.net
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、
恥ずかしくね?
781:卵の名無しさん (ワッチョイ 22c9-+YZw [59.139.198.126])
23/08/06 11:08:31.71 8o4WLPlZ0.net
>>744
言っておくがここはアンタ以外東大卒だからなw
さらに脳内医者もアンタだけw
Rがどうとか言ってるのも脳内医者だけw
782:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 11:10:20.15 sYHzVRqAa.net
358:132人目の素数さん:[sage]:2023/07/28(金) 16:52:04.72 ID:UEEGgPFt
181:卵の名無しさん:[sage]:2023/07/28(金) 16:49:51.11 ID:mTfI1Jyr
>>187
732:卵の名無しさん:[sage]:2021/12/29(水) 19:18:36.08 ID:nTHc+sRn
145 卵の名無しさん[sage] 2021/12/29(水) 18:35:39.75 ID:mlY+9euZ
>>150
「いや、もう火遊びはやめたのよ」
「先生の家はHIでしょ?」
頭のいいナースと仕事できると楽しくて( ・∀・)イイ!!
180:卵の名無しさん:[sage]:2023/07/28(金) 15:23:44.55 ID:UlP+1rs6
>>184
話をしたのが懇意なナースだったので、
オペ台や検査台の上での同意は無効になる可能性があるんだけど、
「**してくれるなら****あげる」というベッドでの同意はどう思う?と聞いたら
「娘によく教えておかなくちゃ」と言っていた。(一部脚色w)
医師板における尿瓶ジジイの看護師に関する気色悪いレスをまとめてみました
783:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 11:15:12.57 sYHzVRqAa.net
>>745
バグだらけなのはアンタのオツムだろ
784:卵の名無しさん (スップー Sd22-fPrT [1.73.23.145])
23/08/06 11:31:03.40 mqqkSGgrd.net
>>745
FDAwww
やっぱお前、医者じゃなくて医薬品とか医療機器扱ってる系の奴だろ
何で医者がFDAに論文出すんだよ
墓穴掘りすぎだろwww
785:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 12:58:26.15 sYHzVRqAa.net
胆汁ドレナージ(商品名を医療用語と勘違い)とか言ってるあたり多分マジだなw
786:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 13:57:48.56 BOo6XOy90.net
'%.%' <- function(x,y){
if(is.complex(x)|is.complex(y)){
a=c(Re(x),Im(x))
b=c(Re(y),Im(y))
}else{
a=x
b=y
}
sum(a*b)
}
ABC2S3d <- function(A,B,C){ cross<-function(x,y)c(x[2]*y[3]-x[3]*y[2],x[3]*y[1]-x[1]*y[3],x[1]*y[2]-x[2]*y[1])
sqrt(sum(cross(A-C,B-C)^2))/2
}
O=c(0,0,0)
A=c(0,0,1)
f=function(x){
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
sum(diff(c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C)))^2) + ((O-A)%.%(B-C))^2 + ((O-B)%.%(C-A))^2 + ((O-C)%.%(A-B))^2
}
optim(runif(6),f)
787:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 13:58:15.60 BOo6XOy90.net
アストラップも商品名。
788:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 13:58:23.70 BOo6XOy90.net
アストラップも商品名。
789:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 14:01:39.42 BOo6XOy90.net
>>750
コロナ関連の論文はRでの解析が多かったぞ。
国立感染症研究所もRのパッケージを使って潜伏期の分布を近似していた。8割おじさんもRとStanで解析。
コードが公開されていたので自分で走らせてみた。
790:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 14:02:19.43 BOo6XOy90.net
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、
哀れだね!
791:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.65.180])
23/08/06 14:13:04.14 t2VeOXVUM.net
これをRで解きたい
四面体OABCは以下の条件を満たす。
(ア)4つの面の面積が全て等しい
(イ)OA⊥BCかつOB⊥CAかつOC⊥AB
この四面体は正四面体か。
792:卵の名無しさん (スップー Sd22-fPrT [1.73.22.213])
23/08/06 14:15:53.37 Q1uPZdHdd.net
>>755
話そらすなよ。何で急にFDAとか出てくるんだ?普通の日本人の医者の発想じゃないんだが。
コロナ関係の論文とか曖昧な表現しないで医者なら�
793:ソゃんとした出典だせや NEJMとかちゃんとしたやつな
794:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 15:59:12.69 BOo6XOy90.net
>>758
FDAはRの解析でも認可するのは常識。
エクセルの統計解析だとだめだけどね。
795:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 16:01:30.33 BOo6XOy90.net
>>758
8割おじさんのスクリプトを走らせたことないの?
ひょっとしたらRを使えないないシリツ医かよ?
796:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.68.57])
23/08/06 16:31:33.80 sYHzVRqAa.net
65:132人目の素数さん:[sage]:2023/08/05(土) 17:16:30.56 ID:AoYbbVsJ
>>61
君たちシダイ卒が求めているようなので貼ってやった。
東大卒による妥当性の検証を希望。
67:132人目の素数さん:2023/08/06(日) 00:26:23.88 ID:5u5eLA7R
ちなみに>>64で歌われてる「数学できない偽医者」ってのが>>65のことね。
こういうイカレた爺さんにならないようにね。>良い子のみなさん
91:132人目の素数さん:[sage]:2023/08/06(日) 13:04:35.89 ID:Nm44ICz6
>>67
残念でした別人です。
俺はシリツと呼んでいるから。
観察力が乏しいね、さてはシリツだな!?
つまり>>65ではなく自分こそが「数学できない偽医者」ですと言いたいようです
797:卵の名無しさん (スップー Sd22-fPrT [1.73.22.213])
23/08/06 17:09:17.76 Q1uPZdHdd.net
>>759
結局、話反らして誤魔化すしか出来ないんだなwww医者なら論文の話をしてるのに普通FDAの話題なんて出さない
お前、どうせ管財課か何かだろ?じゃないと
論文の話してるのに急にFDAの話題なんて出さないよな
FDAって医薬品とか医療機器を仕切ってるアメリカの機関だぞ
Rでオッケーだからどうとかいう問題では無いしそれが理解出来ないなら医者じゃねーよ
医者じゃないから統計ソフトも金のかからんフリーソフト使ってるんだろ?普通はSASとかSPSS使うだろ
798:卵の名無しさん (スップー Sd22-fPrT [1.73.22.213])
23/08/06 17:16:02.67 Q1uPZdHdd.net
>>760
で、結局Rを使用した論文とやたら一切示せないのか?尿瓶ジジイ
で、お前はRを使って論文書いてちゃんとacceptされたのか?脳内医者www
799:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 18:06:24.50 BOo6XOy90.net
O=c(0,0,0)
A=c(0,0,1)
f=function(x){
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
sum(diff(c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C)))^2) + ((O-A)%.%(B-C))^2 + ((O-B)%.%(C-A))^2 + ((O-C)%.%(A-B))^2
}
opt=optim(runif(6),f) ; opt
x=opt$par
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
abs(A-B)
abs(B-C)
abs(C-A)
abs(O-B)
abs(O-C)
c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C))
c((O-A)%.%(B-C), (O-B)%.%(C-A) , (O-C)%.%(A-B))
800:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 18:06:32.07 BOo6XOy90.net
O=c(0,0,0)
A=c(0,0,1)
f=function(x){
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
sum(diff(c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C)))^2) + ((O-A)%.%(B-C))^2 + ((O-B)%.%(C-A))^2 + ((O-C)%.%(A-B))^2
}
opt=optim(runif(6),f) ; opt
x=opt$par
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
abs(A-B)
abs(B-C)
abs(C-A)
abs(O-B)
abs(O-C)
c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C))
c((O-A)%.%(B-C), (O-B)%.%(C-A) , (O-C)%.%(A-B))
801:卵の名無しさん (ワッチョイ 5ec7-lnPD [124.47.76.122])
23/08/06 18:07:30.27 BOo6XOy90.net
>>762
Rの話だが、あんたR使えないシリツかよ?
802:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.148])
23/08/06 19:34:45.47 4URt9vlxa.net
>>766
金なしだから医者が使わないRで発狂してんだろ?
803:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.148])
23/08/06 19:37:09.12 4URt9vlxa.net
>>761続き
98:132人目の素数さん:[sage]:2023/08/06(日) 16:28:05.65 ID:P0hZyEmz
>>93
ということは数学できない偽医者は自分ですってことか?ww
100:132人目の素数さん:[sage]:2023/08/06(日) 18:11:42.36 ID:Nm44ICz6
>>101
さてはシリツだな。
返す言葉もなくて発狂w
804:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.148])
23/08/06 19:38:33.81 4URt9vlxa.net
>>766
R使った論文さっさと出せよ尿瓶クソジジイ
805:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.41])
23/08/06 19:46:54.86 fn2ZpMaBa.net
924:132人目の素数さん:[sage]:2023/08/06(日) 13:09:49.53 ID:zTuxlWyc
>>923
やっぱ医者じゃないんだな
頭悪いのに変だと思ったわ
これが現実だよ尿瓶ジジイ
806:卵の名無しさん (ワッチョイ 5276-Rxlw [223.219.76.206])
23/08/06 21:24:19.99 FizP/rDL0.net
こないだm3を見たら、患者がいきなり出身校を聞いてきて非常に失礼だ!と怒ってる書き込みがあったけど
知りたがる神経質な患者がいるのは当然だし、知られて困るもんでもないので正直に答えてきたんだがやはり底辺私立医か?
807:卵の名無しさん (ワッチョイ a458-fPrT [106.73.2.65])
23/08/06 22:03:21.48 2sD/rCEJ0.net
>>77
808:1 わざわざ出身大学聞かれるってことは、お前の診療が信用ならんってことだろ 普通にやってたら聞かれないよ それか、医者じゃない奴の妄想かどっちかだな
809:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.70.84])
23/08/06 22:45:11.49 0K8PYmyga.net
>>761
>>768
数学できない偽医者さん自分で認めちゃったw
810:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 08:45:02.62 ULIxRH8SM.net
O(0,0,0) A(0,0,100)として
等面積条件は三角形の2辺の外積ベクトルの長さが4つの三角形で等しい、直交条件はベクトルの内積が0として計算。
Σ(外積ベクトルの長さの差)^2 + Σ内積^2 を返す関数を作って
初期値は乱数発生させてL-BFGS-B法にてこの関数の返す値が0になるようなB,Cの座標を探索させて四面体の辺の長さを求める。
これがすべて100であれば正四面体と判定。
811:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 08:46:22.30 ULIxRH8SM.net
重要なclinical pearl
女をみたら妊婦と思え
シリツ医をみたら裏口と思え
812:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.69.93])
23/08/07 09:04:49.93 JgDyrRn4a.net
アホなレスを見かけたら尿瓶ジジイと思え
813:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:11:09.97 ULIxRH8SM.net
>>771
俺も聞かれたので答えたら
「あぁ良かった、自分より馬鹿な医者に診てもらいたくなかったので」と言う返事が返ってきた。
814:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:12:52.15 ULIxRH8SM.net
その話をしたら外来ナースもよく聞かれると言っていた。
シリツ医のときは私達は知らされていませんと答えることにしていると。
815:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:14:48.20 ULIxRH8SM.net
>>771
俺は誇りある母校名を答えられるから歓迎だけどね。
国立大学として新型コロナ診療に後輩が頑張っている姿が報道されていたのも幸い。
816:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:15:40.66 ULIxRH8SM.net
>>769
>774をRに実装してみ!
シリツだとできんだろ。
817:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:17:28.84 ULIxRH8SM.net
>>771
教員、とりわけ校長とか教頭に、偏差値で判断する患者が多いと思う。同僚も同じ印象を持っていた。
818:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:23:10.82 ULIxRH8SM.net
教員だと、あの学力でシリツ医学部合格=裏口とわかっているからだろうな。
819:卵の名無しさん (ブーイモ MM1a-lnPD [49.239.67.161])
23/08/07 11:24:47.07 ULIxRH8SM.net
>>774
問題はこれ
四面体OABCは以下の条件を満たす。
(ア)4つの面の面積が全て等しい
(イ)OA⊥BCかつOB⊥CAかつOC⊥AB
この四面体は正四面体か。
820:卵の名無しさん (ワッチョイ ce12-+YZw [153.137.214.121])
23/08/07 11:39:10.88 7YNAljMk0.net
尿瓶ジジイよりアホな患者なんていないだろ
821:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.69.216])
23/08/07 11:39:58.89 1hB4n964a.net
で、脳内卒業証書と医師免許はまだか?
822:卵の名無しさん (ワッチョイ ce12-+YZw [153.137.214.121])
23/08/07 12:06:49.89 7YNAljMk0.net
>>780
さっさとご自慢のR使った論文出せよ
それとも脳内だからそんなものないのか?
823:卵の名無しさん (ワッチョイ b276-Rxlw [223.218.117.45])
23/08/07 12:57:01.26 BV3idkh10.net
>>785
お前がうpしたら俺も見せるよ
もちろんID付きでうpしろよ
824:卵の名無しさん (スップー Sd22-fPrT [1.73.21.21])
23/08/07 17:52:25.34 x5Kk+yeLd.net
>>787
卒業証明書と医師免許証にIDなんて無いだろwww
持ってませんって言ってるようなもんじゃんwww
825:卵の名無しさん (ワッチョイ ce12-+YZw [153.137.214.121])
23/08/07 17:55:30.85 7YNAljMk0.net
尿瓶ジジイじゃないと思うそいつ
826:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 17:58:09.14 TfKU7SHu0.net
>>769
>774をRに実装してみ!
シリツだとできんだろ。
827:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.197])
23/08/07 17:59:51.15 DXkLQDmpa.net
>>790
で、Rを使った論文まだかよ?
828:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:02:57.87 TfKU7SHu0.net
鯛のアラを使って鯛めしを妻とママごと。
グリルで焼いてから出汁を取ると実に旨くて(・∀・)イイ!!
とうもろこしの青のり塩での味
829:付けも合格点。 まあ、バターコーンの旨さには勝てんが。
830:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:03:48.84 TfKU7SHu0.net
>>791
んであんたR使えないシリツ医だろ?
831:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:04:25.07 TfKU7SHu0.net
重要なclinical pearl
女をみたら妊婦と思え
シリツ医をみたら裏口と思え
832:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:05:21.86 TfKU7SHu0.net
ド底辺シリツ卒でもこんな間違いはしない。
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
833:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:06:58.69 TfKU7SHu0.net
>>791
8割おじさんのペーパーを読んだことないの?
RとStanを使ってたぞ。
834:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.71.197])
23/08/07 18:07:51.00 DXkLQDmpa.net
R使った論文なんかないからそうやって発狂するしか能がないのか
脳内医者はチープだねw
835:卵の名無しさん (ワッチョイ 2458-i2SF [14.13.16.0])
23/08/07 18:09:35.83 TfKU7SHu0.net
>>792
懇意なナースに薦めておこう。
スーパーに鯛のアラがあったらゲットだな。
手間をかけた分だけ旨くなるのが( ・∀・)イイ!!
836:卵の名無しさん (スップー Sd22-fPrT [1.73.21.21])
23/08/07 18:17:14.15 x5Kk+yeLd.net
8割おじさんとかぼろくそに言われてるけどwww
Rとか使ってるから駄目ってことか
読んだことないけど
837:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.69.86])
23/08/07 18:52:32.92 5unrq0o1a.net
お医者さんごっこもあまりにアホすぎると素人にもバレバレなんだねw
838:卵の名無しさん (アウアウウー Sacb-+YZw [106.128.69.86])
23/08/07 18:59:39.76 5unrq0o1a.net
>>796
ソースはネット()
839:卵の名無しさん (ワッチョイ b276-Rxlw [223.218.117.45])
23/08/07 19:33:37.89 BV3idkh10.net
>>788
ここで表示される自分のIDを紙に書いて、現物と一緒に写真に撮るんだよ、本物だと証明するために
840:卵の名無しさん (ワッチョイ a458-fPrT [106.73.2.65])
23/08/07 20:16:30.22 9ukIP9CL0.net
>>802
お前、誰だよ
841:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:23:45.54 9Rchc0KK0.net
'%.%' <- function(x,y){
if(is.complex(x)|is.complex(y)){
a=c(Re(x),Im(x))
b=c(Re(y),Im(y))
}else{
a=x
b=y
}
sum(a*b)
}
ABC2S3d <- function(A,B,C){ cross<-function(x,y)c(x[2]*y[3]-x[3]*y[2],x[3]*y[1]-x[1]*y[3],x[1]*y[2]-x[2]*y[1])
sqrt(sum(cross(A-C,B-C)^2))/2
}
abs = \(z){
if(is.complex(z)){
return(sqrt(Re(z)^2+Im(z)^2))
}else{
return(sqrt(sum(z^2)))
}
}
842:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:24:06.59 9Rchc0KK0.net
a=100
O=c(0,0,0)
A=c(0,0,a)
f=function(x){
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
sum(diff(c(ABC2S3d(O,A,B),ABC2S3d(O,B,C),ABC2S3d(O,C,A),ABC2S3d(A,B,C)))^2) + ((O-A)%.%(B-C))^2 + ((O-B)%.%(C-A))^2 + ((O-C)%.%(A-B))^2
}
opt=optim(runif(6,0,10*a),f, method='L') ; opt
x=opt$par
b1=x[1]
b2=x[2]
b3=x[3]
c1=x[4]
c2=x[5]
c3=x[6]
B=c(b1,b2,b3)
C=c(c1,c2,c3)
abs(A-B)
abs(B-C)
abs(C-A)
abs(O-B)
abs(O-C)
843:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:37:53.39 9Rchc0KK0.net
重要なclinical pearl
女をみたら妊婦と思え
シリツ医をみたら裏口と思え
医師国家試験に出題してもいいくらいだなぁw
844:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:41:12.41 9Rchc0KK0.net
>>777
ちなみにその患者は早稲田卒だった。
早稲田卒に誇りをもっているようであった。
「あぁ良かった」という言葉に
早稲田卒の俺が裏口容疑者に診てもらいたくないという雰囲気を感じたよ。
845:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:42:27.35 9Rchc0KK0.net
真実を語る者が追い詰められる社会だよなぁ。
URLリンク(livedoor.blogimg.jp)
846:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:47:34.76 9Rchc0KK0.net
>>797
こういうのも知らんの?
StanどころかRも使えないんじゃないの?
URLリンク(github.com)
847:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:48:58.17 9Rchc0KK0.net
>>799
あんたにゃ、読める素養すらないだろ。
Stanの使用経験ないと公開コードを走らせることもできんだろうし。
848:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 02:59:02.62 9Rchc0KK0.net
ちょっと調べればでてくるのに
URLリンク(www.okadajp.org)
裏口シリツ医って検索すらできんの?
849:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.2 [上級国民]])
23/08/08 03:01:44.66 9Rchc0KK0.net
>>799
京大教授に栄転したけどね。
RやStanが使える基礎学力がないと読むこともできんだろう。
850:卵の名無しさん (スップ Sd14-fPrT [49.96.233.30])
23/08/08 07:42:16.68 78Ru2Ca7d.net
>>809
アホだなwww
この論文は査読前のacceptされてないやつじゃねーか
論文集めてるサーバーにただ置いてるだけ
医者からしたら何の信頼性もない論文なのに何でドヤ顔で出してるんだよ!
普通の医者なら気づくけどな
やっぱお前医者じゃないだろ
851:卵の名無しさん (スップ Sd14-fPrT [49.96.233.30])
23/08/08 07:53:20.27 78Ru2Ca7d.net
>>811
ちなみにそれもブログレベルの話
852:で、信憑性にかけるだろ そもそも、RをFDAが使ってるとしても論文として出してないんだったら世界レベルで信用されてることにならない だからわざわざ資料って書いてるだろ? やっぱお前どう考えても医者じゃないわ 発想が素人のそれwww
853:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.71.58])
23/08/08 08:47:12.08 o8NGUI8ha.net
>>807
ん?今度は早稲田卒になったのか?
早稲田に医学部なんてあったかなぁ?ww
854:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.71.58])
23/08/08 08:48:14.68 o8NGUI8ha.net
>>809
>>811
さっさとまともな論文もってこいよ尿瓶
855:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.68.204])
23/08/08 16:08:42.68 Wr79/XTba.net
おい尿瓶ジジイ
R使ったまともな論文はどうしたんだよ?さっさと出せよ
856:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.17 [上級国民]])
23/08/08 16:27:39.23 G5rvunoO0.net
ド底辺シリツ卒でもこんな間違いはしない。
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
857:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.17 [上級国民]])
23/08/08 16:27:55.83 G5rvunoO0.net
>774をRに実装することすらできないってド底辺シリツの予感。
ド底辺シリツ未満だな。
ド底辺シリツ卒でもこんな間違いはしない。
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違いを投稿!
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
858:卵の名無しさん (ワッチョイ 0824-swBy [149.50.210.17 [上級国民]])
23/08/08 16:28:54.01 G5rvunoO0.net
バターコーンは高カロリーなので、今日はカレー味にしてみた。
SBの赤缶使用。
859:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.68.21])
23/08/08 16:58:03.59 LcOKde8ta.net
尿瓶またコピペで発狂w
860:卵の名無しさん (JP 0Hc8-fPrT [202.253.111.210])
23/08/08 21:30:35.41 GMfbsjcyH.net
尿瓶ジジイクソ雑魚過ぎてつまんねぇぞ
もう少し頑張れよ・・・
論文ぐらいググれば出てくるだろ!
それすら出来ねぇとかどんだけ能力低いんだよ・・・
料理のことしか書かないなら料理板行けよ
861:卵の名無しさん (スッププ Sdc2-7AGF [49.105.74.104])
23/08/09 08:23:48.91 pY8Seg/xd.net
2014年に福岡県鞍手町で谷尾幸生さん夫妻を殺害した
「ネトゲ実況板のマナ公」が若い女性にストーカー行為を働いてます
以下に該当する女性をご存知の方は危険を伝えるようお願いします
・現在の年齢が28歳前後、名前の読みが「りか」
・約3年前に福岡市周辺の障害者デイケアに参加または勤務していた
・そこで20歳年上の精神障害者の男性に言い寄られてLINEブロックした
0930 既にその名前は使われています 2023/08/08(火) 11:15:17.49
マナちゃんは一生処女だろうね
りかちゃんは俺が奪うけど、これは宿命
逢えるのはいつになるのだろうか…
ID:Sr19Far8(24/25)
マナ公のTikTok(殺人の前歴がありながら堂々と顔出し)
kazra2020
マナ公のTwitter
メドローア最強@DQW
@medlorerika
862:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.69.60])
23/08/09 09:01:37.29 jb61gbTOa.net
>>822
また認知症がひどくなったみたいだね尿瓶ジジイさん
863:卵の名無しさん (ワッチョイ 2458-swBy [14.13.16.0 [上級国民]])
23/08/09 18:11:37.92 IRzyrZyU0.net
ik2ca = \(ai,ak,bi,bk){ # A:緯度ai,経度ak B:緯度bi,経度bk で弧ABの中心角を返す
a=c(cos(ai)*cos(ak),cos(ai)*sin(ak),sin(ai))
b=c(cos(bi)*cos(bk),cos(bi)*sin(bk),sin(bi))
sum(a*b) |> acos()
}
on.sphere =\(n){ # 球面上から一様分布でn個選んだ地点の緯度・経度を返す
ans=NULL
z=runif(n,-1,1)
t=runif(n,-pi,pi)
x=sqrt(1-z^2)*cos(t
864:) y=sqrt(1-z^2)*sin(t) i=asin(z) k=sign(y)*acos(x/sqrt(x^2+y^2)) cbind(i=i,k=k) } f=\(){ x3=on.sphere(3) all(combn(3,2,\(x) ik2ca(x3[x[1],1],x3[x[1],2],x3[x[2],1],x3[x[2],2])) < 1/2) } y=replicate(1e6,f()) mean(y)
865:卵の名無しさん (ワッチョイ dcb1-F7oX [126.169.105.68])
23/08/09 20:23:43.27 G8EIyKgG0.net
医学部医学科の格付け
SSS 東京大学(理科三類)
SS 京都大学
S 大阪大学 慶応義塾大学
--------------- 東大理一の壁 ---------------
A+ 東北大学 東京医科歯科大学 千葉大学 名古屋大学 京都府立医科大学 九州大学
A- 北海道大学 東京大学(理科二類) 横浜市立大学 大阪市立大学 神戸大学 岡山大学
B+ 筑波大学 金沢大学 広島大学 熊本大学 防衛医科大学校
B- 新潟大学 名古屋市立大学 奈良県立医科大学 長崎大学 自治医科大学 東京慈恵会医科大学
--------------- 京大上位学部(薬学部・理学部)の壁 ---------------
C+ 群馬大学 信州大学 岐阜大学 三重大学 和歌山県立医科大学 鹿児島大学 順天堂大学
C 弘前大学 鳥取大学 山口大学 徳島大学 日本医科大学 大阪医科大学
C- 札幌医科大学 福島県立医科大学 山梨大学 浜松医科大学 滋賀医科大学 産業医科大学 関西医科大学
--------------- 京大工学部の壁 ---------------
D+ 山形大学 富山大学 福井大学 香川大学 愛媛大学 大分大学 昭和大学 東京医科大学
D- 旭川医科大学 秋田大学 島根大学 高知大学 佐賀大学 宮崎大学 琉球大学 日本大学 東邦大学
--------------- 京大農学部、東工大の壁 ---------------
E+ 岩手医科大学 東北医科薬科大学 東京女子医科大学 国際医療福祉大学 近畿大学 福岡大学 久留米大学
E- 独協医科大学 埼玉医科大学 北里大学 杏林大学
--------------- 早慶理工の壁 ---------------
F 帝京大学 東海大学 聖マリアンナ医科大学 愛知医科大学 藤田保健衛生大学 金沢医科大学 兵庫医科大学 川崎医科大学
--------------- 上智理工、理科大、同志社理工の壁 ------
866:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.69.173])
23/08/10 07:27:55.47 mV+LJmOYa.net
尿瓶ジジイまた数学板で発狂中
867:卵の名無しさん (JP 0Hf8-swBy [217.138.252.98 [上級国民]])
23/08/10 07:32:06.82 Mskly8ggH.net
直径10cm 高さ5cmの中華饅頭(球を平面で切断した形、球欠 spherical segmentとして計算)を中心から3cm 離れたところで包丁で切ったときの断面積を計算せよ。
答は有効数字3桁でよい。
参考画像 URLリンク(img1.kakaku.k-img.com)
868:卵の名無しさん (アウアウウー Sacb-IPGJ [106.128.69.173])
23/08/10 20:33:56.33 mV+LJmOYa.net
チンパンジーは生ゴミでも食ってろ
869:卵の名無しさん (JP 0Hf8-swBy [217.138.212.122 [上級国民]])
23/08/11 06:15:28.23 WFgOkZ+xH.net
ik2xyz=\(i,k){
x=cos(i)*cos(k)
y=cos(i)*sin(k)
z=sin(i)
c(x,y,z)
}
flg=FALSE
while(!flg){
ABCik=cbind(runif(3,-pi/2,pi/2),runif(3,-pi,pi))
ABC=apply(ABCik,1,\(ik) ik2xyz(ik[1],ik[2]))
P=ABC[1,]/2+ABC[2,]/2
Q=ABC[2,]/2+ABC[3,]/2
R=ABC[3,]/2+ABC[1,]/2
flg = sum(P^2)<1/4 & sum(Q^2)<1/4 & sum(R^2)<1/4
}
870:卵の名無しさん (JP 0Hf8-swBy [217.138.212.122 [上級国民]])
23/08/11 07:40:50.24 WFgOkZ+xH.net
"
半径1の球の表面および内部から無作為に3点A B Cを選ぶ。
球を構成するどの点も選ばれる確率は等しいとする。
三角形ABCの面積をSのとするとき、
(1) Sが1/2より大きい確率を求めよ(有効数字3桁でよい)
(2) Sの期待値を求めよ(有効数字3桁でよい)
(3) Sの分布を図示せよ
・
"
source('toolmini.R')
sim1=\(){
ABC=in.sphere(3)
ABC2S3d(ABC[1,],ABC[2,],ABC[3,])
}
S=replicate(1e6,sim1())
summary(S)
BEST::plotPost(S,compVal = 1/2)
mean(S>1/2)
871:卵の名無しさん (JP 0Hf8-swBy [217.138.212.122 [上級国民]])
23/08/11 17:02:11.89 WFgOkZ+xH.net
どの面の面積も等しく、ある1つの球に4頂点が内接する四面体は正四面体か。
Do=\(par,rotate=FALSE){
opt=optim(c(par),f,method = 'L')
print(opt)
ABC=opt$par
Ai=ABC[1]
Ak=ABC[2]
Bi=ABC[3]
Bk=ABC[4]
Ci=ABC[5]
Ck=ABC[6]
D=c(0,0,1)
A=ik2xyz(Ai,Ak)
B=ik2xyz(Bi,Bk)
C=ik2xyz(Ci,Ck)
c(ABC2S3d(D,A,B),ABC2S3d(D,B,C),ABC2S3d(D,C,A),ABC2S3d(A,B,C)) |> print()
c(abs(A-B),abs(B-C),abs(C-A),abs(D-A),abs(D-B),abs(D-C)) |> print()
ABCD=rbind(A,B,C,D)
round(ABCD,6)
plot3d(ABCD,col=2,size=10,xlab='x',ylab='y',zlab='z') ; aspect3d('iso')
combn(4,2,\(x) segments3d(ABCD[c(x[1],x[2]),1], ABCD[c(x[1],x[2]),2],ABCD[c(x[1],x[2]),3],col=2))
if(rotate) Play3d(15)
invisible(opt$par)
}
872:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/12 18:43:55.37 hE0YxBJe0.net
addon2int=\(I,N){
A=1e8
Total=A+A*I*N/12
int1=ganri(A=A,N=N,x=Total/N)
int2=uniroot(\(i) gankin(A=A,i,N=N)[2] - Total,c(0,1)$root
c(int1,int2)
}
873:卵の名無しさん (JP 0Hb2-FnP7 [217.138.212.122 [上級国民]])
23/08/12 19:00:14.86 NTQGy9PqH.net
addon2int=\(I,N){
A=1e8
Total=A+A*I*N/12
int1=ganri(A=A,N=N,x=Total/N)
# names(int1)=NULL
int2=uniroot(\(i) gankin(A=A,i,N=N)[[2]] - Total,c(0,1),tol = 1e-16)$root
data.frame(元利均等返済=int1,元金均等返済=int2)
}
addon2int(0.05,120)
874:卵の名無しさん (JP 0Hb2-FnP7 [217.138.212.122 [上級国民]])
23/08/12 19:00:29.12 NTQGy9PqH.net
> addon2int(0.05,120)
元利均等返済 元金均等返済
年利 0.08689224 0.09917355
875:卵の名無しさん (JP 0H03-FnP7 [146.70.31.52 [上級国民]])
23/08/13 06:23:40.12 30EpRvIzH.net
長さが2,3,4,5,6,7の6本の線分を組み合わせて最も鋭利な頂点をもつ三角錐を作る。
底面を三角形ABC、頂点をDとし頂点Dが最も鋭利とする。
(1)∠ADB+∠BDC+∠CDAは何度か?
(2)その三角錐の高さを求めよ
(3)その三角錐の体積を求めよ。
(4)その三角錐を図示せよ。
> fn(pm[imin[1],]) ; aspect3d('iso')
$D
[1] -2.500000 -2.452889 3.568380
$Vol
[1] 3.455069
$S
[1] 26.06022
$s
[1] 2.904738
$h
[1] 3.56838
$apex
[1] 74.75361
876:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.23.123])
23/08/13 10:22:31.24 hCmtcqpIa.net
尿瓶ジジイ数学板でももう相手にされてないみたい
日本中を騒がせた某教祖みたいにブツブツ何か言ってるだけ
877:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/13 14:01:44.54 e8Wi33VE0.net
>>837
東大卒は計算してレスしている。
東大合格者でないと答が出せないみたい。
878:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.32.188])
23/08/13 14:2
879:0:50.07 ID:Pticxi9Ga.net
880:卵の名無しさん (RU 0H3b-FnP7 [80.66.69.29 [上級国民]])
23/08/13 16:23:51.80 S2lyrM4lH.net
>>839
進振りの話とか体験者ならではのレスをしていたからね。
んで、あんたの自称は?
どこの国立を落ちたの?
881:卵の名無しさん (RU 0H3b-FnP7 [80.66.69.18 [上級国民]])
23/08/13 16:36:19.86 qn/UONEYH.net
v=rep(1:13,4)
cm=comboGeneral(v,5)
full_house=\(x){
if(length(unique(x))!=2) return(FALSE)
else if(!(4 %in% table(x))) return(TRUE)
return(FALSE)
}
apply(cm,1,full_house) |> sum() -> fh # 3744/2598960
fh/nrow(cm)
fh/nrow(cm) |> MASS::fractions()
"
ジョーカー1枚を含めて53枚のトランプから無作為に5枚取り出したとき
フルハウスである確率を求めよ。
尚、
1,1,2,1,Jokerはフルハウス
1,1,1,2,Jokerはフォーカードとみなしてフルハウスにカウントしない。
"
x=c(0,1,1,2,2) # full house
x=c(0,1,1,1,3) # 4 card
x=c(0,1,1,2,3) # 3 card
x=c(0,1,2,3,4) # one pair
vj=c(0,v)
cmj=comboGeneral(vj,5)
cmj[987654,]
tail(cmj)
full_houseJ =\(x){
if(0 %in% x){
return(all(table(x[x!=0])==2))
}
else return(full_house(x))
}
replicate(1e6,full_houseJ(cmj[sample(nrow(cmj),1),])) |> mean()
882:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.32.188])
23/08/13 17:18:07.36 Pticxi9Ga.net
>>840
東大の進振りなんて証明にならんだろw
アホすぎるから尿瓶ジジイの自称は誰にも相手にされてないから自分が建てたクソスレでしか発狂できないんだろうがw
883:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.32.188])
23/08/13 17:27:15.02 Pticxi9Ga.net
49:132人目の素数さん:[sage]:2023/08/13(日) 14:05:02.12 ID:cES63u1X
マネーリタラシーのある人のための練習問題
「利息制限法の上限金利15%で月々25万ずつ25年ローンで返済するとして借りられる元金はなんぼや?」
どうも尿瓶チンパポンコツフェチのキーキー電卓は金利計算対応していないようだ。
こんなマヌケなレスするやつのどこが医者なん?w
リタラシーww
884:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/13 18:35:10.06 DUnW1R9sH.net
尿瓶おまる洗浄係と違って臨床やっていると実体験に基づいてこういう投稿ができる。
URLリンク(egg.2ch.sc)
支持するレスがついたな。
885:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.32.188])
23/08/13 18:56:20.14 Pticxi9Ga.net
どこが支持だよ脳内医者
886:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 07:16:59.03 lOjCl5CzH.net
flash=\(x) x%/%13 |> unique() |> length() == 1
str8=\(x){
y=x%%13
y[y==0]=13
y |> sort() |> diff() |> unique() |> length() == 1 -> re1
y[y==1]=14
y |> sort() |> diff() |> unique() |> length() == 1 -> re2
re1 | re2
}
887:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 07:21:07.66 lOjCl5CzH.net
>>845
困ったときの右下側臥位。
CS経験ない椰子は知らんだろ、尿瓶おまる洗浄係くん
888:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 07:27:41.79 lOjCl5CzH.net
flashJ=\(x){
x=x[x!=0]
x%/%13 |> unique() |> length() == 1
}
str8=\(x){
y=x%%13
y[y==0]=13
y |> sort() |> diff() |> unique() |> length() == 1 -> re1
y[y==1]=14
y |> sort() |> diff() |> unique() |> length() == 1 -> re2
re1 | re2
}
889:卵の名無しさん (ブーイモ MMee-20OZ [133.159.152.55])
23/08/14 08:41:02.76 gtktrR2OM.net
flashJ=\(x){
x=x[x!=0]
x%/%13 |> unique() |> length() == 1
}
str8J=\(x){
if(0 %in% x){
x=x[x!=0]
y=x%%13
y[y==0]=13
y |> sort() |> diff() -> d1
d1=c(1,1,1) | d1==c(1,1,2) | d1==c(1,2,1) | d1=c(2,1,1) -> re1
y[y==1]=14
y |> sort() |> diff() -> d2
d2=c(1,1,1) | d2==c(1,1,2) | d2==c(1,2,1) | d2=c(2,1,1) -> re2
}else{
y=x%%13
y[y==0]=13
y |> sort() |> diff() |> unique() |> length() == 1 -> re1
y[y==1]=14
y |> sort() |> diff() |> unique() |> length() == 1 -> re2
}
return(re1 | re2)
}
890:卵の名無しさん (ブーイモ MMee-20OZ [133.159.152.55])
23/08/14 11:50:11.82 gtktrR2OM.net
flashJ=\(x){
x=x[x!=0]
x%/%13 |> unique() |> length() == 1
}
f1=\(d){
all(d==c(1,1,1))|all(d==c(1,1,2))|all(d==c(1,
891:2,1))|all(d==c(2,1,1)) } f2=\(y) y |> sort() |> diff() |> unique() |> length() == 1 str8J=\(x){ if(0 %in% x){ x=x[x!=0] y=x%%13 y[y==0]=13 y |> sort() |> diff() |> f1() -> re1 y[y==1]=14 y |> sort() |> diff() |> f2() -> re2 }else{ y=x%%13 y[y==0]=13 y |>f2() -> re1 y[y==1]=14 y |>f2() -> re2 } return(re1 | re2) } str8J(c(1,2,3,0,5))
892:卵の名無しさん (JP 0Hfe-eBN6 [61.199.225.248])
23/08/14 12:55:40.97 8/MmQB6PH.net
684:卵の名無しさん (スッププ Sda2-Xleb [49.105.85.56]):[sage]:2023/08/03(木) 17:32:32.22 ID:69uB5F4Yd
尿瓶ジジイも自分の分をわきまえたのか
間違った脳内医療をここに晒すことはなくなったなwww偉いじゃん
尿瓶ジジイは内視鏡スレには現れるのに、消化器専門医スレには現れない
何故なんでしょうねwww
ヒント:内視鏡はコメディカルでもある程度理解できます!
893:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.27.76])
23/08/14 12:57:46.36 iX2BHMk9a.net
あとそれぞれ専門科があるのにいちいち経験でマウントとろうとするのも愚の骨頂なw
脳内医者は何でもできる設定みたいだが
894:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.27.76])
23/08/14 13:01:59.58 iX2BHMk9a.net
尿瓶リタラシージジイはいつまで恥を晒し続けるのか
895:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 13:38:32.43 lOjCl5CzH.net
> flash=apply(cm,1,flashJ)
> mean(flash)
[1] 0.002541394
> str8=apply(cm,1,str8J)
> mean(str8)
[1] 0.02493828
> str8flash= apply(cm[str8,],1,flashJ)
> mean(str8flash)
[1] 0.009348145
896:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/14 14:12:38.54 0qVMbDv90.net
>>852
Rは専門を問わず使えた方がいい。
897:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/14 14:13:36.77 0qVMbDv90.net
専門が何であれ
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
>当たり前だろ
>当たり前だろ
>当たり前だろ
とドヤ顔で大間違い!
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
898:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/14 14:15:05.22 0qVMbDv90.net
ジョーカーを含めた53枚のトランプから無作為に5枚を選ぶ
(1)フラッシュの確率を求めよ。
(2)ストレートの確率を求めよ。
(3)ストレートフラッシュの確率を求めよ。
899:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 14:48:03.14 lOjCl5CzH.net
>>852
一般外科やってたら、開腹・開胸手術、アンギオ、内視鏡、麻酔くらいできるだろ。
シャント造設から大腿骨頸部骨折手術までできる一般外科医と仕事をしたことがある。
frozen pelvis症例には自分で膀胱鏡下にステントを留置していたなぁ。
んで、あんたは薬屋の売り子以外に何ができるの?
900:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/14 14:51:35.20 lOjCl5CzH.net
マウントをとるとかいう発想はサルだね。
ボノボまで進化するとマウントをとるとかしない。
901:卵の名無しさん (スップ Sda2-iQuL [49.97.8.29])
23/08/14 18:05:03.10 gSPWucFjd.net
>>858
開胸とかアンギオとか麻酔とか、えらいアバウトだな
開胸っていったら心臓血管外科の手術もあるのに本当に出来るんか?
アンギオって具体的にどこをどう造影すんだよ。全部出来るって意味?
麻酔もヘルツのオペとかもあんのに出来るってことか?
ちなみに少し触ったことあるのは出来るって言わないんだよ
責任持ってちゃんと出来ることをできるって表現するんだけど、本当に責任持って仕事してる奴はお前みたいに何でもできるって言わないんだわ
お前の文章は具体性が皆無で素人の脳内妄想にしか見えんわ
902:卵の名無しさん (ワッチョイ bb58-iQuL [106.73.2.65])
23/08/14 22:40:21.94 463kk6ES0.net
209 卵の名無しさん sage 2023/08/14(月) 14:11:14.64 ID:0qVMbDv9
>>211
でもボルンホルム病の記載はなかったな。
ボルンホルム島ってどこにあるのかすら知らなかったがノルドストリーム破壊事件でデンマーク領と知った。
ACSのシマウマ鑑別にあげるとマニアックでいい。
実際に経験したことないが。
↑
これも、尿瓶の投稿だろうけど
内容が支離滅裂過ぎて笑えるわ
言葉の意味とか理解しないでそれらしい単語羅列してるんだろうな
統合失調症発症し過ぎだろ
903:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.21.253])
23/08/14 22:41:23.47 LozjRxFKa.net
言葉のサラダって知ってるか?尿瓶ジジイ
904:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/15 04:51:33.94 JPJRanbe0.net
>>860
んで、あんたは薬屋の売り子以外に何ができんの?
905:卵の名無しさん (ワッチョイ 7f58-e563 [14.13.16.0 [上級国民]])
23/08/15 04:52:42.52 JPJRanbe0.net
>悪性症候群で抗ドーパミン薬止めない
というのができるのは知っている。
とても俺にはできんけどねw
906:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.21.253])
23/08/15 07:05:13.53 LrgH0bqGa.net
尿瓶ジジイ支離滅裂に関しては反論できなくて発狂かよ
907:卵の名無しさん (JP 0H03-FnP7 [146.70.76.185 [上級国民]])
23/08/15 07:11:55.54 rtnMmaAJH.net
>>850 バグ修正
flashJ=\(x){
x=x[x!=0]
(x-0.1)%/%13 |> unique() |> length() == 1
}
f1=\(d){
all(d==c(1,1,1))|all(d==c(1,1,2))|all(d==c(1,2,1))|all(d==c(2,1,1))
}
f2=\(y) y |> sort() |> diff() |> unique() |> length() == 1
str8J=\(x){
if(0 %in% x){
x=x[x!=0]
y=x%%13
y[y==0]=13
y |> sort() |> diff() |> f1() -> re1
y[y==1]=14
y |> sort() |> diff() |> f2() -> re2
}else{
y=x%%13
y[y==0]=13
y |>f2() -> re1
y[y==1]=14
y |>f2() -> re2
}
return(re1 | re2)
}
908:卵の名無しさん (スップ Sda2-iQuL [49.97.9.108])
23/08/15 07:44:41.07 m1xx/2Xtd.net
>>863
結局、具体性に何出来るのか言えないわけね
お前が、素人の癖に適当に発言してることが良くわかったよ
909:卵の名無しさん (ブーイモ MMee-20OZ [133.159.148.208])
23/08/15 08:11:18.81 7fgfJRfYM.net
>>867
アッペ、ヘルニア、ヘモ、胆摘、胃切、腸切、肝切、PD、食道切除、骨盤内臓全摘、乳切、肺切とかやったよ。
ENBDやPTCDで緊急減黄やTAEでの止血も外科医の仕事だった。
緊急帝王切開の麻酔や助手も外科医の仕事。
紫色の新生児に挿管したのは思い出しても冷汗もの。
んで、あんたは
>悪性症候群で抗ドーパミン薬止めない
というのができるののよね!
とても俺にはできないなぁw
910:卵の名無しさん (ブーイモ MMee-20OZ [133.159.148.208])
23/08/15 08:15:37.48 7fgfJRfYM.net
心臓麻酔も数例やらされた。
まあ、ポンプマンがしっかりしていたので
挿管係みたいなものだった。
シャントやグラフト取りも指導医のもとに術者をさせてもらった。
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってwww
とかいう医師はシリツ医にもいなかった。
typoだと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。
911:卵の名無しさん (ブーイモ MMee-20OZ [133.159.148.208])
23/08/15 08:28:58.80 7fgfJRfYM.net
>>861
救急外来で診断がつかない胸痛の原因として知る人ぞ知る疾患。
これは俺も知らなかったw、当たり前ですって。
>悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
912:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.21.253])
23/08/15 08:39:15.66 LrgH0bqGa.net
用語羅列すれば脳内ブラックジャックできると思ってんのか?w
913:卵の名無しさん (アウアウキー Sab7-eBN6 [182.251.21.253])
23/08/15 08:41:24.97 LrgH0bqGa.net
尿瓶ジジイは医者はなんでもできると思ってる素人みたい
どこの専門か聞いてんだよw