The code is the following
Clear["Global`*"]
kpg = 1
pinf = 20000000
kpm = (.005 - um)
sm = .003 /(.005 - um)
um = .002
kmp = .01
kpn = 1.8
snr = 0.08
knn = .01
knp = .1
unr = .12
un = 0.05
eq2 = snr (knn nstar[t] + knp p[t])/(unr + knn nstar[t] + knp p[t]) -
un nstar[t] == 0
b = Solve[eq2, nstar[t]]
It shows an error message saying
Solve was unable to solve the system with inexact coefficients. The answer was obtained by solving a corresponding exact system and numericizing the result. >>
The solution is
{{nstar[t] ->
0.2 (-26. - 25. p[t] -
1. Sqrt[676. + 1700. p[t] + 625. p[t]^2])}, {nstar[t] ->
0.2 (-26. - 25. p[t] + Sqrt[676. + 1700. p[t] + 625. p[t]^2])}}
