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

564 ：卵の名無しさん ：2024/04/04(木) 09:54:04.69 ID:WjlV2HHp0.net

# https://ja.wikipedia.org/wiki/パスカルの蝸牛形
# 実数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()