code: def test_arithmetic():
      a = MutableDenseNDimArray([3 for i in range(9)], (3, 3))
      b = MutableDenseNDimArray([7 for i in range(9)], (3, 3))
  
      c1 = a + b
      c2 = b + a
      assert c1 == c2
  
      d1 = a - b
      d2 = b - a
      assert d1 == d2 * (-1)
  
      e1 = a * 5
      e2 = 5 * a
      e3 = copy(a)
      e3 *= 5
      assert e1 == e2 == e3
  
      f1 = a / 5
      f2 = copy(a)
      f2 /= 5
      assert f1 == f2
      assert (
          f1[0, 0]
          == f1[0, 1]
          == f1[0, 2]
          == f1[1, 0]
          == f1[1, 1]
          == f1[1, 2]
          == f1[2, 0]
          == f1[2, 1]
          == f1[2, 2]
          == Rational(3, 5)
      )
  
      assert (
          type(a)
          == type(b)
          == type(c1)
          == type(c2)
          == type(d1)
          == type(d2)
          == type(e1)
          == type(e2)
          == type(e3)
          == type(f1)
      )
  
      z0 = -a
      assert z0 == MutableDenseNDimArray([-3 for i in range(9)], (3, 3))
  
imports: from copy import copy
  from sympy.tensor.array.dense_ndim_array import MutableDenseNDimArray
  from sympy.core.function import diff
  from sympy.core.numbers import Rational
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.core.sympify import sympify
  from sympy.matrices import SparseMatrix
  from sympy.matrices import Matrix
  from sympy.tensor.array.sparse_ndim_array import MutableSparseNDimArray
  from sympy.testing.pytest import raises
label_desc: Use image node_id 11 for calling a global function with example usage: test_arithmetic() without return types