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.