def make_covmatrix(sigma_x, sigma_y, rho=0):
    """ create covariance matrix from given standard deviations and normalised correlation rho, vectorized """
    return (sigma_x ** 2, rho * sigma_x * sigma_y), (rho * sigma_x *
        sigma_y, sigma_y ** 2)
