3In the original paper [8], equation (154) was written with cos 2θ instead of (2 cos θ)2, but all the results presented in that paper were computed with OOFEM using an implementation based on (154).