Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import random
- from matplotlib import pyplot as plt
- a=[]
- m=0
- k=0
- x,y = [],[]
- n,N=0,1000000
- #a[0] = x
- #a[1] = y
- for _ in range (N):
- for i in range (2):
- a.append(random.uniform(-7.80,7.80))
- if (a[0])**2 + (a[1])**2 <= 60:
- if (a[0])**2 + (a[1])**2 >=15:
- n+=1
- x.append(a[0])
- y.append(a[1])
- else:
- k+=1
- else:
- m+=1
- a=[]
- print(k)
- print(m)
- print(n)
- print(n+m+k)
- print((n+m)/N)
- plt.plot(x,y, linestyle="",marker="o")
- plt.title('PLOT')
- plt.show()
Advertisement
Add Comment
Please, Sign In to add comment