code: def test_legadd(self):
      for i in range(5):
          for j in range(5):
              msg = f"At i={i}, j={j}"
              tgt = np.zeros(max(i, j) + 1)
              tgt[i] += 1
              tgt[j] += 1
              res = leg.legadd([0] * i + [1], [0] * j + [1])
              assert_equal(trim(res), trim(tgt), 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 1 for calling the TestArithmetic obj's underlying member method code with example usage: obj.test_legadd() without return types