1.  double alpha = *(double *) params;

 1.  
 2. #include <stdio.h>
 3. #include <math.h>
 4. #include <gsl/gsl_integration.h>
 5. 
 6. /*
 7. The integrator QAGS will handle a large class of definite integrals. For example, consider the following integral, which has a algebraic-logarithmic singularity at the origin,
 8. 
 9. \int_0^1 a*x^{-1/2} log(x) dx = -a*4
10. 
11. The program below computes this integral to a relative accuracy bound of 1e-7. 
12. */
13. 
14. double f (double x, void * params) {
15.   double alpha = *(double *) params;
16.   double f = log(alpha*x) / sqrt(x);
17.   return f;
18. }
19. 
20. //double_int_f()
21. 
22. int
23. main (void)
24. {
25.   gsl_integration_workspace * w 
26.     = gsl_integration_workspace_alloc (1000);
27.   
28.   double result, error;
29.   double expected = -4.0;
30.   double alpha = 1.0;
31. 
32.   gsl_function F;
33.   F.function = &f;
34.   F.params = alpha
35. 
36.   gsl_integration_qags (&F, 0, 1, 0, 1e-7, 1000,
37.                         w, &result, &error); 
38. 
39.   printf ("result          = % .18f\n", result);
40.   printf ("exact result    = % .18f\n", expected);
41.   printf ("estimated error = % .18f\n", error);
42.   printf ("actual error    = % .18f\n", result - expected);
43.   printf ("intervals =  %d\n", w->size);
44. 
45.   return 0;
46. } 
47. 

1.  
2.   struct f_params * params
3.        = (struct f_params *)p;
4.   double bias = (params->bias);
5.   double temp = (params->temp);
6.   double delta = (params->delta);
7.   double Gnn = (params->Gnn);

 1.  
 2. #include <gsl/gsl_const_mksa.h>
 3. 
 4. struct f_params { double bias; double temp; double delta; double Gnn};
 5. 
 6. ....
 7. 
 8. //fonction a integrer entre 0 et +inf pour avoir I(V)/(Gnn/e)=int_0^1 f(x)dx
 9. double f (double E, void * p) 
10. {
11.   struct f_params * params
12.        = (struct f_params *)p;
13.   double bias = (params->bias);
14.   double temp = (params->temp);
15.   double delta = (params->delta);
16.   double Gnn = (params->Gnn);
17.   return Nsupra(E,delta)*(fermi(E-bias,temp)-fermi(E,temp));
18.   //return N_supra(E,delta)*(fermi(E+bias,temp)-fermi(E-bias,temp));
19.  }