code: def test_singular(self):
      # 3) Singular points in region of integration.
      def myfunc(x):
          if 0 < x < 2.5:
              return sin(x)
          elif 2.5 <= x <= 5.0:
              return exp(-x)
          else:
              return 0.0
  
      assert_quad(
          quad(myfunc, 0, 10, points=[2.5, 5.0]),
          1 - cos(2.5) + exp(-2.5) - exp(-5.0),
      )
  
imports: import sys
  import math
  import numpy
  from numpy import sqrt, cos, sin, arctan, exp, log, pi
  from numpy.testing import assert_, assert_allclose, assert_array_less, assert_almost_equal
  import pytest
  from scipy.integrate import quad, dblquad, tplquad, nquad
  from scipy.special import erf, erfc
  from scipy._lib._ccallback import LowLevelCallable
  import ctypes
  import ctypes.util
  from scipy._lib._ccallback_c import sine_ctypes
  import scipy.integrate._test_multivariate
label_desc: Use image node_id 3 for calling the TestQuad obj's underlying member method code with example usage: obj.test_singular() and returns: sin, exp, int