Statistics for Solving hyperbolic partial differential equations using multi-domain bivariate spectral collocation method