Integrate a bivariate normal distribution over y when x is fixed

I have a bivariate normal distribution f(x,y). I want to integrate f(x=a, y) over y, say from y1 to y2. Is there any direct/simpler way to do this without using the integral function?

Thank you.

Answer 1

[Paul]

One simple approach is to evaluate f(a,y) for vector of y from y1 to y2 and then use cumtrapz() or trapz() on the result.