pde

model = createpde();
%Геометрия задачи
polygon = polyshape([0 0 1 1],[1 0 0 1]);
pol = triangulation(polygon);
nodes = pol.Points';
elements = pol.ConnectivityList';
geometryFromMesh(model,nodes,elements);
%Строим сетку, Hmax - макс. размер элемента
generateMesh(model,'Hmax', 0.025);
figure
pdemesh(model)
%Задаем граничные условия
applyBoundaryCondition(model,"dirichlet","Edge",[1 2 3 4], 'u', 0);
%Задаем коэффициенты уравнения
specifyCoefficients(model, 'm', 0, 'd', 0, 'c', 1, 'a', 0, 'f', 1);
%Решаем уравнение
result = solvepde(model);
%Визуализация решения
u = result.NodalSolution;
pdeplot(model,'XYData',u,'ZData',u,'Mesh','on')
xlabel('x')
ylabel('y')
zlabel('w')
%В виде градиентов
wx = result.XGradients;
wy = result.YGradients;
pdeplot(model, 'FlowData',[wx wy])
xlabel('x')
ylabel('y')