# System of Differential equations on non uniform grid

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.