In [ ]:

    

import numpy as np
import matplotlib.pyplot as plt


    

In [ ]:

    

N0, w, k = 0.5, 0.43, 0.014

x = np.arange(1, 1001)
Ne0 = 2 * N0 * np.cos(k * x) * np.cos(-w * 0)
Ne2 = 2 * N0 * np.cos(k * x) * np.cos(-w * 2)
Ne8 = 2 * N0 * np.cos(k * x) * np.cos(-w * 8)
Ne9 = 2 * N0 * np.cos(k * x) * np.cos(-w * 9)


    

In [ ]:

    

fig, ax = plt.subplots()
ax.plot(Ne0, 'k:')
ax.plot(Ne2, 'b')
ax.plot(Ne8, 'k')
ax.plot(Ne9, 'b:')
ax.set_title(u'Onda estacionária')
ax.set_xlabel('Eixo x [m]')
_ = ax.set_ylabel(r'$\eta$ [m]')