function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',7);
pdetool('snapon','on');
set(ax,'DataAspectRatio',[1 1 1]);
set(ax,'PlotBoxAspectRatio',[1.5 1 1]);
set(ax,'XLim',[-1.5 1.5]);
set(ax,'YLim',[-1 1]);
set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');
pdetool('gridon','on');
% Geometry description:
pderect([0 0.5 -0.19999999999999996 0.80000000000000004],'R1');
pderect([0.5 1 0.80000000000000004 -0.19999999999999996],'R2');
pdecirc(0.5,0.39999999999999991,0.19999999999999996,'C1');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','(R1+R2)+
C1')
% Boundary conditions:
pdetool('changemode',0)
pdesetbd(6,...
'neu',...
1,...
'1',...
'1')
pdesetbd(5,...
'neu',...
1,...
'1',...
'1')
pdesetbd(4,...
'neu',...
1,...
'1',...
'1')
pdesetbd(3,...
'neu',...
1,...
'1',...
'1')
pdesetbd(2,...
'neu',...
1,...
'1',...
'1')
pdesetbd(1,...
'neu',...
1,...
'1',...
'1')
% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));