% programa de cálculo del coeficiente de descarga
function y=descargain(lii)
Cda=[0.55 0.67];
Cdb=[0.57 0.66];
Cdc=[0.66 0.45]; 
li1=[0 0.1];
li2=[0.1 0.175];
li3=[0.175 0.3];
if lii<=0.1
   y=INTERP1(li1,Cda,lii,'spline');
elseif lii<=0.175
   y=INTERP1(li2,Cdb,lii,'spline');  
else
   y=INTERP1(li3,Cdc,lii,'spline');
end
