func5p3.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 17 行

%Integration,using Simpson methom
x=linspace(180,188,9)
y=linspace(180,188,9)
for i=1:9
    y(i)=(sub5p3(x(i))+4*sub5p3((x(i)+194)/2)+sub5p3(194))/6*(194-x(i))
end
k=0
for i=1:10000
    if(6*randn+170>180)
        k=k+1
    end
end
%x=180       181       182       183       184       185       186       187       188
%y=0.0499    0.0351    0.0240    0.0159    0.0103    0.0064    0.0039    0.0023    0.0013
%（1）查表值F(180)=0.0548,F(182)=0.02275,F(188)=0.00145,与计算值相差不大
%（2）据计算结果，P(X>184)=0.0103,故车门在184cm时能满足题目要求，在180cm是有5%的市民不碰头，满足大部分人利益
%（3）多次使用randn函数，算得所占百分比为4.78%,4.72%,4.95%,与理论值相当接近