How can this function as an excel add-in ?

1 次查看（过去 30 天）
Raúl Rivera
Raúl Rivera 2023-10-23
回答： Walter Roberson 2023-10-23
Ran in:
Hello , I want to turn this code into a function , as I want to export it as a excel add-in ,with 3 entries , length, time, and concentration. The code is a pde equation with the pdepe solver , but does not work properly . What I have to change in the code to make it run? Thanks in advance for the help.
modelo_leaching(100,20,1)
Not enough input arguments.

Error in solution>modelo_leaching/pdeic (line 30)
u0 = conc; %Concentration

Error in pdepe (line 225)
temp = feval(ic,xmesh(1),varargin{:});

Error in solution>modelo_leaching (line 6)
sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t,conc);
function a=modelo_leaching(length,tiempo,conc)
x = 0:1:length; %length
t = 0:1:tiempo;%tiempo
m = 0;
sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t,conc);
u1 = sol(:,:,1);
surf(x,t,u1)
zlabel('Mineral conc(x,t)')
xlabel('Distance x')
ylabel('Time t')
function [c,f,s] = pdefun(length,t,u,dudx) % Equation to solve
U=10;
E=0.1e-8;
c = 1;
k=0.0017;
j=1.3703;
b=1.071;
cc=0.1195;
Ox=8;
cn=159;
f = E.* dudx;% - U.*u;
y =k*((u-0.17)^j)*((cn)^b)*((Ox)^cc) ;% ((u-0.17)^1); %(u-0.17)^0.3;
s = -y -U.*dudx ;
end
% ---------------------------------------------
function u0 = pdeic(length,conc) % Initial Conditions
u0 = conc; %Concentration
end
% ---------------------------------------------
function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,tiempo,conc) % Boundary Conditions
Cai=conc;%Concentration
pl = ul-Cai;
ql = 0;
pr = 0;
qr = 1;
end
end

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2023-10-23
Ran in:
modelo_leaching(100,20,1)
function a=modelo_leaching(length,tiempo,conc)
x = 0:1:length; %length
t = 0:1:tiempo;%tiempo
m = 0;
sol = pdepe(m,@pdefun,@(X)pdeic(X,conc),@(xl,ul,xr,ur,tiempo)pdebc(xl,ul,xr,ur,tiempo,conc),x,t);
u1 = sol(:,:,1);
surf(x,t,u1)
zlabel('Mineral conc(x,t)')
xlabel('Distance x')
ylabel('Time t')
function [c,f,s] = pdefun(length,t,u,dudx) % Equation to solve
U=10;
E=0.1e-8;
c = 1;
k=0.0017;
j=1.3703;
b=1.071;
cc=0.1195;
Ox=8;
cn=159;
f = E.* dudx;% - U.*u;
y =k*((u-0.17)^j)*((cn)^b)*((Ox)^cc) ;% ((u-0.17)^1); %(u-0.17)^0.3;
s = -y -U.*dudx ;
end
% ---------------------------------------------
function u0 = pdeic(length,conc) % Initial Conditions
u0 = conc; %Concentration
end
% ---------------------------------------------
function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,tiempo,conc) % Boundary Conditions
Cai=conc;%Concentration
pl = ul-Cai;
ql = 0;
pr = 0;
qr = 1;
end
end

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Introduction to Installation and Licensing 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by