fixed again pb with constant functions. WIP
still fixing Psi & Psi. WIP
still fixing Psi & Phi: applied to all subsections corrections
WIP
Performed speed of Psi/Phi by using tabulated functions