q = integral2(fun, xmin, xmax, ymin, ymax) intègre numériquement la fonction fun(x, y) sur le domaine xmin <= x <= xmax et ymin(x) <= y <= ymax(x).

Les bornes en y ymin et ymax peuvent être des scalaires ou des handles de fonction de x pour décrire un domaine non rectangulaire.

L'intégration utilise une quadrature de Gauss-Kronrod adaptative imbriquée. Les paires nom/valeur RelativeTolerance (défaut 1e-6) et AbsoluteTolerance (défaut 1e-10) contrôlent la précision ; les anciens noms RelTol et AbsTol restent acceptés.

Par défaut (Vectorized à true), fun doit accepter des tableaux et opérer élément par élément. Mettez Vectorized à false lorsque fun n'accepte que des arguments scalaires : elle est alors évaluée point par point, ce qui est plus lent.

Waypoints indique des points d'intérêt du domaine d'intégration, comme des extrema locaux ou des discontinuités, que l'intégrateur utilise dans son maillage initial : un tableau à deux colonnes [x y] de points, ou un tableau de cellules {x y} de vecteurs de grille. L'intervalle en x est découpé aux abscisses x des points de passage, et chaque intervalle en y à leurs ordonnées y. Les points de passage doivent être réels et finis. N'utilisez pas de points de passage pour indiquer des singularités ; découpez plutôt le domaine.