thetaf function

Theta method forecast