static double calcgmuon (int csboson,double mmuon, double mbosino, double msfermion,double coeff1, double coeff2)
/*compute the contribution for the g-2 from 1 graph
The parameters are
csbosons : charge of the vector boson
mmuon : muon's mass
mbosino : sboson's mass
msfermion : sfermion's mass
coeff1 : 1st coefficient (for the term prop to mmuon )
coeff2 : 2nd coefficient (for the term prop to mbosino )
*/
{
double x,result;
x = pow ((mbosino/msfermion),2);
result = mmuon/(16*pow(M_PI,2))*(mmuon/(12*pow(msfermion,2))*coeff1*integr(csboson,x)+(2*mbosino)/(3*pow(msfermion,2))*coeff2*integr(csboson+2,x));
return(result);
}
double test(int n)
{
return integr(square, 0.0, 1.0, n);
}
int main (){
printf ("zzz\n");
printf("int = %f\n", integr(&f1, -1, 1, 5, 1.e-12, 20));
printf ("!!!\n");
}
