double OPF(double f, double Aopf, double f0,double fcm, double fcM,double number_of_pts_in_integral)

{

double x;

double jopf(),j1;

j1=jopf(2.0*pi*f,1.0/(2.0*pi*f0),2.0*pi*fcm,2.0*pi*fcM,number_of_pts_in_integral);

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

return Aopf*j1;

}