Using for loop to change different input while solving silmultaneous equations

I want to solve the simultaneous equation with one different input value. Help me to solve this! 

Below is the command that I used:

clear

clc

%Energy balance in hot stream

%IN-OUT-ACCUMULATION=0

%m_HW_IN*c*T_HW_IN - m_HW_OUT*c*T_HW_OUT - h*A*(T_HW_OUT-T_CW_OUT)=0

%m_HW_IN*c*(T_HW_IN-T_HW_OUT) - h*A*(T_HW_OUT-T_CW_OUT)=0

%Energy balance in cold stream

%IN-OUT+ACCUMULATION=0

%m_CW_IN*c*T_CW_IN - m_CW_OUT*c*T_CW_OUT + h*A*(T_CW_OUT-T_CW_OUT)=0

%m_CW_IN*c*(T_CW_IN-T_CW_OUT) + h*A*(T_HW_OUT-T_CW_OUT)=0

%Constant

m_HW_IN=132.93;%unit=g/s

c=4.186;%unit=kJ/kg*°C

T_HW_IN=60;%unit=°C

T_CW_IN=33;%unit=°C

h=9000;%unit=W/m^2*°C

A=16*2*pi*0.0079375*4.06;%unit=m^2

syms T_CW_OUT T_HW_OUT

T_HW_OUT=zeros(5,1);

T_CW_OUT=zeros(5,1);

sol=zeros(5,1);

for m_CW_IN=[83.08;99.70;116.32;132.93;149.55];

eqn1=m_CW_IN.*c.*(T_CW_IN-T_CW_OUT)+h.*A.*(T_HW_OUT-T_CW_OUT)==0;

eqn2=m_HW_IN.*c.*(T_HW_IN-T_HW_OUT)-h.*A.*(T_HW_OUT-T_CW_OUT)==0;

sol = solve([eqn1,eqn2],[T_CW_OUT,T_HW_OUT]);

T_HW_OUT= sol.T_HW_OUT;

T_CW_OUT= sol.T_CW_OUT;

T_HW_OUT=vpa(T_HW_OUT,4)

T_CW_OUT=vpa(T_CW_OUT,4)

end