def comp_r_bias_prec_nut
      gamma = comp_gamma_bp
      phi   = comp_phi_bp
      psi   = comp_psi_bp
      fj2 = -2.7774e-6 * @jc
      dpsi_ls, deps_ls = compute_lunisolar
      dpsi_pl, deps_pl = compute_planetary
      dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl
      dpsi += dpsi * (0.4697e-6 + fj2)
      deps += deps * fj2
      r = r_z(gamma)
      r = r_x(phi,          r)
      r = r_z(-psi - dpsi,  r)
      r = r_x(-@eps - deps, r)
      return r
    rescue => e
      raise
    end
def comp_r_bias_prec_nut
      gamma = comp_gamma_bp
      phi   = comp_phi_bp
      psi   = comp_psi_bp
      fj2 = -2.7774e-6 * @jc
      dpsi_ls, deps_ls = compute_lunisolar
      dpsi_pl, deps_pl = compute_planetary
      dpsi, deps = dpsi_ls + dpsi_pl, deps_ls + deps_pl
      dpsi += dpsi * (0.4697e-6 + fj2)
      deps += deps * fj2
      r = r_z(gamma)
      r = r_x(phi,          r)
      r = r_z(-psi - dpsi,  r)
      r = r_x(-@eps - deps, r)
      return r
    rescue => e
      raise
    end
=========================================================================
 Bias + Precession + Nutation 変換行列

 * IAU 2006 (Fukushima-Williams 4-angle formulation) 理論

 @param:  <none>
 @return: r  (変換行列)
=========================================================================