function z=log_fact_complex(m,n,b)
     a=sym(m);
     bi=sym(n);
     x=a+bi*i;
%     b
%     ff=bern_sum(x,b)
     z=(x+1/2)*logcomplex(m,n)-x+log(sqrt(2*sym('pi')))+bern_sum(x,b);
%+1/(12*x)-1/(360*x^3)+1/(1260*x^5)-1/(1680*x^7);