RMTDnum(double f, double p, double Armtd, double fcm,double fcM)

{

/****

numerical calculation of the spectral densities as presented in

M. Vilfan,T. Apih, P.J. Sebastão, G. Lahajnar,

S. Zumer, Phys. Rev. E, 76, 51708(2007)

****/

double x;

double sqgaus(),sqromo(),smidpnt(),jrmtd(),A,j1,j2,af;

void w_f_ptr(),clear_struct();

af =pow(1.0/f,p)*jrmtd(f,fcM,fcm,p) + 4*pow(0.5/f,p)*jrmtd(2*f,fcM,fcm,p);

/* printf("af=%le\n",af); */

return Armtd*af;

}