% FEMLAB Model M-file
% Generated 28-Jan-2002 14:30:03 by FEMLAB 2.1.0.127.

flclear fem
% FEMLAB Version
clear vrsn;
vrsn.name='FEMLAB 2.1';
vrsn.major=0;
vrsn.build=127;
fem.version=vrsn;

% Space dimensions
fem.sdim={'x','y','z'};

% Geometry
p=[-0.017500000000000002 -0.017500000000000002 -0.01 -0.01 -0.01 -0.01  ...
-0.0050000000000000001 -0.0050000000000000001 0.0050000000000000001  ...
0.0050000000000000001 0.01 0.01 0.01 0.01 0.017500000000000002  ...
0.017500000000000002;-0.0050000000000000001 0.0050000000000000001 -0.01  ...
-0.0050000000000000001 0.0050000000000000001 0.01 -0.0050000000000000001  ...
0.0050000000000000001 -0.0050000000000000001 0.0050000000000000001 -0.01  ...
-0.0050000000000000001 0.0050000000000000001 0.01 -0.0050000000000000001  ...
0.0050000000000000001];
rb={1:16,[1 1 2 3 3 5 6 7 7 8 9 11 12 13 13 15;2 4 5 4 11 6 14 8 9 10 10 12  ...
15 14 16 16],zeros(3,0),zeros(4,0)};
wt={zeros(1,0),ones(2,16),zeros(3,0),zeros(4,0)};
lr={[NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN],[0 1  ...
0 0 1 0 0 1 0 1 0 1 1 1 0 1;1 0 1 1 0 1 1 0 1 0 1 0 0 0 1 0],zeros(2,0), ...
zeros(2,0)};
tmp=solid2(p,rb,wt,lr);
g1=extrude(tmp,'Distance',0.14999999999999999,'Scale',[1;1],'Displ',[0;0], ...
'Wrkpln',[0 0 0;0 1 0;0 0 1]);

clear s f c p
objs={g1};
name={'EXT1'};
s.objs=objs;
s.name=name;

objs={};
name={};
f.objs=objs;
f.name=name;

objs={};
name={};
c.objs=objs;
c.name=name;

objs={};
name={};
p.objs=objs;
p.name=name;

drawstruct=struct('s',s,'f',f,'c',c,'p',p);
fem.draw=drawstruct;
fem.geom=geomcsg(fem);

% Recorded command sequence

clear appl

% Application mode 1
appl{1}.mode='flpdeht3d(''dim'',{''T''},''sdim'',{''x'',''y'',''z''},''submode'',''std'',''tdiff'',''on'')';
appl{1}.dim={'T'};
appl{1}.form='coefficient';
appl{1}.border='off';
appl{1}.name='ht';
appl{1}.var={};
appl{1}.assign={'C';'C';'Const';'Const';'Q';'Q';'Tamb';'Tamb';'Tinf';'Tinf'; ...
'flux';'flux';'flux_x';'flux_x';'flux_y';'flux_y';'flux_z';'flux_z';'gradT'; ...
'gradT';'gradTx';'gradTx';'gradTy';'gradTy';'gradTz';'gradTz';'h';'h'; ...
'n_flux';'n_flux';'q';'q';'rho';'rho'};
appl{1}.equ.rho={'7800'};
appl{1}.equ.C={'460'};
appl{1}.equ.k={{{'45'}}};
appl{1}.equ.Q={'0'};
appl{1}.equ.ind=1;
appl{1}.usage=1;
appl{1}.bnd.q={'0','0','0','0','0','0'};
appl{1}.bnd.h={'0','0','5','10','0','0'};
appl{1}.bnd.Tinf={'0','0','T0','T0','0','0'};
appl{1}.bnd.Const={'0','0','0','0','0','0'};
appl{1}.bnd.Tamb={'0','0','0','0','0','0'};
appl{1}.bnd.T={'(273+60)','0','0','0','(273+60)-10/15e-2*x','(273+50)'};
appl{1}.bnd.type={'T','q0','q','q','T','T'};
appl{1}.bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 6];
appl{1}.init.sd={{{'T0'}}};
appl{1}.init.ind=1;

fem.appl=appl;

