1.acosԴ?源码?
2.用matlab画积分图像

acosԴ??

       #include<iostream>

       #include<cmath>

       #include<cstdio>

       #define PI acos(-1)

       using namespace std;

       int main()

       {

       double a1 = 4;

       double a2 = 2;

       double x = 0;

       double y = 0;

       x = (pow(a1, 2) - pow(a1, 2) * PI / 4 + pow(a1, 2) * PI / 4 / 3 - a1 - 2 * (sqrt(pow(a1, 2) - pow(a2, 2)) - 2)) / 2;

       y = pow(a2, 2) - pow(a2, 2) * PI / 4 - x;

       double s = fabs(x - y);

       printf("%.4lf\n", s);

       return 0;

       }

用matlab画积分图像

       d0=0:0.1:3.1;

       n=length(d0);

       for i=1:n

       d=d0(i);

       R=0.3*d+1;

       syms x

       eq=2*x/R^2.*exp(-x.^2/R^2).*acos((x.^2+0.1)./(2*x)); %eq=@(x)2*x/R^2.*exp(-x.^2/R^2).*acos((x.^2+0.1)./(2*x));

       y(i)=int(eq,源码千群助手源码0.1,R); %y(i)=integral(eq,0.1,R,'RelTol',1e-8,'AbsTol',1e-);

       end

       plot(d0,y)

       xlabel('d'),ylabel('y');grid on

       运行结果

注：当d0>3.1时，y为复数值。源码cfjt指标源码优化请题主检查一下积分方程是源码否有错。