23/07/24 11:58:38.90 Day8w5dSM.net
救急当番と内視鏡バイトで２０諭吉ゲット！
コロナの救急搬送が立て続けで忙しかった。
625:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.5])
23/07/24 12:01:49.75 Day8w5dSM.net
これより恥ずかしい投稿を医師板ではみたことがない。
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違い！！
626:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.69.134])
23/07/24 12:17:24.25 NlJK5bBea.net
>>602
胆汁ドレナージなんて言ってるほうがよっぽど恥ずかしいよ
脳内だから自覚ないと思うけどw
627:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.118])
23/07/24 12:41:20.13 7Fiupx7ba.net
尿瓶ジジイ、突っ込まれるのが怖いからか前より脳内医療ネタ()は短文w
628:卵の名無しさん (JP 0Hfd-leiM [202.253.111.210])
23/07/24 15:37:34.55 WkpJ9wPeH.net
>>600
胆汁ドレナージは普通には使わねぇよ
さらっと嘘つくなwww
本当にゴミだな
お前、周りから粘着質でキモいって思われてるぞ間違いなく
629:卵の名無しさん (ﾜｯﾁｮｲ c958-XbVB [14.13.16.0 [上級国民]])
23/07/24 16:30:38.18 8cl3Qnmw0.net
>>605
使うよ。
AstrapもMDLも普通に使う。
Homan鈎のフランス後読みは上級者のみ使う。
こういうことは下衆でも言わない。
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
630:卵の名無しさん (JP 0Hd3-zEeY [217.138.212.122 [上級国民]])
23/07/24 16:37:18.50 yZNTF7aqH.net
>>588
これでdat落ちした投稿が読めた。
ｽﾚﾘﾝｸ(math板:157番)
631:卵の名無しさん (ﾜｯﾁｮｲ c958-XbVB [14.13.16.0 [上級国民]])
23/07/24 16:43:47.20 8cl3Qnmw0.net
昔の投稿も見ることができて便利。神アプリだな。
11 １３２人目の素数さん sage 2020/03/10(火) 12:10:41.22 ID:H1fx2jVB
"
SEIR MODEL
dS(t)/dt = mu*(N-S) - b*S(t)*I(t)/N - nu*S(t)
dE(t)/dt = b*S(t)I(t)/N - (mu+sig)*E(t)
dI(t)/dt = sig*E(t) - (mu+g)*I(t)
dR(t)/dt = g*I(t) - mu*R + nu*S(t)
mu
632::自然死亡率 b:感染率(S->I) nu:ワクチン有効率(S->R) sig:発症率(E->I),g:回復率(I->R) " SEIRモデルのパラメータ SEIR2 <- function( # Parameters contact_rate = 10, # number of contacts per day transmission_probability = 0.01, # transmission probability beta = contact_rate * transmission_probability, # tranmission rate infectious_period = 20, # infectious period gamma = 1 / infectious_period, # Prob[infected -> recovered] latent_period = 5, # latent perior sigma = 1/latent_period, # The rate at which an exposed person becomes infective mu = 0, # The natural mortality rate nu = 0 , # vaccination moves people from susceptible to resistant directly, without becoming exposed or infected. Ro = beta/gamma, # Ro - Reproductive number. # Initial values for sub-populations. s = 99, # susceptible hosts e = 0, # exposed hosts i = 1, # infectious hosts r = 0, # recovered hosts # Compute total population. N = s + i + r + e, # Output timepoints. timepoints = seq (0, 365, by=0.5), ... ) 有病率を１％とすると、3000人にクルーズ船でも100人の屋形船でも感染者のピークは変わらないな。 同一時間あたりのcontact_rateとtransmission_probabilityが宴会での方が高いからだろうな。 パラメータを変えてグラフを書いてみた。 https://i.imgur.com/hCfBTyc.png
633:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.68.244])
23/07/24 16:45:15.28 K7qj9qqDa.net
>>606
認めたくないだろうけど使いませんw
胆汁ドレナージなんて言ってる時点で脳内医者確定なんだよなw
634:卵の名無しさん (ﾜｯﾁｮｲ 0b12-TCAR [153.137.214.121])
23/07/24 16:47:08.63 tTpqXxCT0.net
>>606
じゃあなんで教科書でもGoogleでも胆汁ドレナージが出てこないんだよタコ
ソース出せよ
635:卵の名無しさん (JP 0Hd3-zEeY [217.138.212.122 [上級国民]])
23/07/24 18:31:08.42 yZNTF7aqH.net
>>597
抗ドーパミン薬の意味はわかってた？？
636:卵の名無しさん (JP 0Hd3-zEeY [217.138.212.122 [上級国民]])
23/07/24 18:33:50.73 yZNTF7aqH.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
637:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.68.244])
23/07/24 19:12:08.92 K7qj9qqDa.net
抗ドーパミン薬の意味が分かってないのはアンタだろw
638:卵の名無しさん (ｽﾌﾟｰﾌﾟ Sd73-leiM [1.73.153.125])
23/07/24 19:20:37.73 RG128V21d.net
>>611
お前は認めたくないんだろうが>>334で、すでに論破している
これを反論出来ずに永遠とコピペを繰り返すのは負けてることを認めたくないから、レスの物量で誤魔化そうとしてるか
もしくは、どういう意味か理解出来てないからだろwww
個人的には統合失調症で、考えがまとめられないからだと思うが
見えない敵が見えてるのか俺がスレ覗いていない時も永遠とコピペ続けてたよな？やっぱお前は、統合失調症だよ。医者の俺が言うんだから間違いないよ
639:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.68.244])
23/07/25 08:37:24.35 R29t5DOLa.net
尿瓶ジジイ高校生相手には吠えまくってるのにここではダンマリかよ
640:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 08:40:57.23 4Uu2jZLX0.net
"
楕円 x^2/a^2+y^2/b^2=1が通る格子点の座標を返してグラフ化する。
"
rm(list=ls())
source('toolmini.R')
'%=%' <- function(x,y) abs(x-y) <
641: .Machine$double.eps^0.5 ellatice =\(a,b,print=FALSE){ ecl=\(x,y) (x/a)^2 + (y/b)^2 A=ceiling(a) B=ceiling(b) g=as.matrix(expand.grid(-A:A,-B:B)) i=which(apply(g,1,\(z) ecl(z[1],z[2])) %=% 1) ans=g[i,] rownames(ans)=NULL colnames(ans)=c('x','y') if(print){ Plot(-A,A,-B,B) points(g,cex=0.5) Ellip(0i,a,b,col=8) points(ans,col='red',pch=16) } return(list(ans,nrow(ans))) } f=\(a,b){ ans=ellatice(a,b) ans[[1]] |> as.vector() |> prod() !=0 } f=Vectorize(f) outer(1:1000,1000,f) |> sum() ellatice(sqrt(360),sqrt(24),print=TRUE)
642:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 08:42:07.06 4Uu2jZLX0.net
>>609
小児外科医が使っていたと同業者が投稿していたな。
こういうことは下衆でも言わない。
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
643:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 08:43:17.53 4Uu2jZLX0.net
>>614
わざとtypo?恥ずかしい！！！
>>369
ちなみに、抗ドーパミン薬って書いたのも
わざとだよw
少し付け入る隙を与えてやらんとお前の頭悪すぎて、急に黙り込むからな
それでも、気づくのに5日かかったわけだがw
他に何か脳内医療自慢出来ないの？お前を馬鹿にするのがこのスレの趣旨なんだから
無い頭ふり絞ってもう少し頑張れよwww
644:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 08:43:56.77 4Uu2jZLX0.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
645:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 08:44:23.02 4Uu2jZLX0.net
国試だと禁忌肢になるだろな。
646:卵の名無しさん (JP 0Hd3-yMPU [217.138.212.98 [上級国民]])
23/07/25 09:46:05.76 C7I6Uch3H.net
'%=%' <- function(x,y) abs(x-y) < .Machine$double.eps^0.5
ellatice =\(a,b,dx,dy,print=FALSE){
ecl=\(x,y,dx,dy) ((x+dx)/a)^2 + ((y+dy)/b)^2
A=ceiling(a)+1
B=ceiling(b)+1
g=as.matrix(expand.grid(-A:A,-B:B))
i=which(apply(g,1,\(z) ecl(z[1],z[2],dx,dy)) %=% 1)
ans=g[i,]
rownames(ans)=NULL
colnames(ans)=c('x','y')
if(print){
Plot(-A,A,-B,B)
points(g,cex=0.5)
ellip(dx,dy,a,b,col=8)
points(ans,col='red',pch=16)
}
return(list(ans,nrow(ans)))
}
ellatice(2,4,0,0,print=T)
ellatice(sqrt(360),sqrt(24),0,0,print=T)
flowerfrog <- \(a,b){
flg <- ellatice(a,b,0,0)[[2]]!=0
if(flg) ellatice(a,b,0,0,print=TRUE)
while(flg){
dx=runif(1,-1,1)
dy=runif(1,-1,1)
flg <- ellatice(a,b,dx,dy)[[2]]!=0
}
ellip(dx,dy,a,b,col='blue',lwd=2)
return(c(dx,dy))
}
flowerfrog(sqrt(360),sqrt(24))
flowerfrog(5,4)
647:卵の名無しさん (JP 0Hfd-leiM [202.253.111.210])
23/07/25 12:57:37.42 +LRtxEUHH.net
>>618
コメントが「恥ずかしい！」だけかよ
結局、反論出来ないんだろwww
これ以上何か言っても馬鹿晒すだけだぞ
本当に頭悪いんだな
偏差値30ぐらいと奴でも、もう少しちゃんとした反論出来るぞ
それより、サンドラマランって何なの？ちゃんと解説してくれよwww東大合格者なんだろ？
648:卵の名無しさん (ﾜｯﾁｮｲ 0b12-TCAR [153.137.214.121])
23/07/25 13:43:36.25 3W/aXRhE0.net
ちなみに高校生からもボコられてますw
649:卵の名無しさん (ﾜｯﾁｮｲ f324-h0SB [149.50.210.2 [上級国民]])
23/07/25 17:09:32.66 4Uu2jZLX0.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
typoだと言ったり、
わざと間違えたと言ったり、
恥ずかしくね？
650:卵の名無しさん (ｱｳｱｳｷｰ Saf5-TCAR [182.251.24.65])
23/07/25 17:12:44.66 vyAkGQ6Ca.net
高校生にも偽医者とか言われてる尿瓶ジジイのほうが恥ずかしいから安心してw
651:卵の名無しさん (ｽﾌﾟｰﾌﾟ Sd73-leiM [1.73.155.226])
23/07/25 18:31:34.66 8MuMsPUSd.net
>>624
恥ずかしいわけないだろ、お前を馬鹿にするためにわざとやってんだからwww
これだけ、言っても理解出来ないの？
まあ、頭悪くて他に反論出来ないからこんな書き込みが出来るんだろうな
コピペしか出来ないお前の方がよっぽど恥ずかしいわ！www
んで、結局サンドラマランは無視か？
Syndrome malinをサンドラマランとかドヤ顔で
652:書いたら俺なら恥ずかしくて自殺しちゃうかもなwww
653:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 06:31:40.62 0pDrMthqH.net
rm(list=ls())
source('toolmini.R')
library(RcppAlgos)
library(rgl)
ABC=abc2ABC(5,6,4)
A=ABC[1]
B=ABC[2]
C=ABC[3]
Plot(0,6)
pta(A)
pta(B)
pta(C)
Polygon(A,B,C)
A=c(Re(A),Im(A),0)
B=c(Re(B),Im(B),0)
C=c(Re(C),Im(C),0)
pm=permuteGeneral(c(7,8,9),3)
f=\(D,X) rbind(D,X) |> dist()
for(i in 1:6){
fn=\(D) (f(D,A)-pm[i,1])^2 + (f(D,B)-pm[i,2])^2 + (f(D,C)-pm[i,3])^2
optim(c(1,4,7),fn,method='L')$par |> print()
}
i=1
fn=\(D) (f(D,A)-pm[i,1])^2 + (f(D,B)-pm[i,2])^2 + (f(D,C)-pm[i,3])^2
optim(c(1,4,7),fn,method='L')$par |> print()
654:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 07:20:44.08 0pDrMthqH.net
source('toolmini.R')
library(RcppAlgos)
library(rgl)
ABC=abc2ABC(5,6,4)
A=ABC[1]
B=ABC[2]
C=ABC[3]
Plot(0,6)
pta(A)
pta(B)
pta(C)
Polygon(A,B,C)
A=c(Re(A),Im(A),0)
B=c(Re(B),Im(B),0)
C=c(Re(C),Im(C),0)
pm=permuteGeneral(c(7,8,9),3)
f=\(D,X) rbind(D,X) |> dist()
ans=NULL
for(i in 1:6){
fn=\(D) (f(D,A)-pm[i,1])^2 + (f(D,B)-pm[i,2])^2 + (f(D,C)-pm[i,3])^2
ans=rbind(ans,optim(c(1,4,7),fn,method='L')$par)
}
ans
for(i in 1:6){
ABCD2V(A,B,C,ans[i,]) |> print()
}
i=1
fn=\(D) (f(D,A)-pm[i,1])^2 + (f(D,B)-pm[i,2])^2 + (f(D,C)-pm[i,3])^2
D=optim(c(1,4,7),fn,method='L')$par
D
ABCD2V(A,B,C,D)
V=rbind(D,A,B,C)
plot3d(V,type='n',col=2,size=1,xlab='x',ylab='y',zlab='z')
for(i in 1:4) points3d(V[i,],col=i,size=8)
connect <- function(x) segments3d(rbind(V[x[1],],V[x[2],]),col=8,lwd=2)
combn(4,2,connect)
Play3d(10)
655:卵の名無しさん (ﾜｯﾁｮｲ 0b12-TCAR [153.137.214.121])
23/07/26 11:43:30.50 UJTnPc6K0.net
尿瓶ジジイさんざん発狂してたのにもう何も言えないのかw
656:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 12:21:03.87 0pDrMthqH.net
大納言60g＋水400mL 圧力(豆料理5分)+保温15分
657:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 17:48:15.79 0pDrMthqH.net
rm(list=ls())
# 底面:△ABC頂点:Dの三角錐の辺の長さからDの座標と体積を求める
# Bを原点,BCをx軸とする
# ヘロンの公式の四面体版
abcd2ABCD=\(AB,BC,CA,DA,DB,DC){
c=AB
a=BC
b=CA
p=DA
q=DB
r=DC
st=c*exp(1i*acos((c^2+a^2-b^2)/(2*c*a)))
s=Re(st)
t=Im(st)
x0=(a^2+q^2-r^2)/(2*a)
P= p^2-(x0-s)^2
Q=q^2-x0^2
y0=(-P+Q+t^2)/(2*t)
z0=sqrt(Q-(-P+Q+t^2)^2/(4*t^2))
D=c(x0,y0,z0)
A=c(s,t,0)
B=c(0,0,0)
C=c(BC,0,0)
v=rbind(A,B,C,D)
Vol=abs(det(rbind(v[1,]-v[4,],v[2,]-v[4,],v[3,]-v[4,])))/6
list(D=D,Vol=Vol)
}
abcd2ABCD(AB=8,BC=5,CA=7,DA=5,DB=7,DC=8)
abcd2ABCD(3,4,5,6,7,8)
abcd2ABCD(4,5,6,7,8,9)
658:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 18:14:20.86 0pDrMthqH.net
三角形を形成しない数値が与えられたときのエラー処理を追加
# 底面:△ABC頂点:Dの三角錐の辺の長さからDの座標と体積を求める
# Bを原点,BCをx軸とする
abcd2ABCD=\(AB,BC,CA,DA,DB,DC){
c=AB
a=BC
b=CA
p=DA
q=DB
r=DC
abc=sort(c(a,b,c))
if(sum(abc[1:2])<=abc[3]) return(NA)
st=c*exp(1i*acos((c^2+a^2-b^2)/(2*c*a)))
s=Re(st)
t=Im(st)
x0=(a^2+q^2-r^2)/(2*a)
P= p^2-(x0-s)^2
Q=q^2-x0^2
y0=(-P+Q+t^2)/(2*t)
if(Q-(-P+Q+t^2)^2/(4*t^2)<0) return(NA)
z0=sqrt(Q-(-P+Q+t^2)^2/(4*t^2))
D=c(x0,y0,z0)
A=c(s,t,0)
B=c(0,0,0)
C=c(BC,0,0)
v=rbind(A,B,C,D)
Vol=abs(det(rbind(v[1,]-v[4,],v[2,]-v[4,],v[3,]-v[4,])))/6
list(D=D,Vol=Vol)
}
659:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 20:45:58.96 0pDrMthqH.net
直径6cmのシリコン製容器
グラニュー糖15mL
水5mL
突沸飛散防止にペーパーでカバー
電子レンジ600Wで90秒
100秒でビターカラメル
660:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/26 20:51:37.21 0pDrMthqH.net
表面積計算と３D画像作成機能を追加。
abcd2ABCD=\(AB,BC,CA,DA,DB,DC,print=FALSE){
c=AB
a=BC
b=CA
p=DA
q=DB
r=DC
abc=sort(c(a,b,c))
if(sum(abc[1:2])<=abc[3]) return(NA)
st=c*exp(1i*acos((c^2+a^2-b^2)/(2*c*a)))
s=Re(st)
t=Im(st)
x0=(a^2+q^2-r^2)/(2*a)
P= p^2-(x0-s)^2
Q=q^2-x0^2
y0=(-P+Q+t^2)/(2*t)
if(Q-(-P+Q+t^2)^2/(4*t^2)<0) return(NA)
z0=sqrt(Q-(-P+Q+t^2)^2/(4*t^2))
D=c(x0,y0,z0)
A=c(s,t,0)
B=c(0,0,0)
C=c(a,0,0)
v=rbind(A,B,C,D)
Vol=abs(det(rbind(v[1,]-v[4,],v[2,]-v[4,],v[3,]-v[4,])))/6
ABC2S3d=\(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
}
S=ABC2S3d(A,B,C)+ABC2S3d(B,C,D)+ABC2S3d(C,D,A)+ABC2S3d(D,A,B)
if(print){
library(rgl)
open3d()
plot3d(v,type='n',col=2,size=1,xlab='x',ylab='y',zlab='z')
for(i in 1:4) points3d(v[i,],col=i,size=8)
connect <- function(x) segments3d(rbind(v[x[1],],v[x[2],]),col=8,lwd=2)
combn(4,2,connect)
}
return(list(D=D,Vol=Vol,S=S))
}
661:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.207])
23/07/26 21:30:36.23 31kfz9Dha.net
904：卵の名無しさん：[sage]：2023/07/26(水) 19:44:23.02 ID:31kfz9Dh
243：１３２人目の素数さん：[sage]：2023/07/25(火) 19:35:19.66 ID:hrc4XW/3
6つの辺の長さが3,4,5,6,7,8である四面体は（ ア ）種類ある。
269：１３２人目の素数さん：[sage]：2023/07/26(水) 18:45:05.79 ID:sev74d4g
>>250
車輪の再発明の神のお告げによれば、39通り
274：１３２人目の素数さん：[sage]：2023/07/26(水) 19:22:12.39 ID:pzlYX2uz
鏡像を同じとみなすなら四面体の各辺に3～8の数字をあてがう方法は30通りしかない
同じと見做さないなら答えは偶数
尿瓶赤っ恥w
662:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 05:45:05.96 QI6KEtmJH.net
# 底面:△ABC頂点:Dの三角錐の辺の長さからDの座標,体積と表面積、底面積、高さを求める
# Bを原点,BCをx軸とする
abcd2ABCD=\(AB,BC,CA,DA,DB,DC,print=FALSE){
c=AB
a=BC
b=CA
p=DA
q=DB
r=DC
abc=sort(c(a,b,c))
if(sum(abc[1:2])<=abc[3]) return(NA)
st=c*exp(1i*acos((c^2+a^2-b^2)/(2*c*a)))
s=Re(st)
t=Im(st)
x0=(a^2+q^2-r^2)/(2*a)
P= p^2-(x0-s)^2
Q=q^2-x0^2
y0=(-P+Q+t^2)/(2*t)
if(Q-(-P+Q+t^2)^2/(4*t^2)<0) return(NA)
z0=sqrt(Q-(-P+Q+t^2)^2/(4*t^2))
D=c(x0,y0,z0)
A=c(s,t,0)
B=c(0,0,0)
C=c(a,0,0)
v=rbind(A,B,C,D)
Vol=abs(det(rbind(v[1,]-v[4,],v[2,]-v[4,],v[3,]-v[4,])))/6
ABC2S3d=\(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
}
S=ABC2S3d(A,B,C)+ABC2S3d(B,C,D)+ABC2S3d(C,D,A)+ABC2S3d(D,A,B)
ABC2S=\(A,B,C) abs(Im((A-C)*Conj(B-C)))/2
s=ABC2S(s+1i*t,0i,a+0i)
if(print){
library(rgl)
open3d()
plot3d(v,type='n',col=2,size=1,xlab='x',ylab='y',zlab='z')
for(i in 1:4) points3d(v[i,],col=i,size=8)
connect <- function(x) segments3d(rbind(v[x[1],],v[x[2],]),col=8,lwd=2)
combn(4,2,connect)
}
return(list(D=D,Vol=Vol,S=S,s=s,h=Vol/s/3))
}
663:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 05:46:26.75 QI6KEtmJH.net
sides=2:7
pm=permuteGeneral(sides)
f=\(x){
ans=abcd2ABCD(x[1],x[2],x[3],x[4],x[5],x[6])
if(length(ans)==1) { return(NA)
}else{ return(ans$D[[3]])}
}
ans=apply(pm,1,f)
Ans=ans[!is.na(ans)]
summary(Ans)
imax=which(ans==max(Ans))
pm[imax,]
imin=which(ans==min(Ans))
pm[imin,]
Ans |> unique() |> length()
fn=\(x){
abcd2ABCD(x[1],x[2],x[3],x[4],x[5],x[6],print=TRUE)
}
fn(pm[imax[1],]) ; aspect3d('iso')
fn(pm[imax[nrow(pm[imax,])],]) ; aspect3d('iso')
fn(pm[imin[1],]) ; aspect3d('iso')
fn(pm[imin[nrow(pm[imin,])],]) ; aspect3d('iso')
664:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 05:50:33.60 QI6KEtmJH.net
朝飯前の問題
(1)長さが1,2,3,4,5,6の6本の線分を組み合わせて四面体は作れるか？
(2)長さが2,3,4,5,6,7の6本の線分を組み合わせて高さが最も大きくなるような四面体を作りたい。
そのときの高さを求めよ。小数点第3位まででよい。
665:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.69.106])
23/07/27 07:11:52.09 hTYahWTka.net
おなじみ朝飯前の大発狂w
666:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 09:13:35.22 QI6KEtmJH.net
Step 1: Assess the internal consistency of the values using the Henderseon-Hasselbach equation:
[H+]=24*PaCO2/[HCO3-]
If the pH and the [H+] are inconsistent, the ABG is probably not valid.
pH2nmol =\(pH) 10^(9-pH)
pH=seq(7.00,7.65,0.05)
cbind(pH,round(pH2nmol(pH)))
667:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 09:14:00.39 QI6KEtmJH.net
> cbind(pH,round(pH2nmol(pH)))
pH
[1,] 7.00 100
[2,] 7.05 89
[3,] 7.10 79
[4,] 7.15 71
[5,] 7.20 63
[6,] 7.25 56
[7,] 7.30 50
[8,] 7.35 45
[9,] 7.40 40
[10,] 7.45 35
[11,] 7.50 32
[12,] 7.55 28
[13,] 7.60 25
[14,] 7.65 22
668:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:50:54.68 QI6KEtmJH.net
Step 3: Is the disturbance respiratory or metabolic?
What is the relationship between the direction of change in the pH and the direction of change in the PaCO2?
In primary respiratory disorders, the pH and PaCO2 change in opposite directions;
in metabolic disorders the pH and PaCO2 change in the same direction.
pH=7.27
pCO2=98.2
HCO3=44.9
ΔPaCO2=pCO2-40
ΔHCO3=HCO3-25
pH2nmol =\(pH) 10^(9-pH) # pH -> [H+] nmol/L
c('pH2[H+]'= pH2nmol(pH) , Henderseon_Hasselbach = 24*pCO2/HCO3)
primary_disorder=if
669:else((pH-7.40)*(pCO2-40)>0,'metabolic','respiratory') primary_disorder Step 4: Is there appropriate compensation for the primary disturbance? Usually, compensation does not return the pH to normal (7.35 – 7.45).
670:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:51:03.00 QI6KEtmJH.net
str=
c(
'Disorder','Expected compensation','Correction factor',
'Metabolic acidosis','PaCO2 = (1.5 x [HCO3-])+8','± 2',
'Acute respiratory acidosis','Increase in [HCO3-]= ∆ PaCO2/10','± 3',
'Chronic respiratory acidosis (3-5 days)','Increase in [HCO3-]= 3.5(∆ PaCO2/10)','',
'Metabolic alkalosis','Increase in PaCO2 = 40 + 0.6(∆HCO3-)','',
'Acute respiratory alkalosis','Decrease in [HCO3-]= 2(∆ PaCO2/10)','',
'Chronic respiratory alkalosis','Decrease in [HCO3-] = 5(∆ PaCO2/10) to 7(∆ PaCO2/10)',''
)
mat=matrix(str,ncol=3,by=T)
df=data.frame(mat[-1,])
colnames(df)=mat[1,]
pH=7.27
pCO2=98.2
HCO3=44.9
ΔPaCO2=pCO2-40
ΔHCO3=HCO3-25
expctd=c(
1.5*HCO3+8,
ΔPaCO2/10,
3.5*ΔPaCO2/10,
40+0.6*ΔHCO3,
2*ΔPaCO2/10,
6*ΔPaCO2/10 # c(5*(ΔPaCO2/10),7*(ΔPaCO2/10))
)
Expected_compensation=
c(
'PaCO2',
'Increase in [HCO3-]',
'Increase in [HCO3-]',
'Increase in PaCO2',
'Decrease in [HCO3-]',
'Decrease in [HCO3-]')
expctd_dat=data.frame(disorder=mat[2:7,1],Expected_compensation,expctd,value=c(expctd+c(0,25,25,40,25,25)))
expctd_dat
c(pH=7.27,
pCO2=98.2,
HCO3=44.9,
ΔPaCO2=pCO2-40,
ΔHCO3=HCO3-25)
If the observed compensation is not the expected compensation,
it is likely that more than one acid-base disorder is present.
671:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:51:25.37 QI6KEtmJH.net
Step 5: Calculate the anion gap (if a metabolic acidosis exists):
AG= [Na+]-( [Cl-] + [HCO3-] )-12 ± 2
A normal anion gap is approximately 12 meq/L.
In patients with hypoalbuminemia, the normal anion gap is lower than 12 meq/L;
the “normal” anion gap in patients with hypoalbuminemia is
about 2.5 meq/L lower for each 1 gm/dL decrease in the plasma albumin concentration
(for example, a patient with a plasma albumin of 2.0 gm/dL would be approximately 7 meq/L.)
If the anion gap is elevated, consider calculating the osmolal gap in compatible clinical
672:situations. Elevation in AG is not explained by an obvious case (DKA, lactic acidosis, renal failure Toxic ingestion is suspected OSM gap = measured OSM – (2[Na+] - glucose/18 – BUN/2.8) The OSM gap should be < 10 pH=7.40 pCO2=40.7 HCO3=25.4 Na=161 Cl=118 Alb=3.4 BUN=64 Glucose=436 Na-Cl-HCO3+2.5*(4.5-Alb) metabolic_acidosis = \(Na,Cl,HCO3,Alb){ AG=Na-Cl-HCO3 AGc=AG+2.5*(4.5-Alb) # albumin_corrected_anion_Gap 4.5g/dL gap_gap=(AGc-12)/(24-HCO3) # AGc Excess/HCO3 Defecit : Gap/Gap Ratio if(AGc<=12) return("norma AG(hyperchloremic) metabolic acidosis") if(AGc>12 & gap_gap<1) co_existence='norma AG(hyperchloremic) metabolic acidosis' if(AGc>12 & gap_gap>2) co_existence='metabolic alkalosis' return(paste("co_existence :",co_existence)) } metabolic_acidosis(Na,Cl,HCO3,Alb)
673:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:51:41.84 QI6KEtmJH.net
Step 6: If an increased anion gap is present, assess the relationship between
the increase in the anion gap and the decrease in [HCO3-].
Assess the ratio of the change in the anion gap (∆AG )
to the change in [HCO3-] (∆[HCO3-]) i.e. ∆AG/∆[HCO3-]
This ratio should be between 1.0 and 2.0
if an uncomplicated anion gap metabolic acidosis is present.
If this ratio falls outside of this range, then another metabolic disorder is present:
If ∆AG/∆[HCO3-] < 1.0,
then a concurrent non-anion gap metabolic acidosis is likely to be present.
If ∆AG/∆[HCO3-] > 2.0,
then a concurrent metabolic alkalosis is likely to be present.
It is important to remember what the expected “normal” anion gap for your patient should be,
by adjusting for hypoalbuminemia (see Step 5, above.)
674:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:53:24.38 QI6KEtmJH.net
Step 3: Is the disturbance respiratory or metabolic?
What is the relationship between the direction of change in the pH and the direction of change in the PaCO2?
In primary respiratory disorders, the pH and PaCO2 change in opposite directions;
in metabolic disorders the pH and PaCO2 change in the same direction.
675:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:55:41.68 QI6KEtmJH.net
∆AG/∆[HCO3-]
This ratio should be between 1.0 and 2.0
これがMarioのICU Bookの記載と異なる点だな。
Astrap関数の数値を変えておこう。
676:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 13:55:47.98 QI6KEtmJH.net
HHSでのデータ
pH=7.40
pCO2=40.7
HCO3=25.4
Na=161
Cl=118
Alb=3.4
BUN=64
Glucose=436
> metabolic_acidosis(Na,Cl,HCO3,Alb)
[1] "co_existence : norma AG(hyperchloremic) metabolic acidosis"
ICU Bookの
677:卵の名無しさん (JP 0H2d-yMPU [146.70.31.18 [上級国民]])
23/07/27 14:00:31.96 QI6KEtmJH.net
# Analysis of metabolic acidosis
metabolic_acidosis = \(Na,Cl,HCO3,Alb){
AG=Na-Cl-HCO3
AGc=AG+2.5*(4.5-Alb) # albumin_corrected_anion_Gap 4.5g/dL
gap_gap=(AGc-12)/(24-HCO3) # AGc Excess/HCO3 Defecit : Gap/Gap Ratio
if(AGc<=12 | 1<=gap_gap | gap_gap<=2) return("norma AG(hyperchloremic) metabolic acidosis")
if(AGc>12 & gap_gap<1) co_existence='norma AG(hyperchloremic) metabolic acidosis'
if(AGc>12 & gap_gap>2) co_existence='metabolic alkalosis'
return(paste("co_existence :",co_existence))
}
metabolic_acidosis(Na=140,Cl=105,HCO3=24,Alb=4.5)
metabolic_acidosis(Na=140,Cl=105,HCO3=24,Alb=3)
metabolic_acidosis(Na=161,Cl=118,HCO3=25.4,Alb=3.4)
678:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.69.80])
23/07/27 17:55:04.18 kPyKYdkra.net
尿瓶ジジイ数学板でも脳内医者とか言われてて草
679:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.68.47])
23/07/27 20:25:40.02 uVmATNWsa.net
尿瓶ジジイ、素人からも素人臭いと言われるw
323：１３２人目の素数さん：[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 : ｎ番目の変数
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：１３２人目の素数さん：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
■■■■■■■■■■■■■■■■■■■■■■■■■■■■
岩田純一は、Ａ医師に暴力をふるった。それを見ていた証人は、多数いる。
磯辺医師、池之上医師、小室医師、関医師、半井英夫医師など。
岩田純一に殴られたＡ医師が、本富士警察署に訴えてやると言ったところ、
半井英夫が、そんなことをしたら、Ａ医師をクビにすると強迫した。
すべて実名なので、警視庁サイバー犯罪対策課に 通報してください。
■■■■■■■■■■■■■■■■■■■■■■■■■■■■
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を結ぶとき２つの弧の交点の経度、緯度を求めよ。
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をとり
４つの球面三角形PES,PSW,PWN,PNEの面積が等しくなるように４分割して
米英中露が分割統治する。
点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
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
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 １３２人目の素数さん 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
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
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
朝のお勤めｗ
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
終末予言 「数年後 魚が絶滅して、漁業が終了します」！！
近未来「岸田総理」が「放射能汚染水」を海に流した！
「神」が怒り「人類の大量絶滅期」（たいりょうぜつめつき）がやってきた！
「放射能汚染水」が「ヨーロッパ」に流れついた！
「ヨーロッパの全国民」が怒り、「日本政府」に
「賠償金」（ばいしょうきん）「９千兆円」を請求（せいきゅう）した！
「日本人」は「世界中の人間から」「憎まれた」！！
「放射能汚染水」が原因で「海の水」が「血のように赤くなった」！
「魚」が全滅した！！魚が滅（ほろ）びた！全部死んだ！！
「「日本の漁業が数年後に終了した」」！！「ヨーロッパの漁業も終了した」！
「岸田総理」は「全人類」から「恨（うら）まれ、憎まれた」！！
ミカエル 「魚が絶滅し「食糧危機」がおこり「人類の３分の２」が死んだ」！
「「人間が人間を食べる時代」がやってくるだろう・・」
「この世界はもうすぐ終わります・・・・」
「神」と「天使」が人類をさばいた！！「七つの大災害」がおこった！！」
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
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
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：１３２人目の素数さん：[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
話をしたのが懇意なナースだったので、
オペ台や検査台の上での同意は無効になる可能性があるんだけど、
「＊＊してくれるなら＊＊＊＊あげる」というベッドでの同意はどう思う？と聞いたら
「娘によく教えておかなくちゃ」と言っていた。（一部脚色ｗ）
医師板における尿瓶ジジイの看護師に関する気色悪いレスをまとめてみました
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のパッケージを使って潜伏期の分布を近似していた。８割おじさんもRとStanで解析。
コードが公開されていたので自分で走らせてみた。
790:卵の名無しさん (ﾜｯﾁｮｲ 5ec7-lnPD [124.47.76.122])
23/08/06 14:02:19.43 BOo6XOy90.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
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
８割おじさんのスクリプトを走らせたことないの？
ひょっとしたらRを使えないないシリツ医かよ？
796:卵の名無しさん (ｱｳｱｳｳｰ Sacb-+YZw [106.128.68.57])
23/08/06 16:31:33.80 sYHzVRqAa.net
65：１３２人目の素数さん：[sage]：2023/08/05(土) 17:16:30.56 ID:AoYbbVsJ
>>61
君たちシダイ卒が求めているようなので貼ってやった。
東大卒による妥当性の検証を希望。
67：１３２人目の素数さん：2023/08/06(日) 00:26:23.88 ID:5u5eLA7R
ちなみに>>64で歌われてる「数学できない偽医者」ってのが>>65のことね。
こういうイカレた爺さんにならないようにね。＞良い子のみなさん
91：１３２人目の素数さん：[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：１３２人目の素数さん：[sage]：2023/08/06(日) 16:28:05.65 ID:P0hZyEmz
>>93
ということは数学できない偽医者は自分ですってことか？ww
100：１３２人目の素数さん：[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：１３２人目の素数さん：[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
お前がうｐしたら俺も見せるよ
もちろんID付きでうｐしろよ
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
ド底辺シリツ卒でもこんな間違いはしない。
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
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)))
}
}