% Initialize mesh
fem.mesh=meshinit(fem,...
	'Out',    {'mesh'},...
	'jiggle', 'on',...
	'Hcurve', 0.40000000000000002,...
	'Hcutoff',0.01,...
	'Hgrad',  1.3999999999999999,...
	'Hmaxfact',1);

% Dimension
fem.dim={{'T'}};

% Boundary conditions
fem.border=0;

% Usage
fem.usage=1;

% Problem form
fem.form='coefficient';

% Differentiation simplification
fem.simplify='on';

% Differentiation rules
fem.rules={};

% Define variables
fem.variables={'T0', 293};

% Define application mode variables
fem.var={};

% Boundary conditions
clear bnd
bnd.q={{{'0'}},{{'0'}},{{'5'}},{{'10'}},{{'0'}}};
bnd.g={{{'0'}},{{'0'}},{{'5.*T0'}},{{'10.*T0'}},{{'0'}}};
bnd.h={{{'1'}},{{'0'}},{{'0'}},{{'0'}},{{'1'}}};
bnd.r={{{'0'}},{{'0'}},{{'0'}},{{'0'}},{{'333-200/3*x'}}};
bnd.var={'q',{'0','0','0','0','0'},'h',{'0','0','5','10','0'},'Tinf',{'0', ...
'0','T0','T0','0'},'Const',{'0','0','0','0','0'},'Tamb',{'0','0','0','0', ...
'0'},'n_flux','-ncu1'};
bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 1];
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'3588000'}}};
equ.c={{{'45'}}};
equ.al={{{'0';'0';'0'}}};
equ.ga={{{'0';'0';'0'}}};
equ.be={{{'0';'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'rho',{'7800'},'C',{'460'},'Q',{'0'},'gradTx','-Tx','gradTy','-Ty', ...
'gradTz','-Tz','gradT','sqrt(Tx.^2+Ty.^2+Tz.^2)','flux_x','-cu1x','flux_y', ...
'-cu1y','flux_z','-cu1z','flux','sqrt(cu1x.^2+cu1y.^2+cu1z.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'T0'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:0.10000000000000001:1,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

% Save current mesh
restart.mesh=fem.mesh;
restart.usage=fem.usage;

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',0,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'geom',   'on',...
	'geomcol','bginv',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=1  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Dimension
fem.dim={{'T'}};

% Boundary conditions
fem.border=0;

% Usage
fem.usage=1;

% Problem form
fem.form='coefficient';

% Differentiation simplification
fem.simplify='on';

% Differentiation rules
fem.rules={};

% Define variables
fem.variables={'T0', 293};

% Define application mode variables
fem.var={};

% Boundary conditions
clear bnd
bnd.q={{{'0'}},{{'0'}},{{'5'}},{{'10'}},{{'0'}},{{'0'}}};
bnd.g={{{'0'}},{{'0'}},{{'5.*T0'}},{{'10.*T0'}},{{'0'}},{{'0'}}};
bnd.h={{{'1'}},{{'0'}},{{'0'}},{{'0'}},{{'1'}},{{'1'}}};
bnd.r={{{'333'}},{{'0'}},{{'0'}},{{'0'}},{{'333-200/3*x'}},{{'323'}}};
bnd.var={'q',{'0','0','0','0','0','0'},'h',{'0','0','5','10','0','0'}, ...
'Tinf',{'0','0','T0','T0','0','0'},'Const',{'0','0','0','0','0','0'},'Tamb', ...
{'0','0','0','0','0','0'},'n_flux','-ncu1'};
bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 6];
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'3588000'}}};
equ.c={{{'45'}}};
equ.al={{{'0';'0';'0'}}};
equ.ga={{{'0';'0';'0'}}};
equ.be={{{'0';'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'rho',{'7800'},'C',{'460'},'Q',{'0'},'gradTx','-Tx','gradTy','-Ty', ...
'gradTz','-Tz','gradT','sqrt(Tx.^2+Ty.^2+Tz.^2)','flux_x','-cu1x','flux_y', ...
'-cu1y','flux_z','-cu1z','flux','sqrt(cu1x.^2+cu1y.^2+cu1z.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'T0'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:0.10000000000000001:1,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

% Save current mesh
restart.mesh=fem.mesh;
restart.usage=fem.usage;

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',0,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'geom',   'on',...
	'geomcol','bginv',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=1  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Dimension
fem.dim={{'T'}};

% Boundary conditions
fem.border=0;

% Usage
fem.usage=1;

% Problem form
fem.form='coefficient';

% Differentiation simplification
fem.simplify='on';

% Differentiation rules
fem.rules={};

% Define variables
fem.variables={'T0', 293};

% Define application mode variables
fem.var={};

% Boundary conditions
clear bnd
bnd.q={{{'0'}},{{'0'}},{{'5'}},{{'10'}},{{'0'}},{{'0'}}};
bnd.g={{{'0'}},{{'0'}},{{'5.*T0'}},{{'10.*T0'}},{{'0'}},{{'0'}}};
bnd.h={{{'1'}},{{'0'}},{{'0'}},{{'0'}},{{'1'}},{{'1'}}};
bnd.r={{{'333'}},{{'0'}},{{'0'}},{{'0'}},{{'333-200/3*x'}},{{'323'}}};
bnd.var={'q',{'0','0','0','0','0','0'},'h',{'0','0','5','10','0','0'}, ...
'Tinf',{'0','0','T0','T0','0','0'},'Const',{'0','0','0','0','0','0'},'Tamb', ...
{'0','0','0','0','0','0'},'n_flux','-ncu1'};
bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 6];
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'3588000'}}};
equ.c={{{'45'}}};
equ.al={{{'0';'0';'0'}}};
equ.ga={{{'0';'0';'0'}}};
equ.be={{{'0';'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'rho',{'7800'},'C',{'460'},'Q',{'0'},'gradTx','-Tx','gradTy','-Ty', ...
'gradTz','-Tz','gradT','sqrt(Tx.^2+Ty.^2+Tz.^2)','flux_x','-cu1x','flux_y', ...
'-cu1y','flux_z','-cu1z','flux','sqrt(cu1x.^2+cu1y.^2+cu1z.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'T0'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:0.10000000000000001:1,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

% Save current mesh
restart.mesh=fem.mesh;
restart.usage=fem.usage;

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',0,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'geom',   'on',...
	'geomcol','bginv',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=1  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Dimension
fem.dim={{'T'}};

% Boundary conditions
fem.border=0;

% Usage
fem.usage=1;

% Problem form
fem.form='coefficient';

% Differentiation simplification
fem.simplify='on';

% Differentiation rules
fem.rules={};

% Define variables
fem.variables={'T0', 293};

% Define application mode variables
fem.var={};

% Boundary conditions
clear bnd
bnd.q={{{'0'}},{{'0'}},{{'5'}},{{'10'}},{{'0'}},{{'0'}}};
bnd.g={{{'0'}},{{'0'}},{{'5.*T0'}},{{'10.*T0'}},{{'0'}},{{'0'}}};
bnd.h={{{'1'}},{{'0'}},{{'0'}},{{'0'}},{{'1'}},{{'1'}}};
bnd.r={{{'333'}},{{'0'}},{{'0'}},{{'0'}},{{'333-200/3*x'}},{{'323'}}};
bnd.var={'q',{'0','0','0','0','0','0'},'h',{'0','0','5','10','0','0'}, ...
'Tinf',{'0','0','T0','T0','0','0'},'Const',{'0','0','0','0','0','0'},'Tamb', ...
{'0','0','0','0','0','0'},'n_flux','-ncu1'};
bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 6];
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'3588000'}}};
equ.c={{{'45'}}};
equ.al={{{'0';'0';'0'}}};
equ.ga={{{'0';'0';'0'}}};
equ.be={{{'0';'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'rho',{'7800'},'C',{'460'},'Q',{'0'},'gradTx','-Tx','gradTy','-Ty', ...
'gradTz','-Tz','gradT','sqrt(Tx.^2+Ty.^2+Tz.^2)','flux_x','-cu1x','flux_y', ...
'-cu1y','flux_z','-cu1z','flux','sqrt(cu1x.^2+cu1y.^2+cu1z.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'T0'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:0.10000000000000001:1,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

% Save current mesh
restart.mesh=fem.mesh;
restart.usage=fem.usage;

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',0,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'geom',   'on',...
	'geomcol','bginv',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=1  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Dimension
fem.dim={{'T'}};

% Boundary conditions
fem.border=0;

% Usage
fem.usage=1;

% Problem form
fem.form='coefficient';

% Differentiation
fem.diff='off';

% Differentiation simplification
fem.simplify='on';

% Differentiation rules
fem.rules={};

% Define variables
fem.variables={'T0', 293};

% Define application mode variables
fem.var={};

% Boundary conditions
clear bnd
bnd.q={{{'0'}},{{'0'}},{{'5'}},{{'10'}},{{'0'}},{{'0'}}};
bnd.g={{{'0'}},{{'0'}},{{'5.*T0'}},{{'10.*T0'}},{{'0'}},{{'0'}}};
bnd.h={{{'1'}},{{'0'}},{{'0'}},{{'0'}},{{'1'}},{{'1'}}};
bnd.r={{{'333'}},{{'0'}},{{'0'}},{{'0'}},{{'333-200/3*x'}},{{'323'}}};
bnd.var={'q',{'0','0','0','0','0','0'},'h',{'0','0','5','10','0','0'}, ...
'Tinf',{'0','0','T0','T0','0','0'},'Const',{'0','0','0','0','0','0'},'Tamb', ...
{'0','0','0','0','0','0'},'n_flux','-ncu1'};
bnd.ind=[1 2 3 4 3 3 4 4 5 5 5 5 3 3 4 4 4 6];
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'3588000'}}};
equ.c={{{'45'}}};
equ.al={{{'0';'0';'0'}}};
equ.ga={{{'0';'0';'0'}}};
equ.be={{{'0';'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'rho',{'7800'},'C',{'460'},'Q',{'0'},'gradTx','-Tx','gradTy','-Ty', ...
'gradTz','-Tz','gradT','sqrt(Tx.^2+Ty.^2+Tz.^2)','flux_x','-cu1x','flux_y', ...
'-cu1y','flux_z','-cu1z','flux','sqrt(cu1x.^2+cu1y.^2+cu1z.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'T0'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:10,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

% Save current mesh
restart.mesh=fem.mesh;
restart.usage=fem.usage;

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',0,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'geom',   'on',...
	'geomcol','bginv',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=10  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',7,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=3  Slice: temperature (T)  ',...
	'solnum', 4,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off')

% Animate solution
postmovie(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',7,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off',...
	'movie',  'on',...
	'repeat', 1,...
	'fps',    12,...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'renderer','zbuffer',...
	'title',  'Slice: temperature (T)  ',...
	'scenelight','off',...
	'camlight','off',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998])

% Plot solution
postplot(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',7,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'view',   [-37.5 30.000000000000004],...
	'title',  'Time=10  Slice: temperature (T)  ',...
	'solnum', 11,...
	'renderer','zbuffer',...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off',...
	'scenelightpos',[-0.33760150336373573 -0.53771274805604652 0.39131189606246325])

% Animate solution
postmovie(fem,...
	'context','local',...
	'slicedata',{'T','cont','on'},...
	'slicexspacing',5,...
	'sliceyspacing',0,...
	'slicezspacing',7,...
	'slicefacestyle','interp',...
	'sliceedgestyle','none',...
	'slicetag','SlicePatch',...
	'slicemap','jet',...
	'slicemaxmin','off',...
	'slicebar','on',...
	'slicebartag','SliceBar',...
	'lightmodel','phong',...
	'lightreflection','default',...
	'transparency',1,...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998],...
	'scenelight','off',...
	'camlight','off',...
	'movie',  'on',...
	'repeat', 1,...
	'fps',    12,...
	'axisequal','on',...
	'axis',   [0 0.14999999999999999 -0.02 0.02 -0.01 0.01],...
	'axisvisible','on',...
	'grid',   'on',...
	'renderer','zbuffer',...
	'title',  'Slice: temperature (T)  ',...
	'scenelight','off',...
	'camlight','off',...
	'campos', [-0.33760150336373573 -0.53771274805604652 0.39131189606246325],...
	'camprojection','orthographic',...
	'camtarget',[0.074999999999999997 0 0],...
	'camup',  [0 0 1],...
	'camva',  14.613108985511817,...
	'axispos',[-0.024999999999999967 -0.05299999999999995 1.085 1.1409999999999998])
