Advertisement
Guest User

ANDRISH DOLBOEB

a guest
Dec 13th, 2019
113
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 1.08 KB | None | 0 0
  1. #import math
  2. import random
  3. import numpy as np
  4. import matplotlib.pyplot as plt
  5. a = 0
  6. b = 20
  7. x = []
  8. for i in range(20):
  9.     x.append((20-0)*random.random()+0)
  10. print (x)    
  11. bins = [0,2,4,6,8,10,12,14,16,18,20]
  12. hist, bins = np.histogram(x, bins=bins)
  13. width = np.diff(bins)
  14. center = (bins[:-1]+bins[1:])/2
  15. fig, ax = plt.subplots(figsize=(8,3))
  16. ax.bar(center, hist, align='center',width=width)
  17. ax.set_xticks(bins)
  18. plt.show()
  19. sigma = 2
  20. mu = 0
  21. k = 0
  22. for i in range(100):
  23.     x = 0
  24.     for i in range(12):
  25.         x = x+random.random()
  26.     x = x - 6
  27.  
  28.     ksi = mu + sigma*x
  29.     if ((ksi>=0) and (ksi<=2)):
  30.         k += 1
  31. print (k/100)
  32. sigma = 2
  33. mu = 1
  34. k = 0
  35. for i in range(100):
  36.      x.append((100-0)*random.random()+0)
  37. print (x)
  38.     x = 0
  39.     for i in range(12):
  40.          x = x+random.random()
  41.     x = x - 6
  42.     teta = mu+ sigma*x
  43.     ksi=math.exp(teta)
  44.      if ((ksi>=-1) and (ksi<=1)):
  45.         k += 1
  46. print (k/100)
  47. l = 0.5
  48. k = 0
  49. for i in range(100):
  50.     a = random.random()
  51.     x = -math.log(a)/l
  52.     if ((x>=0) and (x<=1)):
  53.         k += 1
  54. print (k/100)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement