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)