examp9_14.m

来自「东北大学 薛定谔的高等应用数学问题MATLAB求解的源代码」· M 代码 · 共 7 行

syms x; syms a lam positive; p=lam^a*x^(a-1)/gamma(a)*exp(-lam*x);
for n=1:5, m=int(x^n*p,x,0,inf), end

syms n; m=simple(int((x)^n*p,x,0,inf))

for n=1:7, s=simple(int((x-1/lam*a)^n*p,x,0,inf)), end