Computing eigenpairs of two-parameter Sturm-Liouville systems using the bivariate sinc-Gauss formula