import sympy


def op_neutral_density_filter(nd):
    """
    Jones matrix operator for a neutral density filter with decadic attenuation.

    Args:
        nd: base ten optical density  [---]
    """
    return sympy.Matrix([[10**-nd / 2, 0], [0, 10**-nd / 2]])