c'est pas cela que tu cherches à faire ?
pi=vector()
for(n in 1:1000){
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
pi[n]=4*p
}
plot(pi,type="l")
c'est pas cela que tu cherches à faire ?
pi=vector()
for(n in 1:1000){
XY=matrix(runif(2*n),ncol=2)
p=mean(XY[,1]^2+XY[,2]^2<1)
pi[n]=4*p
}
plot(pi,type="l")
Vous avez un bloqueur de publicités installé.
Le Club Developpez.com n'affiche que des publicités IT, discrètes et non intrusives.
Afin que nous puissions continuer à vous fournir gratuitement du contenu de qualité, merci de nous soutenir en désactivant votre bloqueur de publicités sur Developpez.com.