Calculate Pi with Monte-Carlo-method

From JSXGraph Wiki
Revision as of 15:21, 29 October 2009 by A WASSERMANN (talk | contribs)
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Here, at construction time each point receives a function pair as coordinates. In each update these functions which return Math.random() are called. Thus in each update each point receives new random coordinates. The 50 points are updated on the onmousemove event.

This can be used to calculate [math]\displaystyle{ \pi }[/math] using statistics. After each update we count the number of points inside the circle (with midpoint [0,0] and radius 1) and outside the circle. The points are restricted to the square with midpoint [0,0] and edges of length 2.

Please, move the mouse pointer over this area:

<div id="jxgbox" class="jxgbox" style="width:400px; height:400px;" onmousemove="board2.update()"></div>
 board2 = JXG.JSXGraph.initBoard('jxgbox', {boundingbox: [-1, 1, 1, -1], axis:true});
 c = board2.createElement('circle', [[0,0], 1]);
 board2.suspendUpdate();
 var p2 = [];
 for (var i=0;i<50;i++) {
   p2[i] = board2.createElement('point',
           [function(){return 2*Math.random()-1;},function(){ return 2*Math.random()-1;}],
           {style:5,name:' '});
 }
 board2.unsuspendUpdate();

 var tin = 0;
 var tout = 0;
 var count = 0;

 function updText() {
   count++;
   var text = '';

   var inp = 0;
   var outp = 0;

   for(var i=0; i<p2.length;i++) {
     var x = p2[i].X();
     var y = p2[i].Y();

     if(x*x+y*y <= 1)
       inp++;
     else
       outp++;
   }

   tin += inp;
   tout += outp;

   text += '<b><u>Current</u></b><br/>In: ' + inp + ', out: ' + outp + ', total: ' + (inp+outp) + '; ratio: ' + (inp/(inp+outp)) + ', ratio*4: ' + (4*inp/(inp+outp)) + '.<br /><b><u>Total</u></b> (' + count + ' updates in total)<br/>In: ' + tin + ', out: ' + tout + ', total: ' + (tin+tout) + '; ratio: ' + (tin/(tin+tout)) + ', ratio*4: ' + (4*tin/(tin+tout));
   document.getElementById('resulttext').innerHTML = text;
 }

 var hookid = board2.addHook(updText);