%función para el cálculo de las perdidas de calor
function y=heat_loses1(T,p,T1,p1,V1,p0,V0,TW,V,A,Vs,Vc,gamma,vp,D)
vmot=(2.28*vp);
pmot=p0*(V0/V)^gamma;
vcomb=0.00324*(Vs-Vc)*(T1/(p1*V1))*(p-pmot);
h=0.00326*(p^0.8)*((vmot+vcomb)^0.8)/((D^0.2)*(T^0.53));
y=h*A*(T-TW);
