使用线法的 Fick 扩散第二定律：使用线法的 PDE 解决方案-matlab开发

%% 线法 D = 扩散率：菲克第二扩散定律% by Prof. Roche C. de Guzman %% 自定义外汇函数 Y1 = F(~,Y,D,nx,dx) c = Y; % 指定浓度为 Y Y1 = 零点(nx,1); % 临时衍生品对于 i = 1:nx-2 % 位置计数器Y1(i+1) = D*(c(i+2)-2*c(i+1)+c(i))/dx^2； %求解导数结尾Y1 = [Y1(2); Y1(2:nx-1); Y1(nx-1)]; % 具有零通量边界的导数结尾