D=0.05; %m
V_Ice=10; %m^3
T_mi=24;%C
Ts=0;%C
f_vi=0.8;%
N=10; %Tubes
L=1.55 ;%m
m_dot=linspace(0.005,0.05); %kg/s % <------- changed here. Earlier it was a scalar
roh_ice=920; %kg/m^3
Q=3.34*10^5; %J/kg
cp=1007; %J/kg*K
h_bar=22.3; %W/m^2*K
T_mo=Ts-(Ts-T_mi)*exp(((-pi*D*L)./(m_dot*cp))*h_bar) ;%C
T_mo(T_mo==0)=[];
figure
plot(m_dot,T_mo)
xlabel('mass flow rate')
xlim([0.005 0.05])
ylim('auto')
ylabel('T_m,o')
title('Out-put Temperature of Air vs Mass Flow Rate of Air')
