I need to solve the following set of ODE’s numerically,

$ \frac{dy}{dx}=Ay+Bz, \ \frac{dz}{dx}=Cy+Dz$

where the independent variable $ x$ is an * non-uniform spaced* array of points and so is A, B, C & D.

What is the best way of achieving this as common NDsolve specifies the independent variable in $ \{x,x_{min},x_{max}\} $ format.