(fitteia's extended models library)
ROUSE(double f,double a,double tauc)
{
double w,af;
w = 2*pi*f;
af = -a/log(w*tauc);
return af;
}