code: def test_reckon_stability(matmap3d):
      eng = matlab_engine()
  
      if matmap3d:
          if not has_matmap3d(eng):
              pytest.skip("Matmap3d not found")
      else:
          if not has_mapping(eng):
              pytest.skip("Matlab Toolbox not found")
  
      dlast, alast = nan, nan
      lon1, lon2 = 0.0, 1.0
      for i in range(20):
          lat1 = lat2 = 10.0 ** (-i)
  
          dist_m, az_deg = pymap3d.vincenty.vreckon(
              lat1, lon1, lat2, lon2
          )
  
          assert dist_m != dlast
          assert az_deg != alast
          dist_matlab, az_matlab = reckon(
              eng, matmap3d, lat1, lon1, lat2, lon2
          )
  
          assert dist_m == approx(dist_matlab)
  
          assert az_deg == approx(az_matlab, rel=0.005)
  
file_path: pymap3d/src/pymap3d/tests/test_matlab_vreckon.py
imports: from __future__ import annotations
  from math import nan
  import pytest
  from pytest import approx
  import pymap3d.vincenty
label_desc: Use image node_id 2 for calling a global function with example usage: test_reckon_stability(matmap3d) without return types