Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # -*- coding: cp1252 -*-
- #leibniz.py
- import math
- def leibniz(n):
- nenner=1
- piLeibniz=0.0
- for counter in range(0,n):
- piLeibniz+=((-1.0)**counter)*4/nenner
- nenner+=2
- return (piLeibniz)
- def wallis(abweichung):
- m=1
- zaehler=2.0
- nenner=3.0
- piWallis_=2.0
- while ((abs(2*piWallis_-pi))>10**-7+abweichung):
- if m%2==0:
- piWallis_*=float(zaehler/nenner)
- nenner+=2
- else:
- piWallis_*=float(zaehler/nenner)
- zaehler+=2
- m+=1
- return (m,piWallis_)
- pi=math.pi
- n=input("Bitte n eingeben: ")
- piLeibniz=leibniz(n)
- (m,piWallis_)=wallis(abs(piLeibniz-pi))
- print "n=",n,", m=",m, ", Abweichung Leinbniz=",abs(piLeibniz-pi), ", Abweichung Wallis=",abs(2*piWallis_-pi)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement