大家快来看看啊，有懂的帮帮啊
<FONT size=7>我积分了怎么出来个这啊<BR>WhittakerM是什么啊</FONT><BR><BR>>> clear;syms x l m A h;<BR>>> int(x^(2*l+2)*exp((-(2*m*A)*x)/((l+1)*h^2)),x)<BR> <BR>ans =<BR> <BR>2^(-3-2*l)*(m*A/(l+1)/h^2)^(-2*l)/m^3/A^3*(l+1)^3*h^6*(2^(l+1)*x^(2*l)*(m*A/(l+1)/h^2)^(2*l)*(l+1)*(m*A*x/(l+1)/h^2)^(-l)*exp(-m*A*x/(l+1)/h^2)*WhittakerM(l,l+1/2,2*m*A*x/(l+1)/h^2)-x^(2*l)*2^l*(m*A/(l+1)/h^2)^(2*l)*(2*l+2+2*m*A*x/(l+1)/h^2)*(m*A*x/(l+1)/h^2)^(-l)*exp(-m*A*x/(l+1)/h^2)*WhittakerM(l+1,l+1/2,2*m*A*x/(l+1)/h^2))<BR> <BR> <BR>>>回复：（sxw189）大家快来看看啊，有懂的帮帮啊
这个要问你啊<BR>WhittakerM应该是自己定义的矩阵或者函数 如果设a=2(l+1)b=-4mA/h^2则积分写成 x^a*exp(b/a*x)结果是-(-b/a)^(-a)/b*a*(a*x^a*(-b/a)^a*gamma(a)*(-b/a*x)^(-a)-x^a*(-b/a)^a*exp(b/a*x)-a*x^a*(-b/a)^a*(-b/a*x)^(-a)*gamma(a,-b/a*x)) 这里有个gamma(a,-b/a*x)这个函数，我试了一下diff gamma(a,b*x)ans = -(b*x)^(a-1)*exp(-b*x)*b
所以它应该是把一般的gamma函数的积分"区间"(0,+inf)换为(0,x)
另外：diff WhittakerM(a,b,x)
ans =(1/2-a/x)*WhittakerM(a,b,x)+(1/2+b+a)/x*WhittakerM(1+a,b,x)
另外:y''+A(a*a-x*x)y=0的解是：C1/x^(1/2)*WhittakerM(1/4*A^(1/2)*a^2,1/4,A^(1/2)*x^2)+_C2/x^(1/2)*Whittak
erW(1/4*A^(1/2)*a^2,1/4,A^(1/2)*x^2) （http://www.mitbbs.cn/article_t/Mathematics/3106820.html）
页:
[1]