No Title   |POLVER05_0
d(CT1)/d(t) = (beta*CT2-(1+beta)*CT1)/alpha/tau #
d(CT2)/d(t) = (beta*CT1-beta*CT2)/(1-alpha)/tau #
beta = 0.1 #
alpha = 0.8 #
tau = 40 #
CTe1 = 2000 - 59.6 * t + 0.64 * t ^ 2 - 0.00146 * t ^ 3 - 1.047 * 10 ^ (-5) * t ^ 4 #
CTe2 = 921 - 17.3 * t + 0.129 * t ^ 2 - 0.000438 * t ^ 3 + 5.6 * 10 ^ (-7) * t ^ 4 #
t1 = t-80 #
Cte = if(t<80)then(CTe1)else(CTe2) #
t(0)=0
CT1(0)=2000
CT2(0)=0
t(f)=160