I don't think Solve
can handle this system. It has transcendental dependencies (use of both ny
and Log[ny]
, for example). So you might have to restrict to numeric rather than symbolic parameters, and either use FindRoot
or specify for Solve
to work over the Reals
(and it might still fail).