24/04/27 14:39:04.67 TD4Hw7I6.net
>>588
Rだと浮動小数点数の誤差調整が必要。
試作品
is.oncircle <- function(a,b,c,d,tol=1e-16){
tric <- function(a,b,c){
a1=Re(a) ; a2=Im(a)
b1=Re(b) ; b2=Im(b)
c1=Re(c) ; c2=Im(c)
p = (a1^2*(-b2) + a1^2*c2 - a2^2*b2 + a2^2*c2 + a2*b1^2 + a2*b2^2 - a2*c1^2 - a2*c2^2 - b1^2*c2 - b2^2*c2 + b2*c1^2 + b2*c2^2)/(2*(-a1*b2 + a1*c2 + a2*b1 - a2*c1 - b1*c2 + b2*c1))
q = -(a1^2*(-b1) + a1^2*c1 + a1*b1^2 + a1*b2^2 - a1*c1^2 - a1*c2^2 - a2^2*b1 + a2^2*c1 - b1^2*c1 + b1*c1^2 + b1*c2^2 - b2^2*c1)/(2*(-a1*b2 + a1*c2 + a2*b1 - a2*c1 - b1*c2 + b2*c1))
Ce=p+1i*q
r=abs(Ce-a)
c(Center=Ce,Radius=r)
}
abs(sum((tric(a,b,c)-tric(b,c,d))^2)) < tol # all(tric(a,b,c)==tric(b,c,d))
}
> is.oncircle(1+0i,1i,cos(1)+1i*sin(1),cos(2)+1i*sin(2))
[1] TRUE