def comp_r_bias
      r = r_x( -5.1 * Const::MAS2R)
      r = r_y(-17.3 * Const::MAS2R, r)
      r = r_z( 78.0 * Const::MAS2R, r)
      return r
    rescue => e
      raise
    end
def comp_r_bias
      r = r_x( -5.1 * Const::MAS2R)
      r = r_y(-17.3 * Const::MAS2R, r)
      r = r_z( 78.0 * Const::MAS2R, r)
      return r
    rescue => e
      raise
    end
=========================================================================
 Bias 変換行列（一般的な理論）

 * 赤道座標(J2000.0)の極は ICRS の極に対して12時（x軸のマイナス側）の方
   向へ 17.3±0.2 mas、18時（y軸のマイナス側）の方向へ 5.1±0.2 mas ズレ
   ているので、変換する。
   さらに、平均分点への変換はICRSでの赤経を78±10 mas、天の極を中心に回
   転させる。
     18時の方向の変換はx軸を-5.1mas回転
                + 1     0      0   +
       R1(θ ) = | 0   cosθ   sinθ  |
                + 0  -sinθ   cosθ  +
     12時の方向の変換はy軸を-17.3mas回転
                + cosθ   0  -sinθ  +
       R2(θ ) = |   0    1     0   |
                + sinθ   0   cosθ  +
     天の極を中心に78.0mas回転
                +  cosθ   sinθ   0 +
       R3(θ ) = | -sinθ   cosθ   0 |
                +    0      0    1 +

 @param:  <none>
 @return: r  (回転行列)
=========================================================================