angel4296
2008-09-16, 17:49
a=9.8427;
b=3.1491;
y= (a*b)^((a+b)/2)*x^((a+b)/2-1)*besseli(a-b,2*(a*b*x)^0.5)*(2/pi^0.5)*exp(-t^2);
Y1=int(int(y,t,500/707*x,inf),x,0,inf)
得到的解为
Warning: Explicit integral could not be found.
> In sym.int at 58
Y1 =
int(1621193981275812845703240236051/590295810358705651712*x^(54959/10000)*besseli(8367/1250,1/140737488355328*2455728123682660124221610393600^(1/2)*x^(1/2))*(-pi^(1/2)*erf(500/707*x)+pi^(1/2)),x = 0 .. Inf)
哪位高手能知道怎么得到Y1的近似解？？谢谢了：）
b=3.1491;
y= (a*b)^((a+b)/2)*x^((a+b)/2-1)*besseli(a-b,2*(a*b*x)^0.5)*(2/pi^0.5)*exp(-t^2);
Y1=int(int(y,t,500/707*x,inf),x,0,inf)
得到的解为
Warning: Explicit integral could not be found.
> In sym.int at 58
Y1 =
int(1621193981275812845703240236051/590295810358705651712*x^(54959/10000)*besseli(8367/1250,1/140737488355328*2455728123682660124221610393600^(1/2)*x^(1/2))*(-pi^(1/2)*erf(500/707*x)+pi^(1/2)),x = 0 .. Inf)
哪位高手能知道怎么得到Y1的近似解？？谢谢了：）