I'm trying to plot phase analysis for the Friedmann's equation, however when running the protocol and trying to plot these phase analysis, it comes up with error that rho wasn't allowed in the set of options!

Here's the code i'm having difficulty with!?!

> SFsys := proc (lambda, W, ics, T, v, u) local sys, flatF, rho0line, plot0, plotDSp, plotDSm, plotE, flatF_Plot, IC, df1, df2; 

sys := eval(DSys, {Lambda = lambda, gamma = W});