12、酸介导肿瘤生长模型的数值模拟与分析

酸介导肿瘤生长模型的数值模拟与分析

1. MOL 例程相关

1.1 稀疏矩阵积分函数

在进行 3n 个常微分方程（ODEs）的稀疏矩阵积分时， jpattern_num_1 函数会被 ode15s 调用。以下是为三个偏微分方程（PDEs）设置该函数的代码：

% Set independent, dependent variables for the calculation
% of the sparsity pattern
tbase=0;
for i=1:3*n
    ybase(i)=0.5;
end
ybase=ybase';
%
% Compute the corresponding derivative vector
ytbase=pde_1(tbase,ybase);
fac=[];
thresh=1e-16;
vectorized='on';
[Jac,fac]=numjac(@pde_1,tbase,ybase,ytbase,thresh,fac,vectorized);

特别要注意的是，在两处使用了 pde_1 来定义 ODE 雅可比矩阵的稀疏模式。

1.2 主程序

主程序 pde_1_main 与之前的主程序有相似之处，以下是其详细代码及解释：

% Three-PDE ATG model
% Clear previ