Elipse

R Markdown

library(ggplot2)
set.seed(666)
va_x <- runif(5000, 2, 10)
va_y <- runif(5000, 2, 10)

library(plotrix)
set.seed(123)
x = runif(5000, 2, 10)
y = runif(5000, 2, 10)
plot(2:10, 2:10,
     type = "n",
     xlab = "eje x",
     ylab = "eje y") # Generamos un plot sin gráficos internos (limpio)
points(x,y,cex=0.01)
draw.ellipse(6,6,2,1, border = 'blue') # Dibuja la elipse
points_quarter = 0

for (i in 1:5000){  #tener en cuenta la ecuacion canonica
  if ((((x[i]-6)^2)/4)+ (((y[i]-6)^2)/1)<=1 ){
    points_quarter = points_quarter + 1
    points(x[i], y[i], col= "red", pch=16, cex=0.3) # Condicional que dibuja 1/4 de la elipse
  }
}

area_calc = pi * 1 * 2; area_calc

## [1] 6.283185

area_aprox = (points_quarter)/75; area_aprox

## [1] 6.613333