FormFactorPrism6 now computed with new analytical expression,

based on contour integral, without apparent singularities.
Validated against old (Rémy/Céline) implementation: agrees,
except for strong deviations near q->0, where the new
implementation is smoother.
The analytical calculations are explained in App. C2 of the manual.
This resolves #1300.