def advance_cov(U, C, x, dx):
    """
    Explicit integration of the covariance term, using scipy's sparse matrix for C
    U - density function
    C - transition matrix
    """
    U = (C * U.reshape(len(x) ** 2)).reshape(len(x), len(x))
    return U
