A Runge-Kutta Method for the Numerical Solution of the Goursat Problem in Hyperbolic Partial Differential Equations