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
