MATLAB平方域上的波动方程动画演示

%% Problem Definition
c = 1;
a = 0;
f = 0;
m = 1;

%% Geometry
numberOfPDE = 1;
model = createpde(numberOfPDE);
geometryFromEdges(model,@squareg);
pdegplot(model,'EdgeLabels','on'); 
ylim([-1.1 1.1]);
axis equal
title 'Geometry With Edge Labels Displayed';
xlabel x
ylabel y

%% Specify PDE Coefficients
specifyCoefficients(model,'m',m,'d',0<