Integral from (0,0) to (x,y)

I wanted to implement some of the integral solutions shown here,
in particular the non-homogeneous Dirichlet integral,

https://en.wikipedia.org/wiki/Heat_equation

Is there some way to do that with int2d?
I guess I could make a matrix to do it but it sounds involved.

Thanks.