code: def test_legmul(self):
      # check values of result
      for i in range(5):
          pol1 = [0] * i + [1]
          val1 = leg.legval(self.x, pol1)
          for j in range(5):
              msg = f"At i={i}, j={j}"
              pol2 = [0] * j + [1]
              val2 = leg.legval(self.x, pol2)
              pol3 = leg.legmul(pol1, pol2)
              val3 = leg.legval(self.x, pol3)
              assert_(len(pol3) == i + j + 1, msg)
              assert_almost_equal(val3, val1 * val2, err_msg=msg)
  
imports: from functools import reduce
  import numpy
  import numpy.polynomial.legendre
  from numpy.polynomial.polynomial import polyval
  from numpy.testing import assert_almost_equal, assert_raises, assert_equal, assert_
label_desc: Use image node_id 4 for calling the TestArithmetic obj's underlying member method code with example usage: obj.test_legmul() without return types