omega = k c A = -i A0 / omega (exp(i k z - i omega t), 0, 0) E = -1/c d(A,t) B = curl(A) E B "Verify Maxwell equations" check(div(E) == 0) check(div(B) == 0) check(curl(E) == -1/c d(B,t)) check(curl(B) == 1/c d(E,t)) "ok" Run