Non-separable utility in OLG
Here we analyse the behaviour of the savings function when utility is non-separable. In particular, let utility be represented by
$$U = u(c,d).$$
We use the substitution method, but in this case we will introduce different alternatives, depending on the derivative we compute. In particular, we introduce different elements of the budget constraint:
$$w = c + s$$ $$d = Rs$$ $$w = c + \frac{d}{R}$$
$c$ and $d$ are normal goods
First, we derive conditions such that $c$ and $d$ are normal goods. Optimality requires:
$$-u_c (w-s,Rs) + R u_d (w-s,Rs) =0.$$
Consumption when young
We are interested in computing $\frac{\partial c}{\partial w}$ and showing it is positive so that $c$ is a normal good. First, we rewrite the first order condition using some substitutions:
$$-u_c (c, R(w-c)) + R u_d (c, R(w-c)) =0.$$
Using implicit derivation we have:
$$\frac{\partial c}{\partial w} = - \frac{ -R u_{cd} + R^2 u_{dd}}{-u_{cc} + u_{cd} R + u_{dc} R -u_{dd} R^2} = \frac{R(u_{cd} - R u_{dd})}{-\underbrace{u_{cc}}_{<0} + 2 R u_{cd} - R^2 \underbrace{u_{dd}}_{<0}} >0$$
$$u_{cd}- R u_{dd} >0.$$
Since, in equilibrium $R = \frac{u_c}{u_d}$, the condition becomes
$$\frac{u_{cd}}{u_c} - \frac{u_{dd}}{u_d} > 0.$$
Consumption when old
In this case, we write the optimality condition as
$$ - u_c (w-\frac{d}{R}, d) + R u_d(w-\frac{d}{R}, d).$$
Hence, we can compute the derivative
$$\frac{\partial d}{\partial w} = - \frac{-u_{cc} + R u_{cd}}{u_{cc} \frac{1}{R} - u_{cd} - u_{cd} + u_{dd} R} = \frac{R (R u_{cd} - u_{cc})}{-\underbrace{u_{cc} +2Ru_{cd} - R^2 u_{dd}}_{>0}} > 0$$
$$R u_{cd} - u_{cc} > 0 \implies \frac{u_{cd}}{u_d}-\frac{u_{cc}}{u_c} > 0.$$
Savings, wages and interest rate
We now show that savings are increasing in wages, and that they may increase or decrease with the interest rate.
First, we compute the relationship between savings and wages using the optimality condition $-u_c(w-s, Rs) + u_d(w-s, Rs)$:
$$\frac{\partial s}{\partial w} = - \frac{- u_{cc} + Ru_{cd}}{u_{cc}-u_{cd}R - u_{cd}R + u_{dd}R^2}= \frac{\overbrace{u_{cc}-Ru_{cd}}^{<0,, \mathrm{normal, goods}}}{\underbrace{u_{cc}-2u_{cd}R+u_{dd}R^2}_{<0}} \gt 0 $$
Using the same optimality condition we proceed with the relationship between savings and the interest rate:
$$\frac{\partial s}{\partial R} = - \frac{-u_{cd}s+Ru_{dd}s+u_d}{u_{cc}-2Ru_{cd}+R^2u_{dd}} = \frac{u_d + s(Ru_{dd}-u_{cd})}{\underbrace{-u_{cc} + 2Ru_{cd}-u_{dd}R^2}_{>0}}.$$
Hence, the sign of the derivative depends on
$$u_d + s(R u_{dd} - u_{cd}) \gtreqqless 0 \implies 1 + \underbrace{s}_{\frac{d}{R}}\left(R \frac{u_{dd}}{u_{d}} - \frac{u_{cd}}{u_{d}}\right) \gtreqqless 0 \implies $$ $$1 + \underbrace{\frac{u_{dd}}{u_{d}}d}_{-\frac{1}{\sigma}} - \frac{u_{cd}}{u_{d}}\frac{d}{R} \gtreqqless 0.$$
The relationship depends on the intertemporal elasticity of substitution ($\sigma$) and on the degree of complementarity between $c$ and $d$ as captured by $u_{cd}.$