24/04/04 09:54:04.69 WjlV2HHp0.net
# URLﾘﾝｸ(ja.wikipedia.org)パスカルの蝸牛形
# 実数x,yが(x^2+y^2-2*x)^2=x^2+y^2を満たすときx+yおよびの取りうる値の範囲を求めよ
rm(list=ls())
f=\(x,y) (x^2+y^2-2*x)^2 - (x^2+y^2)
x=y=seq(-3,3,le=500)
z=outer(x,y,f)
contour(x,y,z)
tol=5e-4
k=7500
xy=NULL
for(x in runif(k,-3,3)){
for(y in runif(k,-2.5,2.5)){
if(abs(f(x,y))<tol) xy=rbind(xy,c(x,y))
}
}
nrow(xy)
plot(xy,pch=16,col=colours(),asp=1)
apply(xy,1,sum) |> summary()
apply(xy,1,prod) |> summary()