 			* (cyl_bessel_i0(y) + (k2*f+k1*(1.-f))*tau
 				* cyl_bessel_i1(y)/y); // Full Expression
 	}
-
-}
-
+}
 
 /*
 	Parameters ripped from Quadpack dqk21.f on 3/20/2017:

	c as evaluated with 80 decimal digit arithmetic by l. w. fullerton,
	c bell labs, nov. 1981.
*/


__constant__ static double wg[5] = {
	0.066671344308688137593568809893332,

	cudaSetDevice(device);
	cudaDeviceProp deviceProp;
	cudaGetDeviceProperties(&deviceProp, device);
	int threads = deviceProp.maxThreadsPerBlock/4.;
	int blocks = (N+threads-1)/threads;

	double * d_d;
