powered by
dr = alpha(beta-r)dt + sigma dW
with market price of risk q(r) = q1 +q2 r. The time scale is in years and the units are percentages.
LogLikVasicek( theta, R, tau, days, n)