DSolve[{Derivative[1][p][t] == -nstar[t] p[t] + (1 - p[t]/10000) p[t] - p[t]/(1 + p[t]),
Derivative[1][nstar][t] == -nstar[t] + (nstar[t] + p[t])/(1 + nsatr[t] + p[t])}, {p, nstar}, t]
The output is
DSolve[{Derivative[1][p][
t] == -nstar[t] p[t] + (1 - p[t]/10000) p[t] - p[t]/(1 + p[t]),
Derivative[1][nstar][t] == -nstar[t] + (nstar[t] + p[t])/(
1 + nsatr[t] + p[t])}, {p, nstar}, t]
The functions are p(t) and nstar(t).
Thank you for your help.