# L'intersezione di  due figure
par( mai = c(0.5,0.5,0.2,0.2) )
x1 <- -3; x2 <- 3; y1 <- -2; y2 <- 4
plot(c(x1,x2),c(y1,y2),type="n",xlab="", ylab="",asp=1)
abline(h=seq(y1,y2,1),v=seq(x1,x2,1),col="blue",lty=3)
abline(h=0,v=0)
symbols(-1, 1, circles=1.5, inches=FALSE, add=TRUE, fg="blue")
f <- function(x) -x+1/2; plot(f,x1,x2,add=TRUE,col="blue")
# genero 20 mila punti nel rettangolo cartesiano e traccio quelli
# che stanno nel cerchio E sotto alla retta:
np <- 2e4; po <- x1+runif(np)*(x2-x1); pv <- y1+runif(np)*(y2-y1)
for(i in 1:np) if( f(po[i]) < pv[i] & (po[i]+1)^2+(pv[i]-1)^2 < 1.5^2)
                  points(po[i],pv[i],pch=".",col="red")
   
# Ovvero:
x1 <- -3; x2 <- 3; y1 <- -2; y2 <- 4
plot(c(x1,x2),c(y1,y2),type="n",xlab="", ylab="",asp=1)
abline(h=seq(y1,y2,1),v=seq(x1,x2,1),col="blue",lty=3)
abline(h=0,v=0)
symbols(-1, 1, circles=1.5, inches=FALSE, add=TRUE, fg="blue")
plot(f,x1,x2,add=TRUE,col="blue")
np <- 2e4; po <- x1+runif(np)*(x2-x1); pv <- y1+runif(np)*(y2-y1)
for(i in 1:np) if( f(po[i]) < pv[i]) points(po[i],pv[i],pch=".",col="red");
np <- 2e4; po <- x1+runif(np)*(x2-x1); pv <- y1+runif(np)*(y2-y1)
for(i in 1:np) if( (po[i]+1)^2+(pv[i]-1)^2 < 1.5^2) points(po[i],pv[i],pch=".",col="green4")