>>573
台形公式でよければ >>212 を少し変更すればおｋ
return 4/(1+x*x);
↓
return sin(2*x)*cos(x);

printf("%f\n", f_integral(0, 1, 1000));
↓
printf("%f<x<%f %d %f\n", 0., M_PI/4, n, f_integral(0., M_PI/4, n));