P(i)=(((md*pi.*R.^2).^n(i))./factorial(n(i))).*exp(-md*pi*R.^2); n(i)替换为i P(i)=(((md*pi.*R.^2).^i)./factorial(i)).*exp(-md*pi*R.^2);
