File Exchange

image thumbnail

Fick's 2nd Law of Diffusion using Method of Lines

version 1.0.0 (1.98 KB) by Roche de Guzman
PDE solution using Method of Lines

1 Download

Updated 26 Apr 2019

View License

%% Method of Lines with D = diffusivity: Fick's 2nd Law of Diffusion
% by Prof. Roche C. de Guzman
%% Custom fx
function Y1 = F(~,Y,D,nx,dx)
c = Y; % assign concentration as Y
Y1 = zeros(nx,1); % temporary derivatives
for i = 1:nx-2 % position counter
Y1(i+1) = D*(c(i+2)-2*c(i+1)+c(i))/dx^2; % solve for derivatives
end
Y1 = [Y1(2); Y1(2:nx-1); Y1(nx-1)]; % derivatives with zero-flux boundary
end

Cite As

Roche de Guzman (2019). Fick's 2nd Law of Diffusion using Method of Lines (https://www.mathworks.com/matlabcentral/fileexchange/71354-fick-s-2nd-law-of-diffusion-using-method-of-lines), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (0)

MATLAB Release Compatibility
Created with R2019a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.


Learn About Live Editor