function jifenjisuan
clear all;clc
format long
u=0.00025;
o=2.4934e-16;
k=[0 1 2 3];
for i=1:4
mk(i)=quadgk(@(x) x.^k(i)/sqrt(pi*o).*exp(- (x-u).^2/2/o^2   ),0,Inf  );
end
mk