底辺私立医大を卒業した医者って頭悪いよね？ Part26

960 ：卵の名無しさん (JP 0H03-irN4 [146.70.31.18 [上級国民]])：2023/01/05(木) 15:35:56.29 ID:jdm5LEPgH.net

三角形abcの内接円をoとし中心をpとする。線分ap,bp,cpと円の交点をそれぞれa'b'c'とおく。
aa'=bb'=1/3 cc'=9/7のときの内接円の直径を求めよ。
f=\(x,y,z,Print=FALSE){
a=x+1i*y
b=0i
c=z+0i
o=incircle(a,b,c)
p=o[1]; names(p)=''
r=Re(o[2]) ; names(r)=''
aa.=abs(p-a)-r
bb.=abs(p-b)-r
cc.=abs(p-c)-r
if(Print){
Plot(0,2.5,bty='l')
Polygon(a,b,c)
pt(a,'a') ; pt(b,'b') ; pt(c,'c') ; pt(p,'p')
Cir(p,r,col=2,lwd=2)
seg(a,p) ; seg(b,p) ; seg(c,p)
}
(aa.-1/3)^2+(bb.-1/3)^2+(cc.-9/7)^2
}
opt=optim(c(0.45,0.8,2), \(x) f(x[1],x[2],x[3]) )
x=opt$par[1]
y=opt$par[2]
z=opt$par[3]
a=x+1i*y
b=0i
c=z+0i
incircle(a,b,c)