import math
#Variable declaration
rho=9*10**9#in nt-m2/coul2
e=1.6*10**-19#coul
r=1*10**-10#in m
#Calculations
k=(rho*e**2)/(r**3)#nt/m
m=9.11*10**-31#in kg
c=3*10**8#in m/s
v=(1/(2*math.pi))*math.sqrt(k/m)
lamda=c/v
#Result
print "Part(b)\nThe wavelength is %.1e m "%lamda
import math
#Variable declaration
N=10**4#in rad, Number of atoms tarversed
#Calculations
theta=(2*10**-2)/math.sqrt(N)
#Result
print "Average deflection %.1e rad "%theta
#Variable declaration
rho=9*10**9#in nt-m2/coul2
m=9.11*10**-31#in kg
e=1.6*10**-19#coul
h=1.05*10**-34#in j-sec
#Calculations
E=-(rho*m*e**4)/(2*h**2)
#Result
print "Binding energy is %.2e Joule "%E
#Answer given in the book is wrong
#Variable declaration
mu=207#207*me
#Calculations&Results
M=1836#183*me
u=(mu*M)/(mu+M)
D=(1./u)*5.3*10**-11
print "Part(a)\nMuon nucleus seperation is %.1e m "%D
E=-u*13.6
print "\nPart(b)\n Binding energy is %.f ev "%E
R=109737#in cm
lamda=(1./u)*(1/0.75)*(1./R)
print "\nPart(c)\n Wavelength is %.1e cm "%lamda
#Calculations&Results
#For Hydrogen atom
R_o=109737#in cm
m=1
M=1836.
RH=(R_o)/(1+(m/M))
print "Spectrum line for Hydrogen occur at %.f /cm "%RH
#For Deuterium atom
R_o=109737#in cm
m=1
M=2*1836.
RD=(R_o)/(1+(m/M))
print "Spectrum line for Deuterium occur at %.f /cm "%RD