Create cal_jpole_bjcj.m

Calculate J-pole (Pade approximation) coefficient bj and cj for Z_J(zeta)=sum(bj/(zeta-cj))