r 第一次存入Rt代码：lesoleilàl'horizo​​n avec une inconsistance dans le traitement cos / sin qui do

Posted 2021-05-24

x = runif(10000, 0, 99)
y = runif(10000, 0, 99)
plan = data.frame(x, y)
plan$dist2=sqrt((plan$x-60)^2 + (plan$y-60)^2)
plan$ang2 = ifelse(plan$y-60>0, acos(abs(plan$x-60)/plan$dist2), asin(abs(plan$x-60)/plan$dist2))
plot(plan$x, plan$y, xlab="", ylab="", xlim=c(0,100), ylim=c(0,100), pch=20, col="darkblue", cex=0.5)
segments(plan$x, plan$y, plan$x+(1 * cos(plan$ang2)), plan$y+(1 * sin(plan$ang2)), col=ifelse(plan$dist2>15, "black", "red"), lwd=2 )