P1 = [4.0554489230596245, 31.525109459989388, 57.16281922464213, 
    44.08050738932008, 14.684956192885803, 2.1866330685079025, 
    -0.1402560791713545, -0.03504246268278482, -0.0008574567851546854]
Q1 = [1, 15.779988325646675, 45.39076351288792, 41.3172038254672, 
    15.04253856929075, 2.504649462083094, -0.14218292285478779, 
    -0.03808064076915783, -0.0009332594808954574]
P2 = [3.2377489177694603, 6.915228890689842, 3.9388102529247444, 
    1.3330346081580755, 0.20148538954917908, 0.012371663481782003, 
    0.00030158155350823543, 2.6580697468673755e-06, 6.239745391849833e-09]
Q2 = [1, 6.02427039364742, 3.6798356385616087, 1.3770209948908132, 
    0.21623699359449663, 0.013420400608854318, 0.00032801446468212774, 
    2.8924786474538068e-06, 6.790194080099813e-09]

def ndtri(y0):
    pass
How to calculate a Normal Distribution percent point function in python