The aim of this work is, commencing from an arbitrary starting point,
to locate regions near minima on the molecular ab initio potential ene
rgy surface. We work within the self-consistent field approximation. W
e use the method of dynamic evolution to optimise the geometry and the
molecular wavefunction, which is itself optimised at every nuclear st
ep. Examples are given for small and medium sized molecules.