return 0.0, mu2, 0.0, (a+2.0)*(a+3.0)/mu2-3.0