Matlab da İntegral Hesaplayan Algoritma

function SONUC=integral(a,b,n)
fmax = 1/3;
k=0;
for i=1:n
Ux=a+(b-a)*unifrnd(0,1);
fux =(1/3)*(x^2);
uy=fmax*unifrnd(0,1);
if uy<=fux
k=k+1;
end
end
Sonuc=(k/n)*(b-a)*fmax;
end

Buradaki fmax değeri integralini alacağınız fonksiyonun maksimum değeridir. Buradaki fonksiyon fux=(1/3)*(x^2) ‘ dir.