Lax-Wendroff scheme for system of partial differential equations

I have a system of equations of the form

The problem is an initial-boundary value problem. Y=f(s,t) is an array of unknowns. 's' is the space dimension and 't' is the time. A and B are constant matrices and c is an array which is a function of the dependent variables contained in Y. Can I use Lax-Wendroff scheme to solve the PDE? If yes please help me to obtain the Lax-Wendroff scheme for the above set of equations, otherwise kindly suggest any other numerical method.