You are setting that population to 500 initially in u0 (so it is not starting at zero), and adding 500 again at t==1000, but by then the population appears to be so large that 500 would not make a visible difference (if I’m lining up the state variables to the plot lines correctly).