function z=theta_r(n,b)
     x=sym(n);
     z=imag(log_fact_complex(-3/4,x/2,b))-x/2*log(sym('pi'));