double CROSSRELAX(double f, double a, double tau, double f0)

{

double af;

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

return af;

}