Das deutsche Python-Forum

Code: Alles auswählen

import math
import random

def intp(x):
    if x>= 1.0:
        return 1
    if x<0.0:
        return -1
    return 0


t=int(input("Schusszahl eingeben:"))
s=0
i=1
while i <=t:
    x =  random.random()
    y = random.random()
    a=pi*random.random()
    b=math.cos(a)/2
    c=math.sin(a)/2
    if intp(y-b)!= intp(y+b):
        return s+1
    if intp(x-c)!= intp(x+c):
        return s+1
    
print('PI: '+str((t/s)*4))

Code: Alles auswählen

import math
import cmath                            # mathematische Konstanten
import random

def intp(x):
    if x>= 1.0:
        return 1
    if x<0.0:
        return -1
    return 0


t=int(input("Schusszahl eingeben:"))
s=0
i=1
while i <=t:
    x = random.random()                 # zufällige Koordinaten der Nadel
    y = random.random()                 #          --""--
    a=cmath.pi*random.random()          # Winkel der Nadel von 0 bis PI
    b=math.cos(a)/2
    c=math.sin(a)/2
    if intp(y-b)!= intp(y+b):
        s=s+1
    if intp(x-c)!= intp(x+c):
        s=s+1
    i=i+1

    print('PI: '+str((t/s)*4.0))
   

Code: Alles auswählen

#!/usr/bin/python
# -*- coding: utf-8 -*-
import math
import random
import sys

# import psyco
# psyco.full()
# ^ Damit ginge es noch schneller

intp = lambda x: x >= 1.0 and 1 or x < 0.0 and -1 or 0
# Warum nicht auch mit einer lambda Anweisung? 
# Du kannst aber natürlich auch deine intp Funktion weiter verwenden, das war nicht das Problem.

def buffon(t):
    s = 0
    for i in xrange(t):
        x = random.random()
        y = random.random()
        a = math.pi * random.random()
        b = math.cos(a) / 2.0
        c = math.sin(a) / 2.0
        if intp(y - b) != intp(y + b):
            s = s + 1
        if intp(x - c) != intp(x + c):
            s = s + 1
    return s


# Damit man Zeitunterschiede der Implementierung messen kann,
# soll das Programm auch ein Argument annehmen koennen.
if len(sys.argv) > 1:
    t = int(sys.argv[1])
else:
    t = int(input('Schusszahl eingeben: '))

s = buffon(t)
print 'PI: ' + str(4.0 * t / s)

Code: Alles auswählen

from __future__ import division

Code: Alles auswählen

from itertools import islice
from math import cos, pi as PI, sin
from random import random


BELOW, WITHIN, ABOVE = 'BELOW', 'WITHIN', 'ABOVE'
HORIZONTAL, VERTICAL = ORIENTATIONS = 'HORIZONTAL', 'VERTICAL'


def check_interval(x, lower=0, upper=1):
    if x < lower:
        return BELOW
    elif x >= upper:
        return ABOVE
    else:
        return WITHIN


class UnitNeedle(object):
    def __init__(self, x, y, angle):
        self.x = x
        self.y = y
        self.angle = angle
    
    @classmethod
    def random(cls):
        return cls(random(), random(), PI * random())

    def crosses_line(self, orientation):
        if orientation not in ORIENTATIONS:
            raise ValueError('orientation not in %r' % ORIENTATIONS)
        extent = (sin if orientation is HORIZONTAL else cos)(self.angle) / 2
        value = self.x if orientation is HORIZONTAL else self.y
        return check_interval(value - extent) != check_interval(value + extent)
    
    @property
    def line_cross_count(self):
        return sum(map(self.crosses_line, ORIENTATIONS))


def main():
    times = int(raw_input('times = '))
    needles = islice(iter(UnitNeedle.random, None), times)
    cross_count = sum(needle.line_cross_count for needle in needles)
    print 'PiDach =', 4.0 * times / cross_count


if __name__ == '__main__':
    main()

Code: Alles auswählen

intp=lambda n:n>1 or-(n<0)or 0

Code: Alles auswählen

>>> 0x7f
127
>>> 0b01111
15
>>> 0o664
436

Code: Alles auswählen

import math as m,random as r