function y=fvalue_ini2(h, c, alpha, theta, k, T)
y=T-alpha*c/((1-theta)*(k/h)^theta)-h;