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.