int(1/sigma^(N+M+1)*exp(-1/(2*sigma^2)* r^2))