%打印正态分布表  
%$作者：plp626 时间：2008-5-23 $
function f()
p=0:9;                                        %第2位小数
fprintf('   z  ')
for i=1:10
fprintf('%-8.1d',p(i))
end
fprintf('\n')
for z=0:0.1:2.0                               %第一位小数
   fprintf('%4.1f',z)
   for i=1:1:10
      fprintf('%8.4f',int_f(z+p(i)/100))
      if i==10 fprintf('\n')
      end

   end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function r=int_f(a)
syms x
r=double(int(1/sqrt(2*pi)*exp(x^2/-2),-inf,a));