>> syms x X b k positive; >> p(x)=x*exp(-x^2/2/b^2)/b^2;P(x)=int(p(x),x,0,X); >> Ex=int(x*p(x),x,0,inf);Dx=int((x-Ex)^2*p(x),x,0,inf); >> vk=int(x^k*p(x),x,0,inf);miuk=int((x-Ex)^k*p(x),x,0,inf);