def inv_gaussian_variance(mu, lam):
    """
    As in https://en.wikipedia.org/wiki/Inverse_Gaussian_distribution
    @param mu: the expectation
    @param lam: lambda in Wikipedia's notation
    @return: Var[X]
    """
    return mu ** 3 / lam