Я не знал как посчитать интеграл от решения PDE, полученного с помощью PDETool на треугольной сетке и запросил помощи в форуме www.MatLab.ru. Ведущий раздела PDE В.Е. Шмелёв любезно ответил на мой вопрос.

Если Вы получили решение в PDETool, то выполните следующие действия:

1. Экспортируйте сетку (переменные p,e,t) и узловое распределение решения (переменная u) в рабочую область MATLAB.
2. Выполните следующие операторы:

   ut=pdeintrp(p,t,u);
   ar=pdetrg(p,t);
   integral_ar=sum(ar.*ut.^2)