QCD at finite quark-/baryon-number density, which describes nuclear
matter, has a sign problem which prevents direct application of standard
simulation methods based on importance sampling. When such finite density is
implemented by the introduction of a quark-number chemical potential $\mu$, this
manifests itself as a complex fermion determinant. We apply simulations using
the Complex Langevin Equation (CLE) which can be applied in such cases. However,
this is not guaranteed to give correct results, so that extensive tests are
required. In addition, gauge cooling is required to prevent runaway behaviour.
We test these methods on 2-flavour lattice QCD at zero temperature on a small
($12^4$) lattice at an intermediate coupling $\beta=6/g^2=5.6$ and relatively
small quark mass $m=0.025$, over a range of $\mu$ values from $0$ to
saturation. While this appears to show the correct phase structure with a phase
transition at $\mu \approx m_N/3$ and a saturation density of $3$ at large
$\mu$, the observables show departures from known values at small $\mu$. We
are now running on a larger lattice ($16^4$) at weaker coupling $\beta=5.7$.
At $\mu=0$ this significantly improves agreement between measured observables
and known values, and there is some indication that this continues to small
$\mu$s. This leads one to hope that the CLE might produce correct results in
the weak-coupling -- continuum -- limit.