In [50]:

    

%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt


    

In [51]:

    

m_e = 9.10938356e-31
eV = 1.60217662e-19
h = 1.0545718001e-34


    

In [55]:

    

#### 6 #####
L = 10e-10
E =  np.linspace(0+1e-15, 4, 1e6) * eV
U = 4 * eV
a = np.sqrt(2 * m_e * (U - E)) / h
b = np.sqrt(2 * m_e * E) / h
bla = ((2 * a * b) / (a ** 2 + b ** 2) * np.sinh(a * L)) ** 2

R = 1 / (bla + 1)
T = 1 - R
plt.plot(E / eV , R, ',-')
plt.show()


    

In [54]:

    

R[:10]


    

    
Out[54]:





array([  9.99999801e-01,   1.25778035e-03,   6.29293001e-04,
         4.19621405e-04,   3.14752611e-04,   2.51820773e-04,
         2.09861811e-04,   1.79888968e-04,   1.57408157e-04,
         1.39922382e-04])

In [ ]: