using solve function and getting 0-by-1

Hello!

I have a system of 5 nonlinear equations(eq1...eq5) where w,Po,An,Mu,S are my variables, and another are costant parametes .I try to solve with solve but I cant, but i get this: " Empty sym: 0-by-1 " for all variable.

in addition i get this message: Warning: {Cannot find explicit solution. > In solve (line 316)}

i am beginner in matlab and cant solve yhis problem, but im sure that the system has answer. Can anybody help me pls?

syms Rho  bo br beta delta gamma c alpha  Phi Eta K %parameters
syms w Po Pr An Al Mu S %decision variables

Pr =(delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho);
Al =(delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2;

eq1 =Phi*S + K*alpha - Po*bo - (Po - c)*(bo + Rho*((2*Rho - 2*Mu*Rho)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + 1) - (Rho*delta^2*gamma*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) + An^(1/2)*beta - (c - w)*(Rho*((2*Rho - 2*Mu*Rho)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + 1) + (br*(2*Rho - 2*Mu*Rho))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (Rho*delta^3*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) + gamma*((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2) - Rho*(Po - (delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)) - (2*Mu*Rho*delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2
  
eq2 =An^(1/2)*beta - alpha*(K - 1) + (Po - c)*((Rho*(delta^2 - 2*Rho - 2*br + 2*Mu*br + 2*Mu*Rho))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) - (delta^2*gamma*(Rho + br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) + (c - w)*((Rho*(delta^2 - 2*Rho - 2*br + 2*Mu*br + 2*Mu*Rho))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (br*(delta^2 - 2*Rho - 2*br + 2*Mu*br + 2*Mu*Rho))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (delta^3*(Rho + br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) + delta*((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2) + Rho*(Po - (delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)) - (br*(delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (2*Mu*delta^2*(Rho + br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2

eq3 =(Po - c)*(beta/(2*An^(1/2)) - (Rho*(beta/An^(1/2) - (Mu*beta)/An^(1/2)))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (beta*delta^2*gamma*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(2*An^(1/2)*((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) - (c - w)*(beta/(2*An^(1/2)) + (Rho*(beta/An^(1/2) - (Mu*beta)/An^(1/2)))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (br*(beta/An^(1/2) - (Mu*beta)/An^(1/2)))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) + (beta*delta^3*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(2*An^(1/2)*((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)) - (Mu*beta*delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta))/(An^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2) - 1
  
eq4 =(Rho*((2*Po*Rho + 2*K*alpha + 2*Rho*w + 2*br*w + 2*An^(1/2)*beta)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) - ((4*Rho + 4*br)*(delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2) - (delta^2*gamma*(4*Rho + 4*br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^3))*(Po - c) + (c - w)*(Rho*((2*Po*Rho + 2*K*alpha + 2*Rho*w + 2*br*w + 2*An^(1/2)*beta)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) - ((4*Rho + 4*br)*(delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2) + (br*(2*Po*Rho + 2*K*alpha + 2*Rho*w + 2*br*w + 2*An^(1/2)*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho) - (br*(4*Rho + 4*br)*(delta^2*w - 2*K*alpha - 2*Rho*w - 2*br*w - 2*An^(1/2)*beta - 2*Po*Rho + 2*Mu*Po*Rho + 2*K*Mu*alpha + 2*Mu*Rho*w + 2*Mu*br*w + 2*An^(1/2)*Mu*beta))/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2 + (delta^3*(4*Rho + 4*br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(((delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2)^(1/2)*(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^3)) - (delta^2*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^2 + (2*Mu*delta^2*(4*Rho + 4*br)*(Po*Rho + K*alpha - Rho*w - br*w + An^(1/2)*beta)^2)/(delta^2 - 4*Rho - 4*br + 4*Mu*br + 4*Mu*Rho)^3
 
eq5 =Phi*(Po - c) - Eta*S

Ans=solve(eq1,eq2,eq3,eq4,eq5,Po,w,An,S,Mu);
Ans.Po
Ans.An
Ans.w
Ans.S
Ans.Mu

[Aman Banthia]

Hi,

Seems like the equation is too heavy for MATLAB to handle. I have tried to ease off some load but nothing seems to work.

Ans=solve(eq1,eq2,eq3,eq4,eq5,Po,w,An,S,Mu,'IgnoreAnalyticConstraints', true,'Real', true);

Thankyou