  Multicycle nova evolution models have been calculated over the past twenty
years, the number being limited by numerical constraints. Here we present a
long-term evolution code that enables a continuous calculation through an
unlimited number of nova cycles for an unlimited evolution time, even up to (or
exceeding) a Hubble time. Starting with two sets of the three independent nova
parameters -- the white dwarf mass, the temperature of its isothermal core, and
the rate of mass transfer on to it -- we have followed the evolution of two
models, with initial masses of 1 and 0.65 solar masses, accretion rates
(constant throughout each calculation) of 1e-11 and 1e-9 solar-masses/yr, and
relatively high initial temperatures (as they are likely to be at the onset of
the outburst phase), through over 1000 and over 3000 cycles, respectively. The
results show that although on the short-term consecutive outbursts are almost
identical, on the long-term scale the characteristics change. This is mainly
due to the changing core temperature, which decreases very similarly to that of
a cooling white dwarf for a time, but at a slower rate thereafter. As the white
dwarf's mass continually decreases, since both models lose more mass than they
accrete, the central pressure decreases accordingly. The outbursts on the
massive white dwarf change gradually from fast to moderately fast, and the
other characteristics (velocity, abundance ratios, isotopic ratios) change,
too. Very slowly, a steady state is reached, where all characteristics, both in
quiescence and in outburst, remain almost constant. For the less massive white
dwarf accreting at a high rate, outbursts are similar throughout the evolution.
