www.gusucode.com > pde 案例源码 matlab代码程序 > pde/FindTheActiveInitialConditionsExample.m

    %% Find the Active Initial Conditions
% This example shows find the active initial conditions for a region.
%%
% Create a PDE model that has a few subdomains.
model = createpde();
geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')
ylim([-1.1,1.1])
axis equal
%%
% Set initial conditions on each pair of regions.
setInitialConditions(model,12,'Face',[1,2]);
setInitialConditions(model,13,'Face',[1,3]);
setInitialConditions(model,23,'Face',[2,3]);
%%
% Check the initial conditions specification for region 1.
ics = model.InitialConditions;
ic = findInitialConditions(ics,'Face',1)