def METHOD_NAME(disk_from_geometry):
    # fmt: off
    Md1 = np.array([[ 32.58973,   0.     ,   0.     ,   0.     ],
                    [  0.     ,  32.58973,   0.     ,   0.     ],
                    [  0.     ,   0.     ,   0.17809,   0.     ],
                    [  0.     ,   0.     ,   0.     ,   0.17809]])
    # fmt: on
    assert_almost_equal(disk_from_geometry.M(), Md1, decimal=5)