In reply to Rick Aultman.

Use the QUAD function to integrate over y, where x is a parameter. If you have questions, post them to the SAS/IML Support Community along with definition of G(x,y) and the values of k.