И другие программы этой серии
Допустимые значения для параметра
order – @quad , @quadl или имя любого определенного пользователем квадратурного метода с таким же вызовом и такими же возвращаемыми параметрами, как у функций quad и quadl. (Например, при проверке старых программ можно использовать @quad8 для большей совместимости с прежними версиями MATLAB). По умолчанию (без параметра order) вызывается @quad, поскольку подынтегральные функции могут быть негладкими.
Пример: пусть m-файл integ1.m описывает функцию 2*y*sin(x)+x/ 2*cos(y), тогда вычислить двойной интеграл от этой функции можно следующим образом:
>> result = dblquad(@integ1,pi,2*pi,0,2*pi) result = -78.9574
В систему MATLAB 6.5 была введена новая функция для вычисления тройных интегралов. Она имеет четыре формы записи:
triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax) triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol) triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol,method) triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol,method,p1,p2,...)
Параметры этой функции уже были определены выше (добавлены только пределы для новой переменной z). Пример:
>> Q = triplequad(\'x+y*z\',0,pi,0,1,-1,1,0.001) Q = 9.8696
К сожалению, проверка данной функции на ряде тройных интегралов показала, что она не всегда обеспечивает вычисления. Возможны ситуации, когда данная функция выводит ряд сообщений об ошибках даже для берущихся интегралов.
8.7. Математические операции с полиномами
8.7.1. Определение полиномов
Полиномы (у нас их принято называть также степенными многочленами) – широко известный объект математических вычислений и обработки данных. В MATLAB полином записывается в виде
p(x) = a1xn + a2xn–1 + … + anx + an+1.
В дальнейшем мы ограничимся рассмотрением полиномов с целыми положительными степенями.
Широкое применение полиномов отчасти обусловлено большими возможностями полиномов в представлении данных, а также их простотой и единообразием вычислений. В этом разделе описаны основные функции для работы с полиномами.
Начало в части 1