syms w;f=(abs(0.5*(1+exp(-j*w))./(1-exp(-j*w))-1/(j*w).*exp(-j*w))).^2;val=int(f,w,0,0.75*pi);vpa(val)