Puesto que estamos interesados en encontrar el coste unitario, dejemos que la producción necesaria sea igual a 1, es decir, $y(i)\overset{set}=1$ Tenemos que resolver: $$\begin{align} \min_{l(i),x(i),h(i)\geq0} \quad & w_Ll(i)+x(i)+w_{H}h(i)\\ \textrm{s.t.} \quad & [l(i)^{\frac{\epsilon - 1}{\epsilon}} +\alpha(i)(\tilde{\gamma}x(i))^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon \beta}{\epsilon-1}}h(i)^{1-\beta}=1\\ \textrm{and} \quad & \alpha(i)=\begin{cases}1 & \text{if } i\text{ is automated}\\ 0 & \text{otherwise}\end{cases} \end{align}$$
Voy a utilizar la función de indicador y resolver dos casos diferentes y luego tratar de combinarlos utilizando el indicador.
Caso 1: Si $i$ no está automatizado, es decir $ \alpha(i)=0$
Si $i$ no está automatizado, entonces el problema de minimización de costes se convierte en uno estándar con una función de producción Cobb-Douglas $$\begin{align} \min_{l,x,h\geq 0} \quad & w_Ll+x+w_{H}h\\ \textrm{s.t.} \quad & l^{\beta}h^{1-\beta}=1\\ \end{align}$$
Esto no es más que un problema estándar de minimización de costes con una función de producción cobb-douglas, así que creo que puedes comprobar tú mismo que el problema anterior da: $l(i)=\left(\frac{w_H\beta}{w_L(1-\beta)}\right)^{1-\beta}, \quad x(i)=0, \quad h(i)= \left(\frac{w_L(1-\beta)}{w_{H} \beta}\right)^{\beta}$
por lo tanto, $$\begin{eqnarray} & c=w_L^\beta\left(\frac{w_H\beta}{1-\beta}\right)^{1-\beta}+w_H^{1-\beta}\left(\frac{w_L(1-\beta)}{\beta}\right)^{\beta}\\\\ \implies & {\boxed{c=w_L^\beta w_H^{1-\beta} \beta^{-\beta}(1-\beta)^{\beta-1}}}\tag{a} \end{eqnarray}$$
Caso 2: Si $i$ está automatizado $$\begin{align} \min_{l,x,h\geq 0} \quad & w_Ll+x+w_{H}h\\ \textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon \beta}{\epsilon-1}}h^{1-\beta}=1\\ \end{align}$$
escribamos el problema anterior como $$\begin{align} {\min_{0\leq h \leq 1} \underbrace{\begin{pmatrix} \underset{l,x}{\min} \quad & w_Ll+w_Hh+x \\ \textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon}{\epsilon-1}}=h^\frac{\beta -1}{\beta} \end{pmatrix}}_{\text{auxiliary problem}}} \tag{1} \end{align}$$
Podemos resolver el problema anterior en dos etapas:
Fase 1: resolver primero el problema auxiliar para $l$ y $x$ en $h$ fijo
Fase 2: Sustituir las soluciones del problema auxiliar $-$ $l(h)$ y $x(h)$ $-$ en $(1)$ y resolver el problema para encontrar el óptimo $h$
Fase 1: que tenemos que resolver: $$\begin{align} \underset{l,x}{\min} \quad & w_Ll+w_Hh+x \\ \textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon}{\epsilon-1}}=h^\frac{\beta -1}{\beta}\\ \text{gives:} \quad & {\left.\begin{matrix} l(h)=\frac{h^\frac{\beta -1}{\beta}}{[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{\epsilon}{\epsilon-1}} \\ x(h)=\frac{w_L^\epsilon \tilde{\gamma}^{\epsilon -1} h^\frac{\beta -1}{\beta}}{[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{\epsilon}{\epsilon-1}} \end{matrix}\right\}}\tag{2} \end{align}$$
Porque se nos da que $\epsilon>1$ Se puede comprobar que las isocuantas serán convexas. Así, pude resolver lo anterior igualando las razones de MP a MC para $l$ y $x$ y sustituyéndolos en la restricción.
Fase 2: sustituyendo $(2)$ en $(1)$ obtenemos la resolución final del problema que dará como resultado óptimo $h$ $$\begin{align} \min_{0\leq h\leq 1} \quad & w_Hh+w_L[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{1}{1-\epsilon} h^\frac{\beta -1}{\beta}\tag{3} \end{align}$$
el objetivo anterior es convexo en $h \; \because \beta<1$ por lo que asumiendo la condición paramétrica obtenemos un punto estacionario que se encuentra en el intervalo $(0,1)$ podemos resolver la h óptima en el problema anterior y sustituirla en el propio problema para obtener $c$ para $\alpha=1$ :
$$h^*=\left(\frac{(1-\beta)w_L[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{1}{1-\epsilon}}{w_H\beta}\right)^\beta$$
sustituyendo $h$ en $(3)$ danos: $$\boxed{c=w_H^{1-\beta}(1-\beta)^{\beta-1}\beta^{-\beta}\left[w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1}\right]^\frac{\beta}{1-\epsilon}}\tag{b}$$
podemos combinar las soluciones de caso 1 y caso 2 es decir, $(a)$ y $(b)$ Por lo tanto, tenemos
$c=\begin{cases} w_H^{1-\beta} \beta^{-\beta}(1-\beta)^{\beta-1}w_L^\beta & \text{if } \alpha(i)=0 \\ w_H^{1-\beta}(1-\beta)^{\beta-1}\beta^{-\beta}\left[w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1}\right]^\frac{\beta}{1-\epsilon} & \text{if } \alpha(i)=1 \end{cases}$
utilizando $\alpha(i)$ lo anterior también se puede escribir como: $$\boxed{c=w_H^{1-\beta}\beta^{-\beta}(1-\beta)^{\beta-1}\left(w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1} \alpha(i) \right)^\frac{\beta}{1-\epsilon}}$$