Changes

*Neumann boundary condition: ∂ψ /∂n = k = const.

In the above, ∂ψ/∂n denotes the normal derivative of the potential at the surface of the domain boundary. [[EM.Ferma]]'s default domain boundary condition for both the electrostatic and magnetostatic solvers is Dirichlet. At the interface between different material media, additional boundary conditions must be applied. These boundary conditions involve electric or magnetic field components. The field components can be expressed as partial derivatives of the potential, i.e. in the form of ∂ψ/∂x, ∂ψ/∂y or ∂ψ/∂z. Using the equivalent finite difference approximations of these derivative, one arrive at fairly complicated difference equations involving the constitutive parameters ε, μ and σ.