area1.m

随机过程中的possion算法在matlab中的实现例程

%求到达一定时刻x时的总人数y

function[y]=area1(x)
if(x<=8)
    y=1/2*(x-5)*[write(x)-200]+200*(x-5);
else
    if(x<=18)
        y=area1(8)+(x-8)*1400;
    else
        if(x<=21)
            y=area1(18)+200*(x-18)+1800-1/2*(21-x)*(-400*x+8600);
        end
    end
end