%% 线性法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)-2c(i+1)+c(i))/dx^2; % 求解导数结尾Y1 = [Y1(2); Y1(2:nx-1); Y1(nx-1)]; % 具有零通量边界的导数结尾