code: def test_sens_ic_vector_3_params(self):
      # Big System with Many Parameters
      def ode_func_5(y, t, p):
          dx = p[0] * (y[1] - y[0])
          ds = y[0] * (p[1] - y[2]) - y[1]
          dz = y[0] * y[1] - p[2] * y[2]
  
          return [dx, ds, dz]
  
      # Instantiate ODE model
      model5 = DifferentialEquation(
          func=ode_func_5, t0=0, times=self.t, n_states=3, n_theta=3
      )
  
      # First three columns are derivatives with respect to ode parameters
      # Last three coluimns are derivatives with repsect to initial condition
      # So identity matrix should appear in last 3 columns
      model5_sens_ic = np.array(
          [[1, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0], [0, 0, 1, 0, 0, 0]]
      )
  
      np.testing.assert_array_equal(
          np.ravel(model5_sens_ic), model5._sens_ic
      )
  
imports: import warnings
  import numpy
  import pytensor
  import pytensor.tensor
  import pytest
  from scipy.stats import norm
  import pymc
  from pymc.ode import DifferentialEquation
  from pymc.testing import fast_unstable_sampling_mode
label_desc: Use image node_id 6 for calling the TestSensitivityInitialCondition obj's underlying member method code with example usage: obj.test_sens_ic_vector_3_params() and returns: list