DESCRIPTION:

[LONGITUDINAL]
input ={foldx, Tmax, Cmax, alpha, FB, beta}

EQUATION:

rho=log(foldx)/Tmax
R0=Cmax/exp(rho*Tmax)

A = (1-FB)* Cmax
B=FB*Cmax

if t<0
	ylin=R0
elseif t<Tmax
	ylin=R0*exp(rho*t)
else
	ylin=A*exp(-alpha*(t-Tmax))+B*exp(-beta*(t-Tmax))
end
logy=log(ylin)


OUTPUT:
output = {logy}