Monty PI Notebook

{
 "metadata": {
  "name": "Untitled0"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "import random",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 3
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "def distance_to_origin(x,y):\n    return (x**2 + y**2)**.5",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 62
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "def run_exp(n):\n    in_cir, exp = 0, 0\n    for each in xrange(int(n)):\n        distance = distance_to_origin(random.random(), random.random())\n        in_cir += 1 if distance <= 1 else 0\n        exp += 1\n    return float(4*in_cir)/exp\n",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 63
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 67,
       "text": "[4.0, 4.0, 0.0, 4.0, 0.0, 4.0, 4.0, 4.0, 4.0, 4.0]"
      }
     ],
     "prompt_number": 67
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E2) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 68,
       "text": "[3.28, 3.24, 2.92, 3.32, 3.24, 3.56, 3.28, 2.96, 3.2, 3.28]"
      }
     ],
     "prompt_number": 68
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E3) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 69,
       "text": "[3.216, 3.068, 3.156, 3.168, 3.08, 3.192, 3.16, 3.244, 3.236, 3.204]"
      }
     ],
     "prompt_number": 69
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E4) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 70,
       "text": "[3.1536,\n 3.1224,\n 3.1472,\n 3.1176,\n 3.1444,\n 3.1404,\n 3.1204,\n 3.1732,\n 3.1408,\n 3.1392]"
      }
     ],
     "prompt_number": 70
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E5) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 71,
       "text": "[3.13468,\n 3.13756,\n 3.15212,\n 3.14168,\n 3.1422,\n 3.14376,\n 3.13936,\n 3.14624,\n 3.14428,\n 3.1398]"
      }
     ],
     "prompt_number": 71
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E6) for _ in xrange(10)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 72,
       "text": "[3.142852,\n 3.14042,\n 3.141908,\n 3.14438,\n 3.142544,\n 3.1411,\n 3.14238,\n 3.145932,\n 3.139852,\n 3.13852]"
      }
     ],
     "prompt_number": 72
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "[run_exp(1E7) for _ in xrange(3)] ",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 74,
       "text": "[3.1409556, 3.1426688, 3.1415524]"
      }
     ],
     "prompt_number": 74
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "",
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}