Ajuste de funciones probabilisticas

Bien, trataré de der lo más claro posible.
Tengo una serie de valores medidos (ya tengo la frecuencia acumulada F(x)), de los cuales quiero ajustar funciones de distribución de probabilidad, (Gumbel, doble Gumbel, Normal, Log-Normal, etc), en este momento estoy atorado con la función de distribución Normal.
Al despejar la variable (considerando que tengo como dato F(x), y mi variable es x), queda una raíz con el radicando negativo, por lo cual no puedo sustituir los valores originales de F(x)
Ya busque si existe alguna forma para ajustarla por método numérico, pero nada.
Por favor si alguien me puede ilustrar se lo agradeceré.
Saludos