function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',3);
set(ax,'DataAspectRatio',[1 0.75 1]);
set(ax,'PlotBoxAspectRatio',[5 3.3333333333333335 1]);
set(ax,'XLim',[0 10]);
set(ax,'YLim',[0 5]);
set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');

% Geometry description:
pderect([1.8 3.6000000000000001 2.6000000000000001 1.3],'R4');
pderect([6.2999999999999998 8.0999999999999996 2.6000000000000001 1.3],'R1');
pdepoly([ 0.59999999999999998,...
 5.0123456790123466,...
 9.5,...
 8.5999999999999996,...
 5.0123456790123466,...
 1.3,...
],...
[ 3,...
 4.6092592592592583,...
 3,...
 3,...
 4.3092592592592585,...
 3,...
],...
 'P1');
pderect([4.2000000000000002 5.7000000000000002 0 2.3999999999999999],'R3');
pderect([1 9 0 3],'R5');
pderect([3.5 4 3 4],'R2');
pderect([5.9000000000000004 6.4000000000000004 3 4],'R6');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','(R5-R1-R3-R4)+(P1+R2+R6)')

% Boundary conditions:
pdetool('changemode',0)
pdesetbd(42,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(41,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(40,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(39,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(38,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(37,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(36,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(34,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(32,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(30,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(28,...
'dir',...
2,...
str2mat('1','0','0','1'),...
str2mat('0','0'))
pdesetbd(27,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(26,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(25,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(23,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(21,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(19,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(17,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(16,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(14,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(13,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(11,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(8,...
'dir',...
2,...
str2mat('1','0','0','1'),...
str2mat('0','0'))
pdesetbd(7,...
'dir',...
2,...
str2mat('1','0','0','1'),...
str2mat('0','0'))
pdesetbd(6,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(5,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(4,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(3,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(2,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))
pdesetbd(1,...
'neu',...
2,...
str2mat('0','0','0','0'),...
str2mat('0','0'))

% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
pdetool('initmesh')
pdetool('jiggle')
pdetool('jiggle')

% PDE coefficients:
pdeseteq(1,...
str2mat('2*((3E10)./(2*(1+(0.2))))+(2*((3E10)./(2*(1+(0.2)))).*(0.2)./(1-(0.2)))','0','(3E10)./(2*(1+(0.2)))','0','(3E10)./(2*(1+(0.2)))','2*((3E10)./(2*(1+(0.2)))).*(0.2)./(1-(0.2))','0','(3E10)./(2*(1+(0.2)))','0','2*((3E10)./(2*(1+(0.2))))+(2*((3E10)./(2*(1+(0.2)))).*(0.2)./(1-(0.2)))'),...
str2mat('0.0','0.0','0.0','0.0'),...
str2mat('0','-9.8*2400/0.1'),...
str2mat('2400','0','0','2400'),...
'0:10',...
'0.0',...
'0.0',...
'[0 100]')
setappdata(pde_fig,'currparam',...
['3E10         ';...
'0.2          ';...
'0            ';...
'-9.8*2400/0.1';...
'2400         '])

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

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

% Solve PDE:
pdetool('solve')