In this paper we introduce a new symbolic Gaussian formula for theevaluation of an integral over the first quadrant in a Cartesianplane, in particular with respect to the weight function$w(x)=\exp(-x^T x-1/x^T x)$, where $x=(x_1,x_2)^T\in \mathbb{R}^2_+$. Itintegrates exactly a class of homogeneous Laurent polynomials withcoefficients in the commutative field of rational functions in twovariables. It is derived using the connection between orthogonalpolynomials, two-point Padé approximants, and Gaussian cubatures.We also discuss the connection to two-point Padé-typeapproximants in order to establish symbolic cubature formulas ofinterpolatory type. Numerical examples are presented toillustrate the different formulas developed in the paper.