code: def is_solenoidal(field):
      """
      Checks if a field is solenoidal.
  
      Parameters
      ==========
  
      field : Vector
          The field to check for solenoidal property
  
      Examples
      ========
  
      >>> from sympy.vector import CoordSys3D
      >>> from sympy.vector import is_solenoidal
      >>> R = CoordSys3D('R')
      >>> is_solenoidal(R.y*R.z*R.i + R.x*R.z*R.j + R.x*R.y*R.k)
      True
      >>> is_solenoidal(R.y * R.j)
      False
  
      """
  
      # Field is solenoidal irrespective of system
      # Take the first coordinate system in the result of the
      # separate method in Vector
      if not isinstance(field, Vector):
          raise TypeError("field should be a Vector")
      if field == Vector.zero:
          return True
      return divergence(field).simplify() is S.Zero
  
imports: from sympy.vector.coordsysrect import CoordSys3D
  from sympy.vector.deloperator import Del
  from sympy.vector.scalar import BaseScalar
  from sympy.vector.vector import Vector, BaseVector
  from sympy.vector.operators import gradient, curl, divergence
  from sympy.core.function import diff
  from sympy.core.singleton import S
  from sympy.integrals.integrals import integrate
  from sympy.simplify.simplify import simplify
  from sympy.core import sympify
  from sympy.vector.dyadic import Dyadic
label_desc: Use image node_id 5 for calling a global function with example usage: is_solenoidal(field) and returns: unknown, True