Difficulty Finding Fixed Points for NPZ Model Using fsolve in Python

I’ve been working on finding the fixed points for the NPZ model described by the Franks equations: