Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import math
- w = 0.0
- rez = 1024
- for v in range(0, rez):
- for u in range(0, rez):
- s = (2.0 * u)/rez - 1.0
- t = (2.0 * v)/rez - 1.0
- tmp = 1.0+s*s+t*t
- w += 4.0/(tmp * math.sqrt(tmp))
- print (4*math.pi/6), w/(rez*rez), (4*math.pi/6)-w/(rez*rez)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement