A Gaussian quadrature method for the numerical solution of the characteristic intitial value problem _