import math


def norm_pdf(x, mu=0, sigma=1): # type: (float, float, float) -> float
    """Calculate a the probability density of a normal distribution"""
    return math.exp(-(((x - mu) ** 2) / (2 * (sigma ** 2)))) / math.sqrt(2 * math.pi * (sigma ** 2))