Steady states and eigenvalues of equation in MATLAB

Question

I tried to solve the given problem above with the following code.

syms y r K theta
eqn = r*y*(y/theta - 1)*(1-y/K) == 0;
S=solv