import pylab import numpy as np #Reconstruccion A=1. a0=1./2 t=np.linspace(-3,3,10000) s_per=0 for k in range(-30,30): #ak=((-1)^k-1)*A/(np.pi^2*k^2) if k!=0: s_per=s_per+((-1)**k-1)*A/(np.pi**2*k**2)*np.exp(1j*k*2*np.pi*t) s_per=s_per+a0 pylab.plot(t,s_per), pylab.show() pylab.grid(True)
No hay comentarios.:
Publicar un comentario