Solving system nonlinear hyperbolic pdes

I have to solve a 1-D nonlinear pde system but I can't use pdepe because the equations are hyperbolic.

I think I can use solvepde (hence pde toolbox), even with coefficients c and a functions of the solution u, but I don't know how to assign the linear domain [0, L].

I have seen a single suggestion that recommends a rectangular domain in any case, with appropriate bc, but I do not know how to apply it in practice.

Otherwise, other solution possibilities in Matlab through already implemented functions? (without having to develop a finite-difference solution on the geometry)

Thank you so much in advance for the help, Giuseppe

[Bill Greene]

I suggest posting the definition of your problem-- PDEs, boundary conditions, and initial conditions.