integrate exp(x^2 + 2*x*y*C) dx for x=0..pi