ЧМ (ЭКТ-3) / Лабы / Всякие сделанные варианты / other / integr

n=500

x1=0

x2=1

h=(x2-x1)/n

S=(log(1+x1)/(1+(x1^2)))+(log(1+x2)/(1+(x2^2)))

xi=x1

for I=1:n-1

xi=xi+h

if rem(I,2)==1

S=S+(4*(log(1+xi)/(1+(xi^2))))

else

S=S+(2*(log(1+xi)/(1+(xi^2))))

end

end

otvet_SIMPSON= (S*h)/3

x=rand(1,n)

y=rand(1,n)

n1=0

for I=1:n

if y(I)<(log(1+x(I))/(1+(x(I)^2)))

n1=n1+1

end

end

otvet_MONTE=n1/n

otvet_SIMPSON

otvet_MONTE

otvet_SIMPSON - otvet_MONTE

x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]

for I=1:size(x,2)

y(I)=(log(1+x(I))/(1+(x(I)^2)))

end