bpo-42323: Fix math.nextafter() for NaN on AIX (GH-24265)

diff --git a/Modules/mathmodule.c b/Modules/mathmodule.c
index 86b64fb..8133d6b 100644
--- a/Modules/mathmodule.c
+++ b/Modules/mathmodule.c
@@ -3473,6 +3473,12 @@
            Bug fixed in bos.adt.libm 7.2.2.0 by APAR IV95512. */
         return PyFloat_FromDouble(y);
     }
+    if (Py_IS_NAN(x)) {
+        return x;
+    }
+    if (Py_IS_NAN(y)) {
+        return y;
+    }
 #endif
     return PyFloat_FromDouble(nextafter(x, y));
 }