Imaginary solution with symbolic toolbox while solving a system of non-linear equations

Hi everyone, I'm trying to use the symbolic toolbox for modeling a non-linear system, but I get an imaginary solution. The physics behind it are verified using valid references but I guess I'm making a mistake in the matlab implementation of it. As an alternative, I tried to use the optimization toolbox (fsolve) to get the solutions, but that gives an imaginary solution as well. I'm 100% confident that the system equations are correct, but I get a problem with the solutions in matlab. Also some values are negative, which isn't supposed to happen either.

I've uploaded the m-file for reference. I've also added the system equations as comments inside the loop, along with the equilibrium conditions that I'm solving for.

Thanks in advance.