from java.util import Random
from java.lang import Math
from jhplot import *

h1 = H2D("PI",100,0.0,1.0,100,0.0,1.0)

nThrows=0
nSuccess=0
for i in range(1000000):
 x,y=Math.random(),Math.random()
 nThrows=nThrows+1
 if (x*x+y*y<=1):
              nSuccess=nSuccess+1;
              h1.fill(x,y)

c1 = HPlot3D("Canvas",600,600)
c1.visible()
c1.setNameX("X")
c1.setNameY("Y")
c1.setRangeZ(0,20)
c1.setDensity()
c1.draw(h1)

guesssed=4*nSuccess/float(nThrows)
print "Pi/4=",nSuccess/float(nThrows)
print "Pi=",guesssed
print "Actual Pi=",Math.PI
print "Difference=",(1-guesssed / Math.PI)*100