Hello, I suspect there are no closed-form solutions for this equation. Given admissible X values, however, the corresponding Y values can be determined, as shown in the plot below:
ListLinePlot[
Table[{x,
y /. NSolve[(4.208 10^-7) (9.98687) y^2/x^2 ==
1 - y - ((y^1.58273)/(x^0.58273)), y, Reals][[1]]}, {x, .1, 20,
0.1}], AxesLabel -> {"x", "y"}, LabelStyle -> Black]
Which will generate
