Запустите приложение PDE Modeler и нарисуйте модульный круг и модульный квадрат.

pdecirc(0,0,1)
pderect([0 1 0 1])

Введите C1-SQ1 в поле Set formula.

Используйте значение по умолчанию граничное условие Дирихле u = 0 для всех контуров. Чтобы проверить граничное условие, переключитесь на граничный режим путем выбора Boundary> Boundary Mode. Используйте Edit> Select all, чтобы выбрать все контуры. Выберите Boundary> Specify Boundary Conditions и проверьте, что граничное условие является условием Дирихле с h = 1, r = 0.

Экспортируйте геометрию и граничные условия к рабочему пространству MATLAB путем выбирания опции Export Decomposed Geometry, Boundary Cond's в меню Boundary.

Разложите экспортируемую геометрию на минимальные области. Результатом является одна минимальная область с пятью сегментами ребра: три круговых сегмента ребра и два сегмента ребра линии.

Запишите получившуюся матрицу граничного условия в файл. Назовите файл boundary.m.

fid = wbound(b,"boundary");