In the VCO model given in the paper replace the last line of the model,
Code:V(out) <+ transition(n ? Vhi : Vlo, 0, tt);
with one that computes a sine wave based on the phase variable already present in the model.
Code:V(out) <+ ampl*cos(phase);
-Ken