      t0 = 6.D0/35.D0/R**2*M*sqrt(pi)*sqrt(35.D0)
                                dxC33(1,1)=t0

      t0 = cmplx(0.D0,-6.D0/35.D0)/R**2*M*sqrt(pi)*sqrt(35.D0)
                                dxC33(1,2)=t0

      t0 = 0
                                dxC33(1,3)=t0

      t0 = 0
                                dxC33(1,4)=t0

      t0 = -6.D0/35.D0/R**2*M*sqrt(pi)*sqrt(35.D0)
                                dxC33(2,2)=t0

      t0 = 0
                                dxC33(2,3)=t0

      t0 = 0
                                dxC33(2,4)=t0

      t0 = 0
                                dxC33(3,3)=t0

      t0 = 0
                                dxC33(3,4)=t0

      t0 = 0
                                dxC33(4,4)=t0

