matlab pdetool for diffusion

我积极编写的matlab M文件用来求解菲克第二定律的扩散

function pdemodel
[pde_figax]=pdeinit;
pdetool（‘appl_cb‘1）;
set（ax‘DataAspectRatio‘[1 0.029999999999999995 1]）;
set（ax‘PlotBoxAspectRatio‘[83.333333333333329 33.333333333333336 166666.66666666666]）;
set（ax‘XLimMode‘‘auto‘）;
set（ax‘YLimMode‘‘auto‘）;
set（ax‘XTickMode‘‘auto‘）;
set（ax‘YTickMode‘‘auto‘）;

% Geometry description:
pderect（[0 0.001 1.0000000000000001e-005 0]‘R1‘）;
set（findobj（get（pde_fig‘Children‘）‘Tag‘‘PDEEval‘）‘String‘‘R1‘）

% Boundary conditions:
pdetool（‘changemode‘0）
pdesetbd（4...
‘dir‘...
1...
‘1‘...
‘0.001‘）
pdesetbd（3...
‘neu‘...
1...
‘0‘...
‘0‘）
pdesetbd（2...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（1...
‘neu‘...
1...
‘0‘...
‘0‘）

% Mesh generation:
setappdata（pde_fig‘Hgrad‘1.3）;
setappdata（pde_fig‘refinemethod‘‘regular‘）;
pdetool（‘initmesh‘）

% PDE coefficients:
pdeseteq（2...
‘0.000000016‘...
‘0.0‘...
‘0‘...
‘1.0‘...
‘0:100‘...
‘0‘...
‘0.0‘...
‘[0 100]‘）
setappdata（pde_fig‘currparam‘...
[‘0.000000016‘;...
‘0.0        ‘;...
‘0          ‘;...
‘1.0        ‘]）

% Solve parameters:
setappdata（pde_fig‘solveparam‘...
str2mat（‘0‘‘1000‘‘10‘‘pdeadworst‘...
‘0.5‘‘longest‘‘0‘‘1E-4‘‘‘‘fixed‘‘Inf‘））

% Plotflags and user data strings:
setappdata（pde_fig‘plotflags‘[1 1 1 1 1 1 1 1 0 0 0 101 1 0 1 0 0 1]）;
setappdata（pde_fig‘colstring‘‘‘）;
setappdata（pde_fig‘arrowstring‘‘‘）;
setappdata（pde_fig‘deformstring‘‘‘）;
setappdata（pde_fig‘heightstring‘‘‘）;

% Solve PDE:
pdetool（‘solve‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1471  2011-05-23 16:51  for fick law.m

----------- ---------  ---------- -----  ----

                 1471                    1