### How to solve normal distribution equation in Sage

I'm practicing statistics and I'm wondered how one can solve the following problem from Sage:

Find a numerical value of $x$ such that $\frac{1}{\sqrt{2\pi}}\int_{-\infty}^x e^{-t^2/2}dt=0.987654321$. I was thinking to different solutions:

1) Is there an inverse of normal distribution cumulative function?

2) Can we write it as $f'(x)=e^{-x^2/2},f(0)=1/2,f(x_0)=0.987654321$ and use some numerical method to solve the differential equation?

but I don't know are those functions implemented in Sage.