%función para el cálculo de las perdidas de calor
function y=heat_loses(T,p,TW,V,A,gamma,vp,D,proceso)
if proceso ==0
   c1=2.28;
else
   c1=6.18;
end
vmot=c1*vp;
h=0.00326*(p^0.8)*((vmot)^0.8)/((D^0.2)*(T^0.53));
y=h*A*(T-TW);
