A = A((div(1/lambda)+1) div(u phi)) for E = omega A sine (omega t)