现有以下程序：
maxab = max(alpha,beta);
absab = abs(alpha-beta);
param2 = s+1;
if alpha < beta
    param2 = r+X;
end   

part1 = (alpha^r*beta^s/gamma(r))*gamma(r+X);
part2 = 1./((alpha+T).^(r+X).*(beta+T).^s);
if absab == 0
   F1=1./((maxab+TX).^(r+s+X));
   F2=1./((maxab+T).^(r+s+X));
   su=(F1-F2);
else
   F1=h2f1(r+s+X,param2,r+s+X+1,absab./(maxab+TX))./...
       ((maxab+TX).^(r+s+X));
   F2=h2f1(r+s+X,param2,r+s+X+1,absab./(maxab+T))./...
       ((maxab+T).^(r+s+X));
end

f = -sum(log(part1.*(part2+(s./(r+s+X)).*(F1-F2))));
[f/1000 param]
g=[];
运行过程中，程序报错为：Warning: Log of zero.也就是说part1.*(part2+(s./(r+s+X)).*(F1-F2))部分出现了零值，本人想写个条件语句避开零值的情况，将后面部分的程序改为：
if (part1.*(part2+(s./(r+s+X)).*(F1-F2)))>0
    f=-sum(log(part1.*(part2+(s./(r+s+X)).*(F1-F2))));
    [f/1000 param]
    g=[];
else if (part1.*(part2+(s./(r+s+X)).*(F1-F2)))<0
        f=-sum(log(part1.*(part2+(s./(r+s+X)).*(F1-F2))));
        [f/1000 param]
        g=[];
    end
end
但还是报错为：One or more output arguments not assigned during call to 'D:\MATLAB\work\pareto_nbd_ll.m (pareto_nbd_ll)'.
我是菜鸟，一直搞不懂为什么不对。请各位指点啊，论坛币奉上。