|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我今天在做一个积分
syms x r E1fx D1fx;
E1fx=15 ;
D1fx=5;
f=1/((2*pi)^(1/2)*D1fx)*exp(-(x-E1fx)^2/(2*D1fx^2))*(1-exp((4619.4*r+28.15*x)/r^2))*(r/4*exp(-r^2/32));
int(int(f,x,(-3-16.16*r),(2100-0.268*r)),r,0,6)
运行的结果是
Warning: Explicit integral could not be found.
> In sym.int at 58
ans =
int(28746820885731655/576460752303423488*r*exp(-1/32*r^2)*2^(1/2)*pi^(1/2)*erf(9/5*2^(1/2)+202/125*r*2^(1/2))-28746820885731655/576460752303423488*r*2^(1/2)*pi^(1/2)*erf(1/1000*2^(1/2)*(1800*r^2+1616*r^3+70375)/r^2)*exp(-1/160*(5*r^6-739104*r^3-1584845-67560*r^2)/r^4)-28746820885731655/576460752303423488*r*exp(-1/32*r^2)*2^(1/2)*pi^(1/2)*erf(-417/2*2^(1/2)+67/2500*r*2^(1/2))+28746820885731655/576460752303423488*r*2^(1/2)*pi^(1/2)*erf(1/5000*2^(1/2)*(-1042500*r^2+134*r^3+351875)/r^2)*exp(-1/160*(5*r^6-739104*r^3-1584845-67560*r^2)/r^4),r = 0 .. 6)
请教高人如何解决,万分感谢等待回复! |
|