code: def test_LUsolve_iszerofunc():
      # taken from https://github.com/sympy/sympy/issues/24679
  
      M = Matrix([[(x + 1) ** 2 - (x**2 + 2 * x + 1), x], [x, 0]])
      b = Matrix([1, 1])
      is_zero_func = lambda e: False if e._random() else True
  
      x_exp = Matrix(
          [1 / x, (1 - (-(x**2) - 2 * x + (x + 1) ** 2 - 1) / x) / x]
      )
  
      assert (x_exp - M.LUsolve(b, iszerofunc=is_zero_func)) == Matrix(
          [0, 0]
      )
  
imports: import pytest
  from sympy.core.function import expand_mul
  from sympy.core.numbers import I, Rational
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol, symbols
  from sympy.core.sympify import sympify
  from sympy.simplify.simplify import simplify
  from sympy.matrices.exceptions import ShapeError, NonSquareMatrixError
  from sympy.matrices import ImmutableMatrix, Matrix, eye, ones, ImmutableDenseMatrix, dotprodsimp
  from sympy.matrices.determinant import _det_laplace
  from sympy.testing.pytest import raises
  from sympy.matrices.exceptions import NonInvertibleMatrixError
  from sympy.polys.matrices.exceptions import DMShapeError
  from sympy.solvers.solveset import linsolve
  from sympy.abc import x, y
label_desc: Use image node_id 3 for calling a global function with example usage: test_LUsolve_iszerofunc() without return types