code: def test_hessian(austourists_model_fit):
      # The hessian approximations are not very consistent, but the test makes
      # sure they run
      austourists_model_fit.model.hessian(austourists_model_fit.params)
      austourists_model_fit.model.hessian(
          austourists_model_fit.params,
          approx_complex_step=False,
          approx_centered=True,
      )
  
imports: from statsmodels.compat.pandas import QUARTER_END
  from statsmodels.compat.platform import PLATFORM_LINUX32, PLATFORM_WIN
  from itertools import product
  import json
  import pathlib
  import numpy
  from numpy.testing import assert_allclose, assert_almost_equal
  import pandas
  import pytest
  import scipy.stats
  from statsmodels.tsa.exponential_smoothing.ets import ETSModel
  import statsmodels.tsa.holtwinters
  import statsmodels.tsa.statespace.exponential_smoothing
label_desc: Use image node_id 31 for calling a global function with example usage: test_hessian(austourists_model_fit) without return types