code: def test_compare_vicenty():
      taz, tsr = 38, 3000
      pyproj = pytest.importorskip("pyproj")
  
      lat2, lon2 = vreckon(10, 20, tsr, taz)
  
      p4lon, p4lat, p4a21 = pyproj.Geod(ellps="WGS84").fwd(
          lon2, lat2, taz, tsr
      )
      assert p4lon == approx(lon2, rel=0.0025)
      assert p4lat == approx(lat2, rel=0.0025)
  
      p4az, p4a21, p4sr = pyproj.Geod(ellps="WGS84").inv(
          20, 10, lon2, lat2
      )
      assert (p4az, p4sr) == approx((taz, tsr))
  
imports: import pymap3d
  import pytest
  from pymap3d.vincenty import vreckon
  from pytest import approx
label_desc: Use image node_id 1 for calling a global function with example usage: test_compare_vicenty() without return types