IM unsure of what this error means.

Error using symfun/privResolveArgs (line 223)

Symbolic function input arguments must match.

Error in sym/privBinaryOp (line 1031)

args = privResolveArgs(A, B);

Error in + (line 7)

X = privBinaryOp(A, B, 'symobj::zipWithImplicitExpansion', '_plus');

Error in ThermalConductionProject (line 222)

THETA = theta_2c+theta_2a+theta_2b + theta_2d +theta_2e+theta_3a+theta_3b+theta_3c+theta_3d+theta_3e

syms f g x y z

% % % % Used To determine Values of Mew % % % % %

% clc

% clear

% syms mew

% Bi =0.0837

% H=0.0015;

% l=0.005;

% syms f(x)

% f(x) =((mew*mew*l)*((tan(mew*l*.5)))+Bi);

% sol = vpasolve(f)

% % % % % Variables % % % % %

k = (0.6569);

h = 11;

L=(.005);

W=(.0015);

H=(.0015);

X=(.051/4);

Y=(.01/2);

Z=(.005/4);

alpha = 2.57*10^(-7);

theta_2b=0;

% % % % % Temperature Formula 2b % % % % %

for n=0:100

t=0;

beta1 = (n*pi)/H;

eta1 = (((2*pi)/L)*(n+0.5));

mew1 = -18.84;

lamda1 = (beta1^2)+(eta1^2)+(mew1^2);

f(x,y,z) = sin(beta1*z)*cos(mew1*y)*sin(mew1*y)*sin(eta1*x);

C_nmp_numerator1 = int(int(int(f,z,[0 H]),y,[0 W]),x,[0 (L/2)]);

g(x,y,z) = ((cos(mew1*y)^2)*sin(mew1*y) - ((h/(k*mew1))*(sin(mew1*y)^2)*cos(mew1*y)));

C_nmp_denominator1 = ((L*H)/8)* (int(g,y,[0 W]));

C_nmp1 = C_nmp_numerator1/C_nmp_denominator1;

theta_2b = theta_2b+(C_nmp1*(sin(beta1*z)*sin(eta1*x)*((cos(mew1*y)-(h/(k*mew1))*sin(mew1*y)))*exp(-(lamda1^2)*alpha*t)));

end

theta_2b = vpa(subs(theta_2b, {x,y,z}, {X,Y,Z}));

theta_2a=0;

% % % % % Temperature Formula 2a % % % % %

for n=0:100

t=0;

beta2 = (n*pi)/H;

eta2 = (n*pi)/W;

mew2 = -18.84;

lamda2= (beta2^2)+(eta2^2)+(mew2^2);

f(x,y,z) = sin(beta2*z)*cos(mew2*x)*sin(eta2*y)*sin(mew2*x);

C_nmp_numerator2 = int(int(int(f,z,[0 H]),y,[0 W]),x,[0 (L/2)]);

g(x,y,z) = ((cos(mew2*x)^2)*sin(mew2*x) - ((h/(k*mew2))*(sin(mew2*x)^2)*cos(mew2*x)));

C_nmp_denominator2 = ((H*W)/4)* (int(g,x,[0 (L/2)]));

C_nmp2 = C_nmp_numerator2/C_nmp_denominator2;

theta_2a = theta_2a+(C_nmp2*(sin(beta2*z)*sin(eta2*y)*((cos(mew2*x)-(h/(k*mew2))*sin(mew2*x)))*exp(-(lamda2^2)*alpha*t)));

end

theta_2a = vpa(subs(theta_2a, {x,y,z}, {X,Y,Z}));

theta_2d=0;

% % % % % % Temperature Formula 2d % % % % %

for n=0:100

t=0;

eta3= (n*pi)/W;

beta3 = (((2*pi)/L)*(n+0.5));

mew3 = -1036;

lamda3 = (beta3^2)+(eta3^2)+(mew3^2);

f(x,y,z) = sin(beta3*x)*cos(mew3*z)*sin(eta3*y)*sin(mew3*z);

C_nmp_numerator3 = int(int(int(f,z,[0 H]),y,[0 W]),x,[0 (L/2)]);

