double LORTZN(double f, double a, double b, double f0)

{

double af;

af= a/(3.14159265*(1.0+pow((f-f0)/b,2.0)));

return af;

}