Numerical solution of the Goursat problem on a triangular domain with mixed boundary conditions

For the Goursat problem, we consider a triangular domain with mixed Dirichlet and impedance boundary conditions imposed on it. We develop an algorithm for its numerical solution mainly based on Runge-Kutta method and trapezoidal formula. Iterative techniques are constructed to compute some data for the nonlinear part of the differential equation and the impedance boundary condition. Error estimates are derived. Examples are presented to illustrate the effectiveness of the method.