diff --git a/src/volume.rs b/src/volume.rs
index 1f69f6c..9286f5c 100644
--- a/src/volume.rs
+++ b/src/volume.rs
@@ -1664,45 +1664,45 @@ pub fn kappa(sa: f64, ct: f64, p: f64) -> Result<f64> {
     let pi: f64 = non_dimensional_p(p);
 
     let v = V000
-        + s * (V010 + s * (V020 + s * (V030 + s * (V040 + s * (V050 + V060 * s)))))
+        + s * (V100 + s * (V200 + s * (V300 + s * (V400 + s * (V500 + V600 * s)))))
         + tau
-            * (V100
-                + s * (V110 + s * (V120 + s * (V130 + s * (V140 + V150 * s))))
+            * (V010
+                + s * (V110 + s * (V210 + s * (V310 + s * (V410 + V510 * s))))
                 + tau
-                    * (V200
-                        + s * (V210 + s * (V220 + s * (V230 + V240 * s)))
+                    * (V020
+                        + s * (V120 + s * (V220 + s * (V320 + V420 * s)))
                         + tau
-                            * (V300
-                                + s * (V310 + s * (V320 + V330 * s))
+                            * (V030
+                                + s * (V130 + s * (V230 + V330 * s))
                                 + tau
-                                    * (V400
-                                        + s * (V410 + V420 * s)
-                                        + tau * (V500 + V510 * s + V600 * tau)))))
+                                    * (V040
+                                        + s * (V140 + V240 * s)
+                                        + tau * (V050 + V150 * s + V060 * tau)))))
         + pi * (V001
-            + s * (V011 + s * (V021 + s * (V031 + s * (V041 + V051 * s))))
+            + s * (V101 + s * (V201 + s * (V301 + s * (V401 + V501 * s))))
             + tau
-                * (V101
-                    + s * (V111 + s * (V121 + s * (V131 + V141 * s)))
+                * (V011
+                    + s * (V111 + s * (V211 + s * (V311 + V411 * s)))
                     + tau
-                        * (V201
-                            + s * (V211 + s * (V221 + V231 * s))
+                        * (V021
+                            + s * (V121 + s * (V221 + V321 * s))
                             + tau
-                                * (V301
-                                    + s * (V311 + V321 * s)
-                                    + tau * (V401 + V411 * s + V501 * tau))))
+                                * (V031
+                                    + s * (V131 + V231 * s)
+                                    + tau * (V041 + V141 * s + V051 * tau))))
             + pi * (V002
-                + s * (V012 + s * (V022 + s * (V032 + V042 * s)))
+                + s * (V102 + s * (V202 + s * (V302 + V402 * s)))
                 + tau
-                    * (V102
-                        + s * (V112 + s * (V122 + V132 * s))
+                    * (V012
+                        + s * (V112 + s * (V212 + V312 * s))
                         + tau
-                            * (V202
-                                + s * (V212 + V222 * s)
-                                + tau * (V302 + V312 * s + V402 * tau)))
+                            * (V022
+                                + s * (V122 + V222 * s)
+                                + tau * (V032 + V132 * s + V042 * tau)))
                 + pi * (V003
-                    + s * (V013 + V023 * s)
-                    + tau * (V103 + V113 * s + V203 * tau)
-                    + pi * (V004 + V014 * s + V104 * tau + pi * (V005 + V006 * pi)))));
+                    + s * (V103 + V203 * s)
+                    + tau * (V013 + V113 * s + V023 * tau)
+                    + pi * (V004 + V104 * s + V014 * tau + pi * (V005 + V006 * pi)))));
 
     let v_p = C000
         + s * (C100 + s * (C200 + s * (C300 + s * (C400 + C500 * s))))