RosettaCodeData/Task/Monte-Carlo-methods/PARI-GP/monte-carlo-methods.parigp

MonteCarloPi(tests)=4.*sum(i=1,tests,norml2([random(1.),random(1.)])<1)/tests;