% programa de cálculo del coeficiente de descarga
function y=descargaes(lii)
Cd=[0.7 0.73 0.75 0.6];
li=[0 0.1 0.2 0.3];
   y=INTERP1(li,Cd,lii,'spline');
