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

{

double af;

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

return af;

}