522:卵の名無しさん (ﾌﾞｰｲﾓ MMdf-/L8v [49.239.65.151])
23/07/16 12:38:48.79 NIpHVD8wM.net
f=\(xy){
if(c==0){return(NA)
}else{
ss=\(x) sum(x^2)
x=xy[1]
y=xy[2]
z=(-d-a*x-b*y)/c
P=c(x,y,z)
((P-A)*ma/ss(P-A) + (P-B)*mb/ss(P-B) + (P-C)*mc/ss(P-C)) |> ss()
}}
opt=optim(colMeans(ABC)[1:2],f)
x=opt$par[1]
y=opt$par[2]
z=(-d-a*x-b*y)/c
c(x,y,z)
523:卵の名無しさん (ﾌﾞｰｲﾓ MMdf-/L8v [49.239.65.151])
23/07/16 12:40:36.04 NIpHVD8wM.net
A(2,0,2)
B(3,0,7)
C(1,6,0)
の座標に惑星があり、質量比は1:2:4である。
A,B,Cの作る平面上の点で３つの惑星からの引力のベクトル和が０になる点の座標を求めよ。
惑星の大きさは考慮しないでよい。有効数字3桁で答えよ。
524:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.68.89])
23/07/16 13:38:51.62 ZSy8rgbma.net
>>503
で、卒業証書と医師免許はどこにあんだよ？
アンタの脳内か？
525:卵の名無しさん (ﾜｯﾁｮｲ cf12-Zr+2 [153.137.214.121])
23/07/16 16:28:43.51 J1dm9s2R0.net
尿瓶ジジイ、あまりのアホさ加減に素人からも脳内医者扱いされるw
526:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 18:15:35.75 xzU/43/d0.net
a=rep(1:6,2)
sim=\() (sample(a)[1:4] |> unique() |> length()) == 2
replicate(1e5,sim()) |> mean()
527:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 18:33:12.25 xzU/43/d0.net
0.030131
528:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 18:33:29.77 xzU/43/d0.net
0.030131
529:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 18:33:51.43 xzU/43/d0.net
a=rep(1:6,2)
sim=\() (sample(a)[1:4] |> unique() |> length()) == 2
replicate(1e6,sim()) |> mean()
530:卵の名無しさん (ﾌﾞｰｲﾓ MM07-/L8v [210.138.178.209])
23/07/16 19:03:06.85 ubGBYOqbM.net
男８人女４人を６個のツインルームに無作為に割り当てる。
男女混合の部屋ができない確率を求めよ。
a=rep(1:6,2)
sim1=\() (sample(a)[1:4] |> unique() |> length()) == 2
k=1e6
replicate(k,sim1()) |> mean()
b=rep(0:1,c(8,4))
sim2=\() (diff(sample(b))[c(1,3,5,7,9,11)]==0) |> all()
replicate(k,sim2()) |> mean()
531:卵の名無しさん (ﾌﾞｰｲﾓ MM07-/L8v [210.138.178.209])
23/07/16 19:33:16.40 ubGBYOqbM.net
library(RcppAlgos)
pm=permuteGeneral(v=0:1,freqs=c(8,4))
f=\(x) (diff(x)[c(1,3,5,7,9,11)]==0) |> all()
apply(pm,1,f) |> mean()
532:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.71.60])
23/07/16 20:36:48.68 RdGkp9Oja.net
尿瓶ジジイダンマリ決め込んでて草
533:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 22:06:31.48 xzU/43/d0.net
BC c1-b1 c2-b2
BP p1-b1 p2-b2
BA a1-b1 a2-b2
BC%o%BP (c1-b1)*(p2-b2) - (c2-b2)*(p2-b2)
BP%o%BA (p1-b1)*(a2-b2) - (p2-b2)*(a1-b1)
inner3 = \(p1,p2,a1,a2,b1,b2,c1,c2){
((c1-b1)*(p2-b2) - (c2-b2)*(p2-b2)) * ((p1-b1)*(a2-b2) - (p2-b2)*(a1-b1)) > 0
&
((a1-c1)*(p2-c2) - (a2-c2)*(p2-c2)) * ((p1-c1)*(b2-c2) - (p2-c2)*(b1-c1)) > 0
}
534:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 22:09:09.75 xzU/43/d0.net
library(RcppAlgos)
pm=permuteGeneral(v=0:1,freqs=c(8,4))
f=\(x) (diff(x)[c(1,3,5,7,9,11)]==0) |> all()
apply(pm,1,f) |> mean()
535:卵の名無しさん (ﾜｯﾁｮｲ 8fc7-/L8v [124.47.76.122])
23/07/16 23:03:40.27 xzU/43/d0.net
inner3 = \(p1,p2,a1,a2,b1,b2,c1,c2){
((b1-a1)*(p2-a2)-(b2-a2)*(p1-a1))*((c1-b1)*(p2-b2)-(c2-b2)*(p1-b1))>0 &
((c1-b1)*(p2-b2)-(c2-b2)*(p1-b1))*((p1-a1)*(c2-a2)-(p2-a2)*(c1-a1))>0
}
inner3(0.3,0.1,0,0,4,0,4,3)
inner3(1,2,0,0,4,0,4,3)
#((b1-a1)*(p2-a2)-(b2-a2)*(p1-a1))
#((c1-b1)*(p2-b2)-(c2-b2)*(p1-b1))
#((p1-a1)*(c2-a2)-(p2-a2)*(c1-a1))
536:卵の名無しさん (ﾌﾞｰｲﾓ MM97-/L8v [210.148.125.249])
23/07/16 23:33:26.12 1XefIVSwM.net
収束時期のシミュレーションなら可能。
SEIRモデルで有病率を１％に固定して、集団のサイズを変化させてシミュレーションしてみたけどピークは変わらないな。
このモデルでは集会規模の大小には影響されないということになるな。
URLﾘﾝｸ(i.imgur.com)
有病率を変化させて流行の変遷をグラフにすると、
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
有病率を４０％くらいに引き上げるとオリンピックのときには流行が収束していることになるw
537:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.68.57])
23/07/17 00:20:35.46 weqbEJola.net
尿瓶ジジイ数学板でご披露しても馬鹿にされるだけw
538:卵の名無しさん (ﾌﾞｰｲﾓ MM97-/L8v [210.148.125.249])
23/07/17 00:55:22.27
539:vAe3YqQlM.net
540:卵の名無しさん (ﾌﾞｰｲﾓ MM97-/L8v [210.148.125.249])
23/07/17 08:30:18.02 vAe3YqQlM.net
URLﾘﾝｸ(www.calculator.net)
541:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.69.141])
23/07/17 14:13:18.39 Opqf8FtRa.net
おい尿瓶ジジイ
素人にすら偽医者とか言われてバカにされるのがそんなに楽しいか？w
542:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.71.171])
23/07/17 20:00:07.17 gHudq7e2a.net
URLﾘﾝｸ(hissi.org)
うんこ製造機、週に1回ほど大発狂
543:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 11:41:41.14 wE+74BzaH.net
> a=sqrt(2)
> a^2
[1] 2
> a^2-2
[1] 4.440892e-16
> a^2==2
[1] FALSE
544:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 11:53:48.34 wE+74BzaH.net
install.packages("Rcpp", dependencies = TRUE, INSTALL_opts = '--no-lock')
545:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 11:59:49.45 wE+74BzaH.net
The easiest way to avoid this issue is before installing any package run the line below
options("install.lock"=FALSE)
546:卵の名無しさん (ﾜｯﾁｮｲ cf12-Zr+2 [153.137.214.121])
23/07/18 15:12:22.67 GM2pBAaY0.net
>>527
おい尿瓶ジジイ
いつになったら医師免許と卒業証書出せるんだよ？
547:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.70.245])
23/07/18 17:43:26.52 33w597Qga.net
バカすぎて高校生からも脳内学歴扱いだね尿瓶ジジイw
548:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 19:42:42.46 wE+74BzaH.net
in4 = \(P,A,B,C,D){
f =\(P,A,B,C){
cross = \(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])
v=cross(B-A,C-A)
a=v[1]
b=v[2]
c=v[3]
d=-a*A[1]-b*A[2]-c*A[3]
sum(v*P) + d
}
f(P,A,B,C)*f(D,A,B,C)>0 &
f(P,B,C,D)*f(A,B,C,D)>0 &
f(P,C,D,A)*f(B,C,D,A)>0 &
f(P,D,A,B)*f(C,D,A,B)>0
}
549:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 19:43:15.42 wE+74BzaH.net
A=c(1,1,1)
B=c(1,-1,-1)
C=c(-1,1,-1)
D=c(-1,-1,1)
ABCD=rbind(A,B,C,D)
open3d()
plot3d(ABCD,xlab='x',ylab='y',zlab='z',type='n',aspect = 1,
xlim=c(-2,2),ylim=c(-2,2),zlim=c(-2,2))
points3d(ABCD[1,],size=3,col=4)
points3d(ABCD[2,],size=3,col=4)
points3d(ABCD[3,],size=3,col=4)
points3d(ABCD[4,],size=3,col=4)
segments3d(ABCD[c(1,2),],col=4)
segments3d(ABCD[c(1,3),],col=4)
segments3d(ABCD[c(1,4),],col=4)
segments3d(ABCD[c(2,3),],col=4)
segments3d(ABCD[c(2,4),],col=4)
segments3d(ABCD[c(3,4),],col=4)
sim=\(){
P=runif(3,-2,2)
if(in4(P,A,B,C,D)){
points3d(P,size=1,col='black')
}else{
# points3d(P,col='red')
}
}
del=replicate(5e4,sim())
550:卵の名無しさん (JP 0H6f-QoCo [217.138.252.98 [上級国民]])
23/07/18 20:02:14.18 wE+74BzaH.net
俺は理Iを蹴って医科歯科に言った口だけど。
医科歯科を蹴って慶応に入学した医師がいても
医科歯科合格の証拠を出せとか言わないな。
東大合格通知の画像を出せとかいう輩って、東大もしくは国立を落ちたのだろうと憶測。
551:卵の名無しさん (JP 0H6f-gMG4 [217.138.252.98 [上級国民]])
23/07/18 20:03:14.58 wE+74BzaH.net
>>528
んで、どこの国立を落ちたの？
シリツ卒なんだろ？
552:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.71.77])
23/07/19 00:23:48.07 iFTHmQL+a.net
>>532
あんたがアホすぎてそれしか証明しようがないからだろw
553:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.69.68])
23/07/19 00:32:20.24 WXv06B+Pa.net
誰も
554:相手にしてないからこそ医者板じゃここでしか発狂できないんだろうけどなw 高校生にすらバカにされる始末だし
555:卵の名無しさん (JP 0H6f-gMG4 [217.138.212.122 [上級国民]])
23/07/19 07:31:36.73 qC76rOplH.net
直径６ｃｍのシリコンカップ
グラニュー糖小さじ３
水小さじ１
カップをペーパーで覆って電子レンジ６００W１分３０秒
追加の水や熱湯はなし
556:卵の名無しさん (ﾜｯﾁｮｲ cf12-Zr+2 [153.137.214.121])
23/07/19 12:44:37.47 7WAeHWTk0.net
>>533
んで、どこの国立受けられなかったの？
557:卵の名無しさん (ﾜｯﾁｮｲ 6724-QoCo [149.50.210.17 [上級国民]])
23/07/20 07:49:05.42 FCPwP+wT0.net
rm(list=ls())
A=c(-2,0,0)
B=c(2,1,1)
f=\(P){
dist(rbind(P,A))+dist(rbind(P,B))
}
sim=\(){
th=runif(1,-pi,pi)
x=runif(1)*cos(th)
y=runif(1)*sin(th)
z=runif(1,0,2)
f(c(x,y,z))
}
d=replicate(1e6,sim())
max(d)
g=\(){
x=runif(1)
y=runif(1)
if(x^2+y^2>1){return(0)
}else{
z=runif(1,0,2)
return(f(c(x,y,z)))
}
}
r=replicate(1e6,g())
BEST::plotPost(r)
558:卵の名無しさん (ﾜｯﾁｮｲ 6724-QoCo [149.50.210.17 [上級国民]])
23/07/20 08:16:55.26 FCPwP+wT0.net
h=\(th,z){
if(abs(th)>pi | z<0 | 2<z) r=0
else{
r=f(c(cos(th),sin(th),z))
}
return(r)
}
h=Vectorize(h)
optim(c(runif(1,-pi,pi),runif(1,0,2)), \(x) h(x[1],x[2]),control = list(fnscale=-1))
h(-pi/2,2)
559:卵の名無しさん (ﾜｯﾁｮｲ cf12-Zr+2 [153.137.214.121])
23/07/20 11:04:12.21 VIVJcW6Q0.net
おい尿瓶ジジイ
いつになったら医師免許と卒業証書出すんだよハゲ
560:卵の名無しさん (ﾜｯﾁｮｲ 6724-QoCo [149.50.210.17 [上級国民]])
23/07/20 11:41:22.96 FCPwP+wT0.net
平面上の4点A,B,C,Dが存在している。
Aから反時計回りにのこりの3点を線分で結んだときに凸四角形ができるようにする
in3 <- \(P,A,B,C){
opc <- \(a,b) Re(a)*Im(b)-Im(a)*Re(b)
sum(opc(B-A,P-A)>0,opc(C-B,P-B)>0,opc(P-A,C-A)>0)%%3==0
}
is.concave=\(A,B,C,D){
c(in3(A,B,C,D),in3(B,C,D,A),in3(C,D,A,B),in3(D,A,B,C)) |> any()
}
Do =\(){
sq=cbind(runif(4),runif(4))
sq=round(sq,2)
A=sq[1,1]+1i*sq[1,2]
B=sq[2,1]+1i*sq[2,2]
C=sq[3,1]+1i*sq[3,2]
D=sq[4,1]+1i*sq[4,2]
ABCD=c(A,B,C,D)
G=mean(ABCD)
if(is.concave(A,B,C,D)){
par(mfrow=c(2,2))
Plot(0,1,zero=FALSE) ; pta(A) ; pta(B) ; pta(C) ; pta(D) ; pt(G,'.')
Plot(0,1,zero=FALSE) ; pta(A) ; pta(B) ; pta(C) ; pta(D) ; pt(G,'.')
Polygon(A,B,C,D)
Plot(0,1,zero=FALSE) ; pta(A) ; pta(B) ; pta(C) ; pta(D) ; pt(G,'.')
Polygon(A,C,B,D,Col=2)
Plot(0,1,zero=FALSE) ; pta(A) ; pta(B) ; pta(C) ; pta(D) ; pt(G,'.')
Polygon(A,B,D,C,Col=3)
}else{
layout(1)
Arg(ABCD-G) |> order() -> ord
tmp=c(ord,ord)
i=which(tmp==1)[1]
ord=tmp[i:(i+3)]
cat(LETTERS[ord],'\n')
Plot(0,1,zero=FALSE) ; pta(A) ; pta(B) ; pta(C) ; pta(D) ; pt(G,'.')
Polygon(ABCD[ord])
}
invisible(ABCD)
}
(Do())
561:卵の名無しさん (ﾜｯﾁｮｲ 6724-QoCo [149.50.210.17 [上級国民]])
23/07/20 11:47:55.38 FCPwP+wT0.net
URLﾘﾝｸ(readcrx-2.github.io)
アップデートされている。
562:卵の名無しさん (JP 0H6f-gMG4 [217.138.212.98 [上級国民]])
23/07/20 12:45:03.03 1lBdubx6H.net
>>537
俺は二期校最後の年の受験だから
東大理Iと医科歯科を受験。
医科歯科の合格発表日が東大合格者の健康診断日より後だったので東大で健康診断を受けたよ。
んで、あんたはどこの国立を落ちたの？
563:卵の名無しさん (JP 0Hdf-Zr+2 [61.199.225.248])
23/07/20 13:29:54.35 iz25aXn4H.net
>>543
あんた5chしかしてない脳内医者だろw
564:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:31:51.28 1lBdubx6H.net
Astrap <- \(pH,pCO2,HCO3){
e=1e-12
pH_l=7.36-e
pH_u=7.44+e
pCO2_l=36-e
pCO2_u=44+e
HCO3_l=22-e
HCO3_u=26+e
expected_CO2=NULL
expected_HCO3=NULL
status=NULL
HH =\(pCO2,HCO3,pK=6.1) round(pK+log10(HCO3/(pCO2*0.03)),2) # Henderson-Hasselbalch
if(pH_l<pH & pH<pH_u){
if(pCO2>pCO2_u){
primary="mixed respiratory aci
565:dosis" additional="and metabolic alkalosis" } if(pCO2<pCO2_l){ primary="mixed respiratory alkalosis" additional="and metabolic acidosis" } } if(pCO2_l<pCO2 & pCO2<pCO2_u){ if(pH>pH_u){ primary="mixed metaboic alkalosis" additional="and respiratory acidosis" } if(pH<pH_l){ primary="mixed metabolic acidosis" additional="and respiratory alkalosis" } }
566:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:32:06.85 1lBdubx6H.net
if(pH<pH_l & pCO2<pCO2_l){
primary="primary metabolic acidosis"
expected_pCO2=40 + 1.2*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis"
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis"
expected_pH=HH(expected_pCO2,HCO3)
}
if(pH>pH_u & pCO2>pCO2_u){
primary="primary metabolic alkalosis"
expected_pCO2=40 + 0.7*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis"
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis"
expected_pH=HH(expected_pCO2,HCO3)
}
if(pH<pH_l & pCO2>pCO2_u){
primary="primary respiratory aciosis"
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.1*(pCO2-40)
if(status=="chronic") expected_HCO3=24+0.4*(pCO2-40)
if(HCO3<expected_HCO3) additional="incomplete renal response"
if(HCO3>expected_HCO3) additional="secondary metabolic alkalosis"
expected_pH=HH(pCO2,expected_HCO3)
}
if(pH>pH_u & pCO2<pCO2_l){
primary="primary respiratory alkalosis"
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.2*(pCO2-40)
if(status=="chronic") expected_HCO3=24+0.4*(pCO2-40)
if(HCO3<expected_HCO3) additional="secondary metabolic acidosis"
if(HCO3>expected_HCO3) additional="incomplete renal response"
expected_pH=HH(pCO2,expected_HCO3)
}
cat(status,primary,'\n')
cat(additional,'\n\n')
if(!is.null(expected_CO2)|!is.null(expected_HCO3)){
return(c(expected_CO2=expected_CO2,expected_HCO3=expected_HCO3,expected_pH=expected_pH))
}
}
567:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.68.80])
23/07/20 17:34:39.60 //M/SQ0Ba.net
こっちではダンマリがジジイw
568:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:35:32.79 1lBdubx6H.net
パソコンだとクラウド上のファイルをgrepして必要なファイルを取り出すのは簡単だけど、
スマホだと面倒。ローカルに置くと場所を取るし。
有用なファイルはwebで公開しておくと、他の人も流用できるし、
自分で探すのも楽で( ・∀・)ｲｲ!!
569:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:36:08.07 1lBdubx6H.net
>>547
んで、あんたはどこの国立を落ちたの？
570:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:46:16.77 1lBdubx6H.net
医師板だのと行数制限小さいので
動脈血ガス分析のスクリプトは１レスでは投稿できないのが不便。
571:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:50:14.79 1lBdubx6H.net
>>544
月曜日は祝日だったので内視鏡バイトなしだった。
明日はEGD11件の予定。
俺は理Iを蹴って医科歯科に言った口だけど。
医科歯科を蹴って慶応に入学した医師がいても
あっそう、程度にしか思わんな。医科歯科合格の証拠を出せとか言わないね。
東大卒業していたら東大蹴って他大学に行ったという人がいてもなんとも思わんだろうに。
んで、どこの国立を落ちたの？
572:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/20 17:52:06.71 1lBdubx6H.net
>>547
あんたはどこの国立を落ちたの？
にはダンマリかよ？
573:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.71.42])
23/07/20 17:58:25.88 goI3Ka80a.net
>>552
尿瓶ジジイの脳内学歴()に受けるもクソもないだろタコw
574:卵の名無しさん (ﾜｯﾁｮｲ cf12-Zr+2 [153.137.214.121])
23/07/20 17:59:27.26 VIVJcW6Q0.net
>>551
そこまで言うなら医科歯科の卒業証書と医師免許さっさと出せよ
なんで出せないんだ？
575:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:33:51.72 Mz3Ba9iAH.net
>>554
あんたがアップロードする番だぞ。
俺は手間をかけて東大合格通知のコピー画像をアップロードしたぞ。
どこの国立を落ちたの？
576:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:34:14.39 Mz3Ba9iAH.net
>>547
あんたはどこの国立を落ちたの？
にはダンマリかよ？
577:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.68.164])
23/07/20 18:37:52.73 IF4YK+Pma.net
>>556
全部上げない限り脳内学歴の脳内医者でいいのか？
そういう結論になったから医者板の他のスレで発狂できなくなったんだろ？
578:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:41:09.80 Mz3Ba9iAH.net
>>557
んで、どこの国立を落ちたの？
579:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:42:38.75 Mz3Ba9iAH.net
東大合格通知って健康診断の受診票を兼ねていた。
大きさはハガキ大で公印すら押されていなくてありがたみがなかった。
580:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:44:08.44 Mz3Ba9iAH.net
>>557
んで、どこの国立を落ちたの？
にはだんまりかよ？
東大合格通知の書式すら知らなかったから
東大非合格は確定しているが。
シリツ卒なんだろ？
母校に誇りはないのかよ？？
581:卵の名無しさん (JP 0H6f-mROX [217.138.252.98 [上級国民]])
23/07/20 18:44:45.68 Mz3Ba9iAH.net
>>554
どこの国立を落ちたの？に何で答えられないんだ？
582:卵の名無しさん (ｱｳｱｳｳｰ Sa17-Zr+2 [106.128.68.164])
23/07/20 18:56:16.25 IF4YK+Pma.net
>>560
あんたの脳内学歴なんか知らんがな
583:w
584:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/21 06:23:16.54 8CtH7l/hH.net
tetragon <- function(A=NULL,B=NULL,C=NULL,D=NULL,print=TRUE){
in3 <- function(P,A,B,C){
opc <- function(a,b) Re(a)*Im(b)-Im(a)*Re(b)
sum(opc(B-A,P-A)>0,opc(C-B,P-B)>0,opc(P-A,C-A)>0)%%3==0
}
A=ifelse(is.null(A),runif(1)+1i*runif(1),A)
B=ifelse(is.null(B),runif(1)+1i*runif(1),B)
C=ifelse(is.null(C),runif(1)+1i*runif(1),C)
D=ifelse(is.null(D),runif(1)+1i*runif(1),D)
if(print){
par(mfrow=c(2,3))
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,B,C,D,Col=8)
Polygon(0,1,1+1i,1i)
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,B,D,C,Col=2)
Polygon(0,1,1+1i,1i)
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,C,B,D,Col=3)
Polygon(0,1,1+1i,1i)
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,C,D,B,Col=4)
Polygon(0,1,1+1i,1i)
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,D,B,C,Col=5)
Polygon(0,1,1+1i,1i)
Plot(0,1,axes=FALSE,zero=FALSE)
pt(A,'A') ; pt(B,'B') ; pt(C,'C') ; pt(D,'D')
Polygon(A,D,C,B,Col=6)
Polygon(0,1,1+1i,1i)
layout(1)
}
585:卵の名無しさん (JP 0H6f-QoCo [217.138.212.98 [上級国民]])
23/07/21 06:24:22.90 8CtH7l/hH.net
if(any(c(in3(A,B,C,D),in3(B,C,D,A),in3(C,D,A,B),in3(D,A,B,C)))){
ABCD=ADCB=ABC2S(A,B,C)+ABC2S(A,C,D)
ABDC=ACDB=ABC2S(A,B,C)+ABC2S(B,D,C)
ACBD=ADBC=ABC2S(A,C,D)+ABC2S(C,B,D)
if(print){
print(data.frame(ABCD=ABCD,ABDC=ABDC,ACBD=ACBD))
print(data.frame(ADCB=ADCB,ACDB=ACDB,ADBC=ADBC))}
invisible(c(ABCD,ABDC,ACBD))
}else{
v=c(A,B,C,D)
m=mean(v)
i=order(Arg(v-m))
S=ABC2S(v[i[1]],v[i[2]],v[i[3]]) +
ABC2S(v[i[3]],v[i[4]],v[i[1]])
re=LETTERS[i]
re2=rep(re,2)
j=min(which(re2=='A'))
ir=order(Arg(v-m),decreasing = TRUE)
rer=LETTERS[ir]
rer2=rep(rer,2)
jr=min(which(rer2=='A'))
lang=str2lang(paste0('data.frame(',paste0(re2[j:(j+3)],collapse=''),'=S)'))
if(print) print(eval(lang))
langr=str2lang(paste0('data.frame(',paste0(rer2[jr:(jr+3)],collapse=''),'=S)'))
if(print) print(eval(langr))
invisible(S)
}
}
L=10
k=1e5
replicate(k,L^2*tetragon(print=FALSE) |> max()) -> S
hist(S)
mean(S>1/4*L^2)
mean(S)
median(S)
586:卵の名無しさん (JP 0H6f-QoCo [217.138.212.122 [上級国民]])
23/07/21 19:44:17.26 stxMOQIEH.net
> replicate(k,tetragon(print=FALSE)|>length()) -> len
> mean(len==1)
[1] 0.694586
587:卵の名無しさん (ｱｳｱｳｳｰ Sa17-tdBp [106.128.70.115])
23/07/21 20:47:56.37 qBzKS7WRa.net
おい尿瓶クソジジイ
さっさと卒業証書と医師免許出せや
588:卵の名無しさん (JP 0Hd3-zEeY [217.138.252.98 [上級国民]])
23/07/22 05:23:29.31 FYFO7qrzH.net
識者の尿瓶チンパポンコツフェチへの評価
61 １３２人目の素数さん sage 2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57
お前が一番の馬鹿だよ
相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく
純粋に高校数学の話だけをしなさい
589:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 05:46:22.35 54NVUurG0.net
fn <- function(
ai=45, # 北緯°
ak=145, # 東経°
bi=45, # 北緯°(南緯はー°)
bk=-125){ # 東経°(西経はー°)
R=6378 # 地球の半径km
u=pi/180
A=R*c(cos(ai*u)*cos(ak*u),cos(ai*u)*sin(ak*u),sin(ai*u))
B=R*c(cos(bi*u)*cos(bk*u),cos(bi*u)*sin(bk*u),sin(bi*u))
AB=sqrt(sum((A-B)^2))
asin(AB/(2*R))*2*R
}
fn()
590:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 07:05:52.70 54NVUurG0.net
子曰、
高校数学の質問スレ Part430
ｽﾚﾘﾝｸ(math板:61番)
61 名前：１３２人目の素数さん[sage] 投稿日：2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57
お前が一番の馬鹿だよ
相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく
純粋に高校数学の話だけをしなさい
591:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.115])
23/07/22 08:08:31.15 SRPzw5T1a.net
自分の気に食わないレスは全員同じに見える病気みたいだね
尿瓶ジジイもいちいち発狂してる分際でw
592:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.115])
23/07/22 08:46:35.79 SRPzw5T1a.net
尿瓶ジジイ、ウクライナの首都を間違えたことを指摘されて発狂中ww
593:卵の名無しさん (JP 0Hd3-zEeY [217.138.252.98 [上級国民]])
23/07/22 12:56:19.29 FYFO7qrzH.net
東大合格者から国立不合格者へのありがたい助言がこれ！
高校数学の質問スレ Part430
ｽﾚﾘﾝｸ(math板:61番)
61 名前：１３２人目の素数さん[sage] 投稿日：2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57 お前が一番の馬鹿だよ 相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく 純粋に高校数学の話だけをしなさい
595:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 15:54:26.67 54NVUurG0.net
s=c('東京','35.681382','139.76608399999998',
'有楽町','35.675069','139.763328',
'新橋','35.665498','139.75964',
'浜松町','35.655646','139.756749',
'田町','35.645736','139.74757499999998',
'高輪ゲートウェイ','35.6355406','139.7407245',
'品川','35.630152','139.74044000000004',
'大崎','35.6197','139.72855300000003',
'五反田','35.626446','139.72344399999997',
'目黒','35.633998','139.715828',
'恵比寿','35.64669','139.710106',
'渋谷','35.658517','139.70133399999997',
'原宿','35.670168','139.70268699999997',
'代々木','35.683061','139.702042',
'新宿','35.690921','139.70025799999996',
'新大久保','35.701306','139.70004399999993',
'高田馬場','35.712285','139.70378200000005',
'目白','35.721204','139.706587',
'池袋','35.728926','139.71038',
'大塚','35.731401','139.72866199999999',
'巣鴨','35.733492','139.73934499999996',
'駒込','35.736489','139.74687500000005',
'田端','35.738062','139.76085999999998',
'西日暮里','35.732135','139.76678700000002',
'日暮里','35.727772','139.770987',
'鶯谷','35.720495','139.77883700000007',
'上野','35.713768','139.77725399999997',
'御徒町','35.707438','139.774632',
'秋葉原','35.698683','139.77421900000002',
'神田','35.69169','139.77088300000003')
596:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 15:55:02.21 54NVUurG0.net
mat=matrix(s,ncol=3,by=TRUE)
mat[,2] |> as.numeric() -> lat
lat=round(lat,4)
mat[,3] |> as.numeric() -> lon
lon=round(lon,4)
df=data.frame(station=mat[,1],lat=lat,lon=lon)
d=cbind(lat,lon)
d
Hubeny=\(ai,ak,bi,bk){
u=pi/180
Dy=(ai-bi)*u
Dx=(ak-bk)*u
P=(ai/2+bi/2)*u
Rx=6378.137 # 長半径（赤道半径）
Ry=6356.752314245 # 長半径(極半径)
f=1/298.257223563 # 扁平率 （Rx-Ry)/Rx
E2=(Rx^2-Ry^2)/Rx^2 # 離心率^2
W=sqrt(1-E2*sin(P)^2)
M=Rx*(1-E2)/W^3 # 子午線曲率半径
N=Rx/W # 卯酉線曲率半径
D=sqrt((Dy*M)^2 + (Dx*N*cos(P))^2)
D
}
f=\(xy){
ans=0
for(i in 1:nrow(d)){
ans=ans+Hubeny(xy[1],xy[2],d[i,1],d[i,2])
}
ans
}
opt=optim(colMeans(d),f)
opt$par
597:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.180])
23/07/22 16:15:46.60 f5YfabY7a.net
尿瓶ジジイ、またコピペの発作で発狂中w
598:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 18:32:29.13 54NVUurG0.net
道具箱に追加しておこう。
緯度・経度がらみの計算に使えそう。
# 緯度経度から楕円地球上の距離を計算
Hubeny=\(ai,ak,bi,bk){# ai:A地点の緯度 ak:A地点の経度
u=pi/180
Dy=(ai-bi)*u
Dx=(ak-bk)*u
P=(ai/2+bi/2)*u
Rx=6378.137 # 長半径（赤道半径）
Ry=6356.752314245 # 長半径(極半径)
f=1/298.257223563 # 扁平率 （Rx-Ry)/Rx
E2=(Rx^2-Ry^2)/Rx^2 # 離心率^2
W=sqrt(1-E2*sin(P)^2)
M=Rx*(1-E2)/W^3 # 子午線曲率半径
N=Rx/W # 卯酉線曲率半径
D=sqrt((Dy*M)^2 + (Dx*N*cos(P))^2)
D
}
599:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 18:33:26.70 54NVUurG0.net
東大合格者から国立不合格者へのありがたい助言がこれ！
高校数学の質問スレ Part430
ｽﾚﾘﾝｸ(math板:61番)
61 名前：１３２人目の素数さん[sage] 投稿日：2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57
お前が一番の馬鹿だよ
相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく
純粋に高校数学の話だけをしなさい
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
ですって、大切なことなので繰り返しましたｗ
600:卵の名無しさん (ﾜｯﾁｮｲ f324-yMPU [149.50.210.17 [上級国民]])
23/07/22 19:35:59.67 54NVUurG0.net
G7の首都と緯度経度は次の通りとする。
オタワ -カナダ 緯度: 45.4208, 経度: -75.690
パリ -フランス 緯度: 48.8032, 経度: 2.3511
ベルリン-ドイツ 緯度: 52.5186, 経度: 13.4081
ローマ -イタリア 緯度: 41.8905, 経度: 12.4942
東京 - 日本 緯度: 35.6894, 経度: 139.692
ロンドン-英国 緯度: 51.5072, 経度: -0.1275
ワシントン-米国 緯度: 38.9041, 経度: -77.0171
７カ国から集まるのに移動距離の総和が最小になる地点の緯度と経度を計算せよ。
その地点に最も近い国はどこか述べよ。
601:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.2])
23/07/22 21:20:02.90 r8HmaS4wa.net
>>577
尿瓶ジジイって自分の気に食わないレス全員が同じに見える病気なんだなw
602:卵の名無しさん (JP 0Hd3-zEeY [217.138.212.98 [上級国民]])
23/07/23 06:14:06.10 W8dTuzCcH.net
東大合格者から国立不合格者へのありがたい助言がこれ！
高校数学の質問スレ Part430
ｽﾚﾘﾝｸ(math板:61番)
61 名前：１３２人目の素数さん[sage] 投稿日：2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57
お前が一番の馬鹿だよ
相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく
純粋に高校数学の話だけをしなさい
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
603:卵の名無しさん (JP 0Hd3-zEeY [217.138.212.98 [上級国民]])
23/07/23 06:15:17.44 W8dTuzCcH.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
抗パ剤を高パ剤と誤変換で投稿するのは単純ミスだが
抗ドーパミン薬と投稿するのはタイプミスでも誤変換でもないよな。
604:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.2])
23/07/23 08:42:37.98 SFdfjNW1a.net
尿瓶ジジイまともな反論は何一つできないからなw
605:卵の名無しさん (ｽﾌﾟｰﾌﾟ Sd73-leiM [1.73.154.206])
23/07/23 09:40:19.03 Y55+EhIqd.net
>>581
久しぶり見に来たらコイツまだコレ言ってるのかwww
頭悪いお前にはこれしか言えることないもんな
これもわざとtypoしたって論破されてるのにwwwジメジメした性格してて本当気持ち悪い奴だよお前は
統合失調症はまだ治らんのか？
606:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.70.2])
23/07/23 10:47:31.84 SFdfjNW1a.net
現在は下手に突っ込まれないようここでお得意の脳内医療ネタ()もろくに書き込めず数学板で高校生相手に発狂してはバカにされるの繰り返し、憂さ晴らしにここでコピペで発狂といった具合w
607:卵の名無しさん (ﾜｯﾁｮｲ b9c7-2oYC [124.47.76.122])
23/07/23 13:22:50.59 UToIyy7I0.net
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
608:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.123])
23/07/23 13:39:30.06 LHFf3ysQa.net
尿瓶ジジイは気に食わないレスは全員同じに見える病気にかかっていて全く煽りになってないことに気づかないホームラン級のバカw
609:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.123])
23/07/23 13:42:46.07 LHFf3ysQa.net
替え歌荒らしにすら偽医者だのクズ呼ばわりされる始末w
610:卵の名無しさん (ﾜｯﾁｮｲ b9c7-2oYC [124.47.76.122])
23/07/23 16:22:10.87 UToIyy7I0.net
APIMateR apimater:URLﾘﾝｸ({$host})([25]ch\.net|bbspink\.com)$]}
611:卵の名無しさん (ﾜｯﾁｮｲ b9c7-2oYC [124.47.76.122])
23/07/23 16:23:08.21 UToIyy7I0.net
東大合格者から国立不合格者へのありがたい助言がこれ！
高校数学の質問スレ Part430
ｽﾚﾘﾝｸ(math板:61番)
61 名前：１３２人目の素数さん[sage] 投稿日：2023/07/21(金) 19:18:29.80 ID:KcuTGlEs
>>57
お前が一番の馬鹿だよ
相手しなけりゃいいのに、お前自身が1日だけで30以上も書き込みして、余計にこのスレが腐っていく
純粋に高校数学の話だけをしなさい
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
＞お前が一番の馬鹿だよ
612:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.68.155])
23/07/23 16:57:50.41 mZ+wxT/Ea.net
尿瓶チンパンジジイコピペしか能がないのかw
613:卵の名無しさん (ﾜｯﾁｮｲ 9158-leiM [106.73.2.65])
23/07/23 17:06:09.02 a+W/K+U40.net
>>589
コイツが何で東大合格者だと思うんだよ
全く根拠ないだろ
というか、この替え歌君にボコボコにされて悔しかったのか？www
ボコボコにされた相手にはコピペで対抗
いかにも頭の弱いお前の発想だよ
ちなみに一番の馬鹿はお前だぞwww
614:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.150])
23/07/23 17:33:24.68 lqQfOF4fa.net
言い返せないとコピペで発狂
マジでチンパンジーw
615:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.18])
23/07/23 22:53:38.75 33hKJ21QM.net
救急搬送と新入院でインセンティブ4.5諭吉。
久々に両側にGGOのあるコロナ肺炎のCTをみた。
616:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.18])
23/07/23 22:54:15.86 33hKJ21QM.net
>>591
東大合格に速攻で発狂！！
617:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.18])
23/07/24 05:52:45.61 y
618:a9IubfJM.net
619:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.18])
23/07/24 06:27:59.94 ya9IubfJM.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) return("norma AG(hyperchloremic) metabolic acidosis")
if(AGc>12 & gap_gap<1) co_existence='norma AG(hyperchloremic) metabolic acidosis'
if(AGc>12 & gap_gap>1) 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)
620:卵の名無しさん (ｽﾌﾟｰﾌﾟ Sd73-leiM [1.73.153.28])
23/07/24 07:50:00.52 w6igBrqid.net
>>593
こういう、無駄に略語使うとこが素人臭いんだよなwww
意味がわかって使ってないだろ？
だからいつも突っ込まれるんだぞ
621:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.190])
23/07/24 09:21:55.47 Ldvz/BNaa.net
胆汁ドレナージで素人バレバレなのでw
というか医療事務じゃないと出てこないw
622:卵の名無しさん (ｱｳｱｳｳｰ Sabd-TCAR [106.128.71.126])
23/07/24 10:22:09.16 zPtCn08Ta.net
言い返せないとしばらくダンマリ決め込んでコピペで発狂
死ぬまで続ける気か？
623:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.5])
23/07/24 11:56:46.19 Day8w5dSM.net
胆汁ドレナージは普通に使うが。
こういうことを言う医師はまともではないね。
＞悪性症候群で抗ドーパミン薬止めないなんて当たり前だろ。
ですってｗｗｗ
＞当たり前だろ
＞当たり前だろ
＞当たり前だろ
とドヤ顔で大間違いを投稿！
624:卵の名無しさん (ﾌﾞｰｲﾓ MMad-2oYC [210.138.6.5])
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だと言ったり、
わざと間違えたと言ったり、見苦しいよなぁ。