Second order pole fix started

In case of second order poles in the 4-point correlation function one needs to use different expressions.
This starts fixing this by distinguishing the cases, without having all the formulas implemented.