A matricial formalism to solve multi-dimensional initial boundary values problems for hyperbolic equations
written in quasi-linear based on the lambda scheme approach is presented. The derivation is carried out for nonorthogonal,
moving systems of curvilinear coordinates. A uniform treatment of the integration at the boundaries, when the boundary
conditions can be expressed in terms of combinations of time or space derivatives of the primitive variables, is also presented. The methodology is validated against a two-dimensional test case, the supercritical flow through the Hobson
cascade n.2, and in three-dimensional test cases such as the supersonic flow about a sphere and the flow through a plug
nozzle