a = 1.5;
b = 3.0;
x = a:b;

function y = func(x)
y=1 / (x.*log(x));
endfunction

funcprot(0);
v = inttrap(x, func(x));

disp("Решение интеграла методом трапеций: ");
disp(v);

function s = simpson(a, b, n)
    s = 0;
    h = (b-a) ./ (n-1);
    s = func(a) + func(b);
    mn = 4;
    for m = 1:1:n-2
        s = s + mn .* func(a + h .* m);
        if (mn == 4) then mn = 2;
            else mn = 4;
         end
    end
    s = s .*h ./3;
endfunction

n = 2;
r = simpson(a, b, n);
r2 = 0;
while (abs(r - r2)>0.001)
    r2 = r;
    n = n + 2;
    r = simpson(a, b, n);
end

disp("Решение интеграла методом Симпсона: ");
disp(r);

disp("Решение интеграла методом Симпсона(встр. ф-ция): ");
disp(intsplin(x,func(x)));