fixed some bugs, still to check with the Patrick Ciarlet known analytic solution
fixed several bugs. TODO: instead of solely Bc, Ec in the expressions for E3, B3, use Ec+Eg, Bc+Bg, this must give perfect results