g(x,y,z) = ((cos(mew3*z)^2)*sin(mew3*z) - ((h/(k*mew3))*(sin(mew3*z)^2)*cos(mew3*z)));

C_nmp_denominator3 = ((L*W)/8)* (int(g,z,[0 H]));

C_nmp3 = C_nmp_numerator3/C_nmp_denominator3;

theta_2d = theta_2d+(C_nmp3*(sin(beta3*x)*exp(-(lamda3^2)*alpha*t)*sin(eta3*y)*((cos(mew3*z)-(h/(k*mew3))*sin(mew3*z)))));

end

theta_2d = vpa(subs(theta_2d, {x,y,z}, {X,Y,Z}));

theta_2e=0;

% % % % % Temperature Formula 2e % % % % %

for n=0:100

t=0;

beta4 = (((2*pi)/L)*(n+0.5));

eta4 = (n*pi)/W;

mew4 = 0;

lamda4 = (beta4^2)+(eta4^2)+(mew4^2);

f(x,y,z) = sin(beta4*x)*sin(mew4*z)*sin(eta4*y);

C_nmp_numerator4 = int(int(int(f,z,[0 H]), y,[0 W]), x, [0 (L/2)]);

C_nmp_denominator4 = ((H*L*W)/16);

C_nmp4 = C_nmp_numerator4/C_nmp_denominator4;

theta_2e = theta_2e+(C_nmp4 * (sin(beta4*x)*sin(mew4*z)*sin(eta4*y)*exp(-(lamda4^2)*alpha*t)));

end

theta_2e = vpa(subs(theta_2e, {x,y,z}, {X,Y,Z}));

theta_3a=0;

% % % % % Temperature Formula 3a % % % % %

for n=0:100

t=0;

beta5 = (n*pi)/H;

eta5 = (n*pi)/W;

mew5 = -81;

lamda5 = (beta5^2)+(eta5^2)+(mew5^2);

f(x,y,z) = sin(beta5*z)*cos(mew5*x)*sin(eta5*y);

C_nmp_numerator5 = int(int(int(f,z,[0 H]), y,[0 W]), x, [0 (L/2)]);

C_nmp_denominator5 = ((H*L*W)/16);

C_nmp5 = C_nmp_numerator5/C_nmp_denominator5;

theta_3a = theta_3a+(C_nmp5 * (sin(beta5*z)*cos(mew5*x)*sin(eta5*y) * exp(-(lamda5^2)*alpha*t)));

end

theta_3a = vpa(subs(theta_3a, {x,y,z}, {X,Y,Z}));

theta_3b=0;

% % % % % Temperature Formula 3b % % % % %

for n=0:100

t=0;

beta6 = (n*pi)/H;

eta6 = (((2*pi)/L)*(n+0.5));

mew6 = -1036;

lamda6 = (beta6^2)+(eta6^2)+(mew6^2);

f(x,y,z) = sin(beta6*z)*cos(eta6*x)*sin(mew6*y)*cos(mew6*y);

C_nmp_numerator6 = int(int(int(f,z,[0 H]),y,[0 W]),x,[0 (L/2)]);

g(x,y,z) = ((cos(mew6*y)^2)*sin(mew6*y) - ((h/(k*mew6))*(sin(mew6*y)^2)*cos(mew6*y)));

C_nmp_denominator6 = ((L*H)/8)* (int(g,y,[0 W]));

C_nmp6 = C_nmp_numerator6/C_nmp_denominator6;

theta_3b = theta_3b+(C_nmp6*(sin(beta6*z)*cos(eta6*x)*((cos(mew6*y)-(h/(k*mew6))*sin(mew6*y)))*exp(-(lamda6^2)*alpha*t)));

end

theta_3b = vpa(subs(theta_3b, {x,y,z}, {X,Y,Z}));

theta_3d=0;

% % % % % Temperature Formula 3d % % % % %

for n=0:100

t=0;

beta7 = (n*pi)/W;

eta7 = (((2*pi)/L)*(n+0.5));

mew7 = -303;

lamda7 = (beta7^2)+(eta7^2)+(mew7^2);

