Consider a household consisting of a woman and a man, with preferences over leisure and consumption given by:
$U(\overrightarrow{c},\overrightarrow{l}) = \ln{c} + \ln{l^F} + \ln{l^M}$
where $\overrightarrow{l} = (l^F,l^M)$ is the vector of leisure times for the female and male, respectively.
On the other hand, consumption $c$ is defined by the household's preferences over consumption vectors $\overrightarrow{c} = (c^N,c^V)$ given as
$c = \frac{2}{3} c^N + \frac{1}{3} c^V$
where $c^N$ is market goods (which I'd think of as groceries) and $c^V$ is consumption generated by domestic labor (which I'd think of as pure satisfaction of having a clean home)
Consumption of domestic labor is given by (a kind of production function)
$c^V = (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}$
where $V = (V^F,V^M)$ are the domestic labor times for the female and male, respectively.
Consumption of market goods is given by the budget constraint
$c^N = w^F N^F + w^M N^M + \Pi$
where $N = (N^F,N^M)$ are the paid labor times for the female and male, respectively; $\Pi$ is the firm's profits, and $(w^F,w^M)$ are the wages for the female and male, respectively.
Here both individuals have a time endowment of $1$ unit:
$l^i + V^i + N^i = 1$
Find the female's and male's paid labor supply and domestic labor times and compare.
Here is what I have done:
Plugging $c^N, c^V$ into the preferences for $c$,
$c = \frac{1}{3} [2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}]$
$\ln{c} = - \ln{3} + \ln{[2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}]}$
From the time constraint I can eliminate the leisure variables by rewriting $l^i = 1 - N^i - V^i$
Plugging my new expressions for $c$ and the $l^i$ into the utility function:
$U(N^F,N^M,V^F,V^M) = - \ln{3} + \ln{[2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}]} + \ln{(1-N^F-V^F)} + \ln{(1-N^M-V^M)}$
From here I'd get the first order conditions by differentiating with respect to each of the four variables and setting $= 0$.
$\frac{\partial U}{\partial N^F} = \frac{2 w^F}{2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}} - \frac{1}{1-N^F-V^F} = 0$
$\frac{\partial U}{\partial N^M} = \frac{2 w^M}{2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}} - \frac{1}{1-N^M-V^M} = 0$
$\frac{\partial U}{\partial V^F} = \frac{\frac{1}{2} (V^F)^{-\frac{1}{2}} (V^M)^{\frac{1}{2}}}{2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}} - \frac{1}{1-N^F-V^F} = 0$
$\frac{\partial U}{\partial V^M} = \frac{\frac{1}{2} (V^F)^{\frac{1}{2}} (V^M)^{-\frac{1}{2}}}{2w^F N^F + 2w^M N^M + 2 \Pi + (V^F)^{\frac{1}{2}} (V^M)^{\frac{1}{2}}} - \frac{1}{1-N^M-V^M} = 0$
However, solving this system of $4$ equations seems rather complicated.
I tried pairing each two equations with the same negative term and get an equation from each pair, to reduce it to a system of two equations.
However, when I try to solve the new system, I get $w^F w^M = \frac{1}{16}$, which doesn't seem very helpful or even make sense.
Is this method of substituting the consumption equations at the start wrong here, and I should rather try to form a (very long) Lagrangian? Is there a more clever way to solve this problem?
I think I may be able to use marginal rates of substitution but this seems very confusing to me as in this situation, there are now $3$ possible ways the household agents could spend their time, rather than $2$; aside from the fact that the household consists of $2$ agents rather than one.