Лекции по Matlab

     >> xn=[1:0.1:1.5,1.6:0.2:2];

     >> yn=log(xn);

     >> trapz(xn,yn)

     ans =

         0.3856

     Более точная формула численного интегрирования с равноотстоящими узлами – формула  Симпсона:

     Она является точной для многочленов степени не выше 3 (из вида остаточного члена).

     Для вычисления интеграла по формуле  Симпсона в MATLAB используется функция quad. Она вычисляет интеграл с заданной точностью. Для этого использует рекурсивные вызовы самой себя с дроблением шага разбиения.

     Работает  только для аналитически заданной функции. Форма обращения: q=quad(fun,a,b,e), где fun – указатель на подынтегральную функцию, a,b – пределы интегрирования, e –точность (по умолчанию 10^-6).

     >> quad('log(x)',1,2)

     ans =

         0.3863