f(x,y,z) = sin(beta7*y)*cos(eta7*x)*sin(mew7*z)*cos(mew7*z);

C_nmp_numerator7 = int(int(int(f,z,[0 H]),y,[0 W]),x,[0 (L/2)]);

g(x,y,z) = ((cos(mew7*z)^2)*sin(mew7*z) - ((h/(k*mew7))*(sin(mew7*z)^2)*cos(mew7*z)));

C_nmp_denominator7 = ((L*W)/8)* (int(g,z,[0 H]));

C_nmp7 = C_nmp_numerator7/C_nmp_denominator7;

theta_3d = theta_3d+(C_nmp7*(sin(beta7*y)*cos(eta7*x)*(cos(mew7*z)-(h/(k*mew7))*sin(mew7*z))*exp(-(lamda7^2)*alpha*t)));

end

theta_3d = vpa(subs(theta_3d, {x,y,z}, {X,Y,Z}));

theta_3e=0;

% % % % % Temperature Formula 3e % % % % %

for n=0:100

t=0;

beta8 = (n*pi)/W;

eta8 = (((2*pi)/L)*(n+0.5));

mew8 = 57;

lamda8 = (beta8^2)+(eta8^2)+(mew8^2);

f(x,y,z) = sin(beta8*y)*cos(eta8*x)*sin(mew8*z);

C_nmp_numerator8 = int(int(int(f,z,[0 H]), y,[0 W]), x, [0 (L/2)]);

C_nmp_denominator8 = ((H*L*W)/16);

C_nmp8 = C_nmp_numerator8/C_nmp_denominator8;

theta_3e = theta_3e+(C_nmp8 * (sin(beta8*y)*cos(eta8*x)*sin(mew8*z) * exp(-(lamda8^2)*alpha*t)));

end

theta_3e = vpa(subs(theta_3e, {x,y,z}, {X,Y,Z}));

theta_3c=0;

% % % % % Temperature Formula 3-c % % % % %

for n=0:100

t=0;

beta9 = (n*pi)/H;

eta9 = (((2*pi)/L)*(n+0.5));

mew9 = -3145;

lamda9 = (beta9^2)+(eta9^2)+(mew9^2);

f(x,y,z) = sin(beta9*z)*cos(eta9*x)*sin(mew9*y);

C_nmp_numerator9 = int(int(int(f,z,[0 H]), y,[0 W]), x, [0 (L/2)]);

C_nmp_denominator9 = ((H*L*W)/16);

C_nmp9 = C_nmp_numerator9/C_nmp_denominator9;

theta_3c = theta_3c+(C_nmp9 * (sin(beta9*z)*cos(eta9*x)*sin(mew9*y) * exp(-(lamda9^2)*alpha*t)));

end

theta_3c = vpa(subs(theta_3c, {x,y,z}, {X,Y,Z}));

theta_2c=0;

% % % % % Temperature Formula 2c % % % % %

for n=0:100

t=0;

eta10 = (((2*pi)/L)*(n+0.5));

beta10 = (n*pi)/H;

mew10 = 0;

lamda10 = (beta10^2)+(eta10^2)+(mew10^2);

f(x,y,z) = (sin(beta10*z)*sin(eta10*x)*sin(mew10*y));

C_nmp_numerator10 = int(int(int(f,z,[0 H]), y,[0 W]), x, [0 (L/2)]);

C_nmp_denominator10 = ((H*L*W)/16);

C_nmp10 = C_nmp_numerator10/C_nmp_denominator10;

theta_2c = theta_2c+(C_nmp10 * (sin(beta10*z)*sin(eta10*x)*sin(mew10*y)*exp(-(lamda10^2)*alpha*t)));

end

theta_2c = vpa(subs(theta_2c, {x,y,z}, {X,Y,Z}));

THETA = theta_2c+theta_2a+theta_2b + theta_2d +theta_2e+theta_3a+theta_3b+theta_3c+theta_3d+theta_3e

[clayton holmes]

When each temperature formulation is run seperatly I do not get an error, but when i combine so i can add themm all up that is when i have an issue.