double CED_poli(double f, double Aecd, double fc, double R, double v0)

{

double ECD_poli();

return ECD_poli(f, Aecd, 0.0, fc, R, v0);

}

double ECD_poli(double f, double Aecd, double fcm, double fcM, double R, double v0)

{

double J1odfcolpoli();

double af,wr,wcm,wcM,w;

wcm = 2*pi*fcm;

wcM = 2*pi*fcM;

w = 2*pi*f;

wr = w/wcM;

af = Aecd*J1odfcolpoli(wr,R,v0);

return af;

}