Добавил:
NLapina
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:контрольные работы 7 вариант / 9 / mont_karlo
.m 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
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
Соседние файлы в папке 9