      t0 = -4.D0/15.D0*sqrt(pi)*M*sqrt(5.D0)/R
                                 Y20(1,1)=t0

      t0 = 0
                                 Y20(1,2)=t0

      t0 = 0
                                 Y20(1,3)=t0

      t0 = 0
                                 Y20(1,4)=t0

      t0 = -4.D0/15.D0*sqrt(pi)*M*sqrt(5.D0)/R
                                 Y20(2,2)=t0

      t0 = 0
                                 Y20(2,3)=t0

      t0 = 0
                                 Y20(2,4)=t0

      t0 = 8.D0/15.D0*sqrt(pi)*M*sqrt(5.D0)/R
                                 Y20(3,3)=t0

      t0 = 0
                                 Y20(3,4)=t0

      t0 = 0
                                 Y20(4,4)=t0

