diff --git a/py3nj/wigner.py b/py3nj/wigner.py
index f9d15c1..01fd012 100644
--- a/py3nj/wigner.py
+++ b/py3nj/wigner.py
@@ -162,11 +162,11 @@ def wigner9j(two_l1, two_l2, two_l3, two_l4, two_l5, two_l6,
         two_l1, two_l2, two_l3, two_l4, two_l5, two_l6, two_l7, two_l8, two_l9)
 
     # calculat summation directory from 6-j symbol
-    x, sixj = drc6j(np.stack([two_lf, two_lh, two_la]),
-                    np.stack([two_lb, two_lb, two_lj]),
+    x, sixj = drc6j(np.stack([two_lf, two_lb, two_la]),
+                    np.stack([two_lb, two_lf, two_lj]),
                     np.stack([two_lc, two_le, two_lg]),
-                    np.stack([two_la, two_lf, two_lh]),
-                    np.stack([two_lj, two_ld, two_ld]))
+                    np.stack([two_la, two_ld, two_lh]),
+                    np.stack([two_lj, two_lh, two_ld]))
     neg_phase = (2 * x) % 4 - 1
     return np.sum(-neg_phase * (x + 1) * sixj[0] * sixj[1] * sixj[2], axis=-1)
 

wigner 9j disagrees with online calculators
I want to know what following 9j symbol evaluates to
{ 3/2 3/2 1
  2    2    2
1/2 1/2 1 } 

A couple of online calculators say the value is around 0.048:

- http://www-stone.ch.cam.ac.uk/cgi-bin/wigner.cgi?symbol=9j&a=3%2F2&b=3%2F2&c=1&d=2&e=2&f=2&g=1%2F2&h=1%2F2&i=1

- http://dentonwoods.nfshost.com/3j6j9j.html

- py3nj says the value is zero

- sympy.physics.wigner throws an error Value error saying the values must be integer or half integer and fulfill the triangle relation.

I don't know who's right?
diff --git a/tests/test_9j.py b/tests/test_9j.py
index 24df4bf..78ba8bc 100644
--- a/tests/test_9j.py
+++ b/tests/test_9j.py
@@ -32,3 +32,10 @@ def test_wigner9j():
     expected = np.array([value for _, value in NINE_J]).T
     actual = wigner9j(*nine_j)
     assert np.allclose(actual, expected)
+
+
+def test_issue2():
+  # https://github.com/fujiisoup/py3nj/issues/2
+  actual = wigner9j(0, 0, 0, 2 * 1, 0, 2 * 1, 2 * 1, 0, 2 * 1)
+  expected = 1.0 / 3.0
+  assert np.allclose(actual, expected)
\ No newline at end of file

pip install -e .
parse_log_pytest
numpy>=1.16.0
pytest
apt-get update apt-get install -y gfortran
3.9
pytest --no-header -rA --tb=line --color=no -p no:cacheprovider -W ignore::DeprecationWarning
exceptiongroup==1.2.2
iniconfig==2.1.0
numpy @ file:///croot/numpy_and_numpy_base_1736283260865/work/dist/numpy-2.0.2-cp39-cp39-linux_x86_64.whl#sha256=3387e3e62932fa288bc18e8f445ce19e998b418a65ed2064dd40a054f976a6c7
packaging==24.2
pluggy==1.5.0
-e git+https://github.com/fujiisoup/py3nj.git@df2d4391a9869e61a6c48ea1b1bab0625bc5291c#egg=py3nj
pytest==8.3.5
tomli==2.2.1

name: py3nj
channels:
  - defaults
  - https://repo.anaconda.com/pkgs/main
  - https://repo.anaconda.com/pkgs/r
  - conda-forge
dependencies:
  - _libgcc_mutex=0.1=main
  - _openmp_mutex=5.1=1_gnu
  - blas=1.0=openblas
  - ca-certificates=2025.2.25=h06a4308_0
  - ld_impl_linux-64=2.40=h12ee557_0
  - libffi=3.4.4=h6a678d5_1
  - libgcc-ng=11.2.0=h1234567_1
  - libgfortran-ng=11.2.0=h00389a5_1
  - libgfortran5=11.2.0=h1234567_1
  - libgomp=11.2.0=h1234567_1
  - libopenblas=0.3.21=h043d6bf_0
  - libstdcxx-ng=11.2.0=h1234567_1
  - ncurses=6.4=h6a678d5_0
  - numpy=2.0.2=py39heeff2f4_0
  - numpy-base=2.0.2=py39h8a23956_0
  - openssl=3.0.16=h5eee18b_0
  - pip=25.0=py39h06a4308_0
  - python=3.9.21=he870216_1
  - readline=8.2=h5eee18b_0
  - setuptools=72.1.0=py39h06a4308_0
  - sqlite=3.45.3=h5eee18b_0
  - tk=8.6.14=h39e8969_0
  - tzdata=2025a=h04d1e81_0
  - wheel=0.45.1=py39h06a4308_0
  - xz=5.6.4=h5eee18b_1
  - zlib=1.2.13=h5eee18b_1
  - pip:
      - exceptiongroup==1.2.2
      - iniconfig==2.1.0
      - packaging==24.2
      - pluggy==1.5.0
      - pytest==8.3.5
      - tomli==2.2.1
prefix: /opt/conda/envs/py3nj

tests/test_9j.py::test_issue2
tests/test_9j.py::test_wigner9j
py3nj/wigner.py
py3nj/wigner.py