code: def test_issue_17247_expression_blowup_32():
      M = Matrix(
          [
              [x + 1, 1 - x, 0, 0],
              [1 - x, x + 1, 0, x + 1],
              [0, 1 - x, x + 1, 0],
              [0, 0, 0, x + 1],
          ]
      )
      with dotprodsimp(True):
          assert M.LUsolve(ones(4, 1)) == Matrix(
              [
                  [(x + 1) / (4 * x)],
                  [(x - 1) / (4 * x)],
                  [(x + 1) / (4 * x)],
                  [1 / (x + 1)],
              ]
          )
  
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 4 for calling a global function with example usage: test_issue_17247_expression_blowup_32() without return types