def METHOD_NAME(bF0, n, RadPHI0, PHI1):
    return float(bF0 * (((1 + n + ((5 / 4) * pow(n, 2)) + ((5 / 4) * pow(n, 3))) * (PHI1 - RadPHI0))
                  - (((3 * n) + (3 * pow(n, 2)) + ((21 / 8) * pow(n, 3))) * (sin(PHI1 - RadPHI0)) * (cos(PHI1 + RadPHI0)))
                  + ((((15 / 8) * pow(n, 2)) + ((15 / 8) * pow(n, 3))) * (sin(2 * (PHI1 - RadPHI0))) * (cos(2 * (PHI1 + RadPHI0))))
                  - (((35 / 24) * pow(n, 3)) * (sin(3 * (PHI1 - RadPHI0))) * (cos(3 * (PHI1 + RadPHI0))))))