Вычисление вероятности случайной переменной в распределении в Python

Учитывая среднее значение и стандартное отклонение, определяющее нормальное распределение , как бы вы вычислили следующие вероятности в чистом Python (т.е. без Numpy / Scipy или других пакетов, не входящих в стандартная библиотека)?

1. Вероятность случайной величины r, где r
2. Вероятность случайной величины r, где r> x или r> = x.
3. Вероятность того, что случайная величина r, где x> r> y.

Я нашел несколько библиотек, например Pgnumerics , которые предоставляют функции для их вычисления, но лежащая в основе математика мне непонятна.

Править : Чтобы показать, что это не домашнее задание, ниже приведен мой рабочий код для Python

from math import *
import unittest

def erfcc(x):
    """
    Complementary error function.
    """
    z = abs(x)
    t = 1. / (1. + 0.5*z)
    r = t * exp(-z*z-1.26551223+t*(1.00002368+t*(.37409196+
        t*(.09678418+t*(-.18628806+t*(.27886807+
        t*(-1.13520398+t*(1.48851587+t*(-.82215223+
        t*.17087277)))))))))
    if (x >= 0.):
        return r
    else:
        return 2. - r

def normcdf(x, mu, sigma):
    t = x-mu;
    y = 0.5*erfcc(-t/(sigma*sqrt(2.0)));
    if y>1.0:
        y = 1.0;
    return y

def normpdf(x, mu, sigma):
    u = (x-mu)/abs(sigma)
    y = (1/(sqrt(2*pi)*abs(sigma)))*exp(-u*u/2)
    return y

def normdist(x, mu, sigma, f):
    if f:
        y = normcdf(x,mu,sigma)
    else:
        y = normpdf(x,mu,sigma)
    return y

def normrange(x1, x2, mu, sigma, f=True):
    """
    Calculates probability of random variable falling between two points.
    """
    p1 = normdist(x1, mu, sigma, f)
    p2 = normdist(x2, mu, sigma, f)
    return abs(p1-p2)