function z=rem_theta(n,b)
     x=sym(n);
     z=theta_r(x,b)+x/2-x/2*log(x/2/sym('pi'))+sym('pi')/8;