relacion de dispersion del TE1

 cell = mp.Vector3(16, 8, 0)
 geometry = [mp.Block(mp.Vector3(mp.inf, 1, mp.inf),center=mp.Vector3(),material=mp.Medium(epsilon=4),)]
 sources = [mp.Source(mp.ContinuousSource(frequency=0.4), component=mp.Ez, center=mp.Vector3(-7, +0.25),amplitude=+1),
            mp.Source(mp.ContinuousSource(frequency=0.4), component=mp.Ez, center=mp.Vector3(-7, -0.25),amplitude=-1)]
 pml_layers = [mp.PML(1.0)]
 resolution = 10
 sim = mp.Simulation(cell_size=cell,boundary_layers=pml_layers,geometry=geometry,sources=sources,resolution=resolution,)
 sim.run(until=200)
 ez_data = sim.get_array(center=mp.Vector3(), size=cell, component=mp.Ez)

 plt.imshow(ez_data.transpose())
 plt.show()
 plt.plot(ez_data[:,35])   
 plt.show()

 for ic in range(NQx):
  Qx = Qx_i+ic*dQx
  QxV[ic]=Qx
  Exr=0
  Exi=0
  for i in range(I):
   Exr=Exr+ez_data[i_0+i,40]*math.cos(2.0*math.pi*Qx*((i_0+i)/resolution))
   Exi=Exi+ez_data[i_0+i,40]*math.sin(2.0*math.pi*Qx*((i_0+i)/resolution))
  EQx[ic] = math.sqrt(Exr*Exr+Exi*Exi)