Квадрату́рные фо́рмулы, формулы приближённого вычисления определённого интеграла от функции $f(x)$ по её значениям в конечном числе точек. Погрешность вычисления интегралов по квадратурным формулам зависит от гладкости функции $f(x)$. Обычно квадратурные формулы имеют вид $$\int\limits_a^b f(x)dx \approx p_1f(x_1)+ \ldots+p_nf(x_n),$$где $x_1, \ldots, x_n$ – точки отрезка $[a,b]$, их называют узлами квадратурных формул, а $p_1, \ldots, p_n$ – числа, их называют коэффициентами (или весами) квадратурных формул. Часто отрезок $[a,b]$ разбивают на конечное число отрезков, на каждом из которых используют какую-либо простую квадратурную формулу. Простейшими квадратурными формулами являются формула прямоугольников $$\int\limits_a^bf(x)dx \approx(b-a)f \left(\frac{a+b}{2}\right) \int f(x)dx \approx (b-a)f \left(\frac{a+b}{2}\right)$$и формула трапеций $$\int\limits_a^bf(x)dx\approx \frac {b-a}{a}(f(a)+f(b));$$они дают точное значение интеграла для линейных функций. Формула Симпсона $$\int\limits_a^bf(x)dx \approx \frac{b-a}{6} \left(f(a)+4f \left(\frac{a+b}{2} \right)+f(b) \right)$$даёт точное значение интеграла для многочленов 3-й степени.

Для приближённого вычисления кратных интегралов используются аналогичные формулы, которые называют кубатурными.

Для вычисления, в первую очередь кратных интегралов, используют также формулы со случайным выбором узлов и коэффициентов (см. Метод Монте-Карло).