#include #include double cal(double r1,double r2,double dx) { return sqrt(r1*r1-dx*dx)*sqrt(r2*r2-dx*dx); } double fun(double r1,double r2) { int n=1000000; double dx=(r1