where, u=1/2(ε0E2+B2/μ0), is the energy density, and P→=(E→×B→)/μ0 is the poynting vector/energy flux
since B→=1ωk→×E→: