In reply to Haider Mannan.

I am unable to provide personal help, but I will point out that the integral of f over the interval [a,b] equals the integral on (-infinity, b) minus the integral from (-infinity, a). I will also mention that polyserial correlation is available in SAS by using PROC CORR. Good luck.