  Ultracool dwarfs have emerged as key targets for searches of transiting
exoplanets. Precise estimates of the host parameters (including mass, age, and
radius) are fundamental to constrain the physical properties of orbiting
exoplanets. We have extended our evolutionary code CLES (Code Li\'egeois
d'Evolution Stellaire) to the ultracool dwarf regime. We include relevant
equations of state for H, He, as well as C and O elements to cover the
temperature-density regime of ultracool dwarf interiors. For various
metallicities, we couple the interior models to two sets of model atmospheres
as surface boundary conditions. We show that including C and O in the EOS has a
significant effect close the H-burning limit mass. The typical systematic error
associated with uncertainties in input physics in evolutionary models is $\sim
0.0005 M_\odot$. We test model results against observations for objects whose
parameters have been determined from independent techniques. We are able to
reproduce dynamical mass measurements of LSPM J1314+1320AB within $1\sigma$
with the condition of varying the metallicity (determined from calibrations) up
to $2.5\sigma$. For GJ 65AB, a $2\sigma$ agreement is obtained between
individual masses from differential astrometry and those from evolutionary
models. We provide tables of ultracool dwarf models for various masses and
metallicities that can be used as reference when estimating parameters for
ultracool objects.
