контрольные работы 7 вариант / 9 / mont_karlo

function z=mont_karlo(f,a,b,N)
    Ymin=0;
    Ymax=feval(f,b);
    x=a+(b-a)*rand(N,1);
    y=Ymin+(Ymax-Ymin)*rand(N,1);
    s=0;
    for i=1:N
        if y(i)<=feval(f,x(i))
            s=s+1;
        end;
    end;
    z=s*(b-a)*(Ymax-Ymin)/N;
end