function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',1);
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:
pdecirc(-0.5,-0.20,0.40,'C1');
pderect([-1 1 0 -0.40],'R1');
pderect([0.5 1 0.80 -0.60],'R2');
pdecirc(0.5,0.20,0.40,'C2');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','(C1+R1+R2)-C2')

% Boundary conditions:
pdetool('changemode',0)
pdetool('removeb',[8 10 11 13 16 20 23 ]);
pdesetbd(20,'neu',1,'0','-5')
pdesetbd(19,'neu',1,'0','-5')
pdesetbd(18,'neu',1,'0','-5')
pdesetbd(17,'neu',1,'0','-5')
pdesetbd(16,'neu',1,'0','-5')
pdesetbd(15,'neu',1,'0','-5')
pdesetbd(14,'neu',1,'0','-5')
pdesetbd(13,'neu',1,'0','-5')
pdesetbd(12,'dir',1,'1','0')
pdesetbd(11,'dir',1,'1','0')
pdesetbd(10,'dir',1,'1','0')
pdesetbd(9,'dir',1,'1','0')
pdesetbd(8,'dir',1,'1','0')
pdesetbd(7,'dir',1,'1','0')
pdesetbd(6,'dir',1,'1','0')
pdesetbd(5,'dir',1,'1','0')
pdesetbd(4,'dir',1,'1','0')
pdesetbd(3,'dir',1,'1','0')
pdesetbd(2,'dir',1,'1','0')
pdesetbd(1,'dir',1,'1','0')

% Mesh generation:
setuprop(pde_fig,'Hgarad',1.3);
setuprop(pde_fig,'refinemethod','regular');
pdetool('initmesh')
pdetool('refine')
pdetool('refine')

%PDE coefficients:
pdeseteq(1,'1.0','0.0','10.0','1.0','0.10','0.0','0.0','[0 100]')
setuprop(pde_fig,'currparam',['1.0';'0.0';'10.0';'1.0'])

% Solve parameters:
setuprop(pde_fig,'solveparam',...
str2mat('0','2880','10','pdeadworst','0.5','longest','0','1E-4','','fixed','Inf'))


%Plotflags and user data strings :
setuprop(pde_fig,'plotflags',[1 1 1 1 1 1 1 1 0 0 0 1 1 0 1 0 0 1]);
setuprop(pde_fig,'colstring','');
setuprop(pde_fig,'arrowstring','');
setuprop(pde_fig,'deformstring','');
setuprop(pde_fig,'heightstring','');

% Solve PDE:
pdetool('solve')