``````

In [17]:

import numpy as np
import scipy
import matplotlib.pyplot as plt

c = 299792458
h = 6.62607004e-34
hr = h / (2 * np.pi)
m_e = 9.10938356e-31
m_p = 1.6726219e-27
e = 1.60217662e-19
e0 = 8.85418782e-12
pi = np.pi
R = m_e * e ** 4 / ((4 * pi * e0) ** 2 * 4 * pi * hr ** 3 * c)

``````
``````

In [61]:

#exercicio 2
D = (2 * 79 * e ** 2) / (4 * pi * e0 * (7.7 * e * 1e4))
b = D / (2 * np.tan(pi / 3))
print("D(m): %e, b(SI):%2e" % (D, b))

``````
``````

D(m):2.954732e-12, b(SI):8.529578e-13

``````
``````

In [72]:

#exercicio 3
e_k = m_e * e ** 4 / (8 * e0 ** 2 * h ** 2)
delta_e = e_k * (1 - 1/16)
delta_v = np.sqrt(delta_e * 2 / m_p)
print("dE(eV): %f, dV(m/s): %e" % (delta_e / e, delta_v))

``````
``````

dE(eV): 12.755337, dV(m/s): 4.943303e+04

``````
``````

In [81]:

#exercicio 4

np.sqrt(63.3e9 / (R * 0.75))

``````
``````

Out[81]:

87.698882117932513

``````
``````

In [40]:

#exercicio 6
atom_sizes = np.array([1, 2, 3, 4, 6, 7])
red_mass = atom_sizes / (atom_sizes + m_e)
w = 0.75 / (R * red_mass * atom_sizes ** 2)

``````
``````

In [48]:

w[0]

``````
``````

Out[48]:

6.8345028952029222e-08

``````
``````

In [49]:

#exercicio 7
R/c * 0.75

``````
``````

Out[49]:

0.027453321247497935

``````
``````

In [25]:

c / 2.3e-47

``````
``````

Out[25]:

1.3034454695652174e+55

``````
``````

In [ ]:

``````
``````

In [ ]:

# exercicio 7

``````