x_inf=670000
x_sup=672000
y_inf=9886000
y_sup=9888000
puntos_a_generar=5000
diametro=x_sup-x_inf
radio=diametro/2
origen_x=x_inf+radio
origen_y=y_inf+radio
 
datos=matrix(data=NA,nrow=0,ncol=3,byrow=FALSE,dimnames=NULL)
datos=as.data.frame(datos)
colnames(datos)=c("x","y","evento")
 
n = 0
while(n < puntos_a_generar) {
    x_ale=(x_inf+((seq(0,1,length=100)[sample(1:100,1)])*diametro))
    y_ale=(y_inf+((seq(0,1,length=100)[sample(1:100,1)])*diametro))
    if(sqrt((x_ale-origen_x)^2+(y_ale-origen_y)^2)<radio)
    {
        nuevafila=data.frame(x=c(x_ale),y=c(y_ale),evento=NA)
        datos=rbind(datos,nuevafila)
        n=length(datos$x)
    }
}
write.csv(datos,'datos.cvs')

eF9pbmY9NjcwMDAwCnhfc3VwPTY3MjAwMAp5X2luZj05ODg2MDAwCnlfc3VwPTk4ODgwMDAKcHVudG9zX2FfZ2VuZXJhcj01MDAwCmRpYW1ldHJvPXhfc3VwLXhfaW5mCnJhZGlvPWRpYW1ldHJvLzIKb3JpZ2VuX3g9eF9pbmYrcmFkaW8Kb3JpZ2VuX3k9eV9pbmYrcmFkaW8KCmRhdG9zPW1hdHJpeChkYXRhPU5BLG5yb3c9MCxuY29sPTMsYnlyb3c9RkFMU0UsZGltbmFtZXM9TlVMTCkKZGF0b3M9YXMuZGF0YS5mcmFtZShkYXRvcykKY29sbmFtZXMoZGF0b3MpPWMoIngiLCJ5IiwiZXZlbnRvIikKCm4gPSAwCndoaWxlKG4gPCBwdW50b3NfYV9nZW5lcmFyKSB7CiAgICB4X2FsZT0oeF9pbmYrKChzZXEoMCwxLGxlbmd0aD0xMDApW3NhbXBsZSgxOjEwMCwxKV0pKmRpYW1ldHJvKSkKICAgIHlfYWxlPSh5X2luZisoKHNlcSgwLDEsbGVuZ3RoPTEwMClbc2FtcGxlKDE6MTAwLDEpXSkqZGlhbWV0cm8pKQogICAgaWYoc3FydCgoeF9hbGUtb3JpZ2VuX3gpXjIrKHlfYWxlLW9yaWdlbl95KV4yKTxyYWRpbykKICAgIHsKICAgICAgICBudWV2YWZpbGE9ZGF0YS5mcmFtZSh4PWMoeF9hbGUpLHk9Yyh5X2FsZSksZXZlbnRvPU5BKQogICAgICAgIGRhdG9zPXJiaW5kKGRhdG9zLG51ZXZhZmlsYSkKICAgICAgICBuPWxlbmd0aChkYXRvcyR4KQogICAgfQp9CndyaXRlLmNzdihkYXRvcywnZGF0b3MuY3ZzJyk=