Sorry, i can`t :( So, if fsolve can not find solutions it means that there are no real solutions?

Take a look please, can`t solve this problem "Warning, solutions may have been lost" Tk:=2000; > Ta:=800; > Tcs:=458.377; > p:=6.65; > js:=4.044*10^(-3); > f1:=3.925; > f2:=1.76; > T:=(Ta+Tk)/2; > d:=10^(-3); > g1:=1; > g0:=2; > gion:=1; > q:=1.6*10^(-19); > k:=1.38*10^(-23); > M:=133*1.66*10^(-27); > m:=9.11*10^(-31); > vi:=(8*k*T/(3.14*M))^(1/2); > N0:=p/k*T; > Eion:=3.89; > E1:=23.15; > U:=3; > eqn1:=f2-fi2-f1+fi1-U=0; eqn2:=js-1/4*q*n*ve*exp((-1)*q*fi1/(k*Te))-1/4*q*n*ve*exp((-1)*q*fi2/(k*Te))+0.61*q*n*sqrt(k*Te/M)*exp((-1)*q*fi1/(k*Te))-q*G*d=0; eqn3:=js/q*(q*fi1+2*k*Te)-1/4*n*ve*exp((-1)*q*fi1/(k*Te))*(2*k*Te-q*fi1)-1/4*n*ve*exp((-1)*q*fi2/(k*Te))*(2*k*Te-q*fi2)-G*d*(Eion-f1)=0; eqn4:=0.61*n*sqrt(k*Te/M)*exp(-q*fi1/(k*Te))-1/4*vi*n-G*d=0; eqn5:=n*ve*sigma1*N0*g1/g0*exp(-E1/(k*Te))+n*ve*sigma0*N0*exp(-Eion/(k*Te))-G=0; eqn6:=sigma1-gion/g1*32*3.14^2/(9*ve^2)=0; eqn7:=sigma0-gion/g0*32*3.14^2/(9*ve^2)=0; eqn8:=ve-(8*k*Te/(3.14*m))^(1/2)=0; > s:=solve({eqn1,eqn2,eqn3,eqn4,eqn5,eqn6,eqn7,eqn8},{fi2,fi1,n,Te,ve,G,sigma0,sigma1}); Warning, solutions may have been lost s :=