diffusion with difference finite metod
显示 更早的评论
Good morning.
I do not understand the error in this script:
param.Cafeed=0.014158;% concentration [mol/dm3]
param.v=0.1;
param.L=1.1245; % led leight [dm]
param.D=0.001;
nx = 100; % to define the dimension of the initial condition
x = linspace(0,param.L,nx);
h = x(2)-x(1);
y0 = zeros(nx,1);
y0(1) = param.Cafeed;% Boundary condition
tspan = linspace(0,30,nx);
[t,y] = ode15s(@(t,y)fun(t,y,param,h),tspan,y0);%To solve the time dependence
function dy=fun(t,y,param,h)
M = length(y)+1;
dy = zeros(M+1,1);
dy(1)= 0;
for i=2:M
dy(i) = -param.v*(y(i)-y(i-1))/h+(param.D/h^2*(-2*y(i)+y(i-1)+y(i+1)));
end
thank you
1 个评论
回答（1 个）
bianchi fiammetta
2019-2-6
0 个投票
1 个评论
Torsten
2019-2-6
By the way:
This kind of equation can be solved using PDEPE.
类别
在 帮助中心 和 File Exchange 中查找有关 Eigenvalue Problems 的更多信息
2019-2-6
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
