Calculating covariances of random functions linked by poisson equation in R2