Solar heat collector simulation

Hello, I'm trying to simulate the heat transfert in solar heat collector,i have 3 equastions with 3 variables, i used the Solve command to find the solutions of x,y and z, but the solutions are not right and illogic, so is the used programme right and how can i use the Solve command is this situation??

clear; clc;
v_vent=5;
h=0.01;%l'epaisseur de l'air confiné
Tamb=273+30;
T_abs(1)=Tamb + 15;
T_v(1)=Tamb + 10;
n=10;
Al=1/h;
alfa=pi/4;
h_cv_air = 5.7+3.8*v_vent;
for i= 1:n
    Gr(i)=10*0.0037*(T_abs(i)-T_v(i))*h^3/((15*10^-6)^2);
    Nul1(i)=1+1.89*10^-4*Gr(i)*Al^(-3/4);
    h_ab_v1(i) = Nul1(i) * 0.025/h;
    Ral(i)=Gr(i)*0.7;
    Nul2(i) = 1+1.44*(1-(1708/(Ral(i)*cos(alfa))))*(1-(1708*sin(1.8*alfa)^(1.6)/(Ral(i)*cos(alfa))))*Al^(-1/4)+(((Ral(i)*cos(alfa)/5830)^(1/3))-1);
    h_abs_v2(i)=Nul2(i)*0.025/h;
    h_cv_abs_v(i)=(h_ab_v1(i)+h_abs_v2(i))/4;
    h_ryt(i) =5.67*10^-8*(T_abs(i)^2+T_v(i)^2)*(T_abs(i)+T_v(i))/((1-0.97)/0.97+1+(1-0.93)/0.93);
    h_abs_lin(i) = h_cv_abs_v(i)+h_ryt(i)
    h_ray_air(i)=5.67*10^-8*0.93*(T_v(i)^2+Tamb(i)^2)*(T_v(i)+Tamb(i));
    h_v_air(i) = h_cv_air(i)+h_ray_air(i)
    syms x y z
    e1(i)=h_v_air(i)*1.5*(x-303)-1.5*0.02*1000-1.5*0.02*5.67*10^-8*(0.0552*303^1.5)^4-h_cv_abs_v(i)*1.5*(y-x)==0; 
    e2(i)=1.5*0.97*0.83*1000-1.5*h_cv_abs_v(i)*(y-x)-2*100*(y-(z+303)/2)==0;
    e3(i)=1.5*0.97*0.83*1000-0.027*4280*(z-303)-0.9991*1.5*((z+303)/2-303)-h_v_air(i)*1.5*(x-303)==0;
    eqs=[e1(i), e2(i), e3(i)];
    [x,y,z]=solve(eqs,[x,y,z])
    %T_abs(i+1)=y;
    %T_v(i+1)=x;
    %T_f_s(i)=z;
end