[sfepy] How to calculate electrostatic field using sfepy?