[求助]下面的程序错在哪里？

gina-yanke

zuhe.m文件：
function equal=zuhe(l,n,d)
h=log(n)/log(d);
s=0;
for q=0:h-1;
t(q+1)=1;
for p=0:n/d^q-1;
t(q+1)=t(q+1)*(1-l/(n-p));
end
s=s+d^(q+1)*(1-t(q+1));
end;
equal=s;

zuhe2.m文件：
function big=zuhe2(l,j,n,d)
h=log(n)/log(d);
s=0;
for q=0:h;
t(q+1)=1;
for p=l-j+1:n;
t(q+1)=t(q+1)*(1-n/d^q/p);
end
s=s+d^q*t(q+1);
end;
big=s;

zuhe3.m文件：
function small=zuhe3(l,j,d)
small=fix(d*(j-l)/(d-1));
三个文件单独编译都没有问题，但是运行文件drawz.m就有问题了
n=1024;d=2;
[L,J]=meshgrid(0:100:1000)
if L>J
C=zuhe(J,n,d)-zuhe2(L,J,n,d);
else if L<J
C=zuhe(J,n,d)+zuhe3(L,J,d);
else
C=zuhe(L,n,d);
end
end
mesh(L,J,C)
提示的错误在zuhe.m文件

suffer

简单看了一下程序很乱

子程序为标量运算，但传入的变量中又有向量

另外向量的个数也有问题，没有仔细看