使用直接刚度有限元法求解Fick\'s 2nd Law of Diffusion以获得给定一维空间x和时间t下的浓度c，开发环境为MATLAB

D =扩散率的%%有限元方法应用于菲克第二扩散定律：dc / dt = D * d ^ 2c / dx ^ 2 % by Prof. Roche C. de Guzman 清除; clc; 关闭所有'）; %% 给定xi = 0; xf = 0.6; dx = 0.04; % x 范围和步长 = dx [m] xL = 0; xU = 0.1; % 初始值 x 下限和上限 [m] ti = 0; tf = 0.05; dt = 4e-4; % t 范围和步长 = dt [s] ci = 2; % 初始浓度值 [ng/L] cLU = 8; x 上下限范围内的 % 初始浓度值 [ng/L] D = 1.5； % 扩散率或扩散系数 [m^2/s] %% 计算％自变量：x和t X = xi:dx:xf; nx = 数字（X）； T = ti:dt:tf; nt = numel(T); % x 和