dxdn[x_, y_] := (-1 + x)*(g*x - a + 3*y);
dydn[x_, y_] := y/x*(-1 + x)*(g*x - a + 3*y);
Solve[{dxdn[x, y] == 0, dydn[x, y] == 0, y}, {x, y}]