(fitteia's extended models library)
double ODFN(double f, double Aodf, double fcm, double fcM, double delta)
{
double odfN();
double wcm,wcM,w,af;
wcm = 2*pi*fcm;
wcM = 2*pi*fcM;
w = 2*pi*f;
af= Aodf*odfN(wcm,wcM,w,delta);
return af;
}