In [21]:

    

import numpy as np
import matplotlib.pyplot as plt_


    

In [31]:

    

G = 0.01
B_0 = 1


    

In [32]:

    

x = np.linspace(-0.02, 0.02, 100)
omega = 40*(B_0+G*x)


    

In [33]:

    

intensity = np.zeros_like(omega)


    

In [34]:

    

intensity[np.logical_and(x>0.005, x<0.015)] = 2
intensity[np.logical_and(x>-0.015, x<-0.005)] = 1


    

In [50]:

    

xticks = np.arange(-0.015, 0.016, 0.005)
yticks = np.array([0,1,2])
xticks_label = [str(i) for i in 40*(B_0+G*xticks)]
fig, axes = plt.subplots()
axes.plot(x, intensity)
axes.set_xlim(-0.02, 0.02)
axes.set_ylim(-0.5, 2.5)
axes.set_xlabel('Frquency [MHz]', fontsize='x-large')
axes.set_ylabel('Intensity', fontsize='x-large')
axes.set_xticks(xticks)
axes.set_yticks(yticks)
axes.set_xticklabels(xticks_label)


    

    
Out[50]:





[<matplotlib.text.Text at 0x10e71b9d0>,
 <matplotlib.text.Text at 0x10e715910>,
 <matplotlib.text.Text at 0x10e7414d0>,
 <matplotlib.text.Text at 0x10e741c10>,
 <matplotlib.text.Text at 0x10e745390>,
 <matplotlib.text.Text at 0x10e745ad0>,
 <matplotlib.text.Text at 0x10e748250>]

In [17]:

    

%matplotlib


    

    




Using matplotlib backend: MacOSX

In [12]:

    

np.logical_and(omega>1, omega<-1)


    

    
Out[12]:





array([False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False,
       False, False, False, False, False, False, False, False, False, False], dtype=bool)

In [46]:

    

xticks


    

    
Out[46]:





array([ -1.50000000e-02,  -1.00000000e-02,  -5.00000000e-03,
         3.46944695e-18,   5.00000000e-03,   1.00000000e-02,
         1.50000000e-02])

In [47]:

    

3.25*1.5


    

    
Out[47]:





4.875

In [49]:

    

1e-5*1.5/(4*np.pi*1e-7)


    

    
Out[49]:





11.936620731892152

In [ ]: