ROUSE(double f,double a,double tauc)

{

double w,af;

w = 2*pi*f;

af = -a/log(w*tauc);

return af;

}