code: def test_work_ind(self):
      def func(par, x):
          b0, b1 = par
          return b0 + b1 * x
  
      # generate some data
      n_data = 4
      x = np.arange(n_data)
      y = np.where(x % 2, x + 0.1, x - 0.1)
      x_err = np.full(n_data, 0.1)
      y_err = np.full(n_data, 0.1)
  
      # do the fitting
      linear_model = Model(func)
      real_data = RealData(x, y, sx=x_err, sy=y_err)
      odr_obj = ODR(real_data, linear_model, beta0=[0.4, 0.4])
      odr_obj.set_job(fit_type=0)
      out = odr_obj.run()
  
      sd_ind = out.work_ind["sd"]
      assert_array_almost_equal(
          out.sd_beta, out.work[sd_ind : sd_ind + len(out.sd_beta)]
      )
  
imports: import tempfile
  import shutil
  import os
  import numpy
  from numpy import pi
  from numpy.testing import assert_array_almost_equal, assert_equal, assert_warns
  import pytest
  from pytest import raises
  from scipy.odr import Data, Model, ODR, RealData, OdrStop, OdrWarning, multilinear, exponential, unilinear, quadratic, polynomial
label_desc: Use image node_id 24 for calling the TestODR obj's underlying member method code with example usage: obj.test_work_ind() and returns: unknown