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画积分图像
matlab 积分图像可以来实现。其方法利用for循环语句求解每一个(d,源码哲程crm源码y)的源码龙ol免费源码对应值,然后用plot绘出其积分图像。源码360查看页面源码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为复数值。源码粉粉说源码请题主检查一下积分方程是源码否有错。
源码源码