cosinor.lm.default function

Fit cosinor model