code: def numerical_convolve(func1, func2, interval, dt):
      mni, mxi = interval
      time = np.arange(mni, mxi, dt)
      vec1 = func1(time).astype(float)
      vec2 = func2(time).astype(float)
      value = np.convolve(vec1, vec2) * dt
      min_s = min(time.size, value.size)
      time = time[:min_s]
      value = value[:min_s]
      return time, value
  
imports: import re
  import numpy
  import pytest
  import sympy
  from numpy.testing import assert_almost_equal, assert_array_almost_equal, assert_array_equal
  from sympy import DiracDelta, Dummy, Function, Symbol
  from sympy.utilities.lambdify import implemented_function, lambdify
  from nipy.algorithms.statistics.formula import Term
  from ..None import hrf
  from ..utils import Interp1dNumeric, TimeConvolver, blocks, convolve_functions, define, events, interp, lambdify_t, linear_interp, step_function
label_desc: Use image node_id 7 for calling a global function with example usage: numerical_convolve(func1, func2, interval, dt) and returns: time,  value