Advertisement
Guest User

Untitled

a guest
Nov 30th, 2015
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.05 KB | None | 0 0
  1. from numpy import random
  2. from math import sqrt
  3.  
  4. L=input('Enter L: ') #in cm
  5.  
  6. NumberOfFissions=1
  7. PositionOnLine=L*random.random()
  8. print PositionOnLine
  9. R=sqrt(71.4) #in cm
  10.  
  11. #number of fissions produced in fission is 2
  12. #distance to next fission R=sqrt(71.4)
  13.  
  14. for i in range(1,101):
  15. if NumberOfFissions!=0:
  16. for e in range(1,NumberOfFissions+1):
  17. Decider=random.random()
  18. if Decider>0.5: #moves neutron right
  19. PositionOnLine+=R #if in 0 to L, fission occurs
  20. if PositionOnLine<=L or PositionOnLine>=0:
  21. NumberOfFissions+=1
  22. if PositionOnLine>L or PositionOnLine<0: #if not in 0 to L, neutron escapes
  23. NumberOfFissions-=1
  24. if Decider<=0.5: #moves neutron left
  25. PositionOnLine-=R
  26. if PositionOnLine<=L or PositionOnLine>=0:
  27. NumberOfFissions+=1
  28. if PositionOnLine>L or PositionOnLine<0:
  29. NumberOfFissions-=1
  30. print NumberOfFissions
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement