An internal function to compute m and gamma based on continuous-time model (Yashin et. al., 2007)
func1(tt, y, a, f1, Q, f, b, theta)
list(m, gamma) Next values of m and gamma (see Yashin et. al, 2007)
tt - time
y
a (see Yashin et. al, 2007)
f1 (see Yashin et. al, 2007)
Q (see Yashin et. al, 2007)
f (see Yashin et. al, 2007)
b (see Yashin et. al, 2007)
theta