diff --git a/changelog.md b/changelog.md
--- a/changelog.md
+++ b/changelog.md
@@ -2,6 +2,9 @@
 
 ## WIP
 
+## [0.4.0.1] - 2022-09-06
+  - Raise upper bound on ad
+
 ## 0.4 - 2020-09-29
   - Hyperbolic anomaly support
   - Several new utility functions
diff --git a/default.nix b/default.nix
--- a/default.nix
+++ b/default.nix
@@ -17,8 +17,6 @@
       overrides = self: super:
         {
           exact-real = markUnbroken (dontCheck (doJailbreak super.exact-real));
-          # reanimate = self.callCabal2nix "" ../../src/reanimate {};
-          # reanimate-svg = self.callCabal2nix "" ../../src/reanimate-svg {};
         } // pkgs.lib.optionalAttrs hoogle {
           ghc = super.ghc // { withPackages = super.ghc.withHoogle; };
           ghcWithPackages = self.ghc.withPackages;
diff --git a/orbits.cabal b/orbits.cabal
--- a/orbits.cabal
+++ b/orbits.cabal
@@ -1,13 +1,13 @@
 cabal-version: 1.24
 
--- This file has been generated from package.yaml by hpack version 0.33.1.
+-- This file has been generated from package.yaml by hpack version 0.34.5.
 --
 -- see: https://github.com/sol/hpack
 --
--- hash: b0d20db4cd07be0ae62542dc01fb9fddb137451e34eb3a68439405d4619e86d9
+-- hash: 30b47e80ba31dcb494b9c514befc874eb979af8102ca6646a21a2e39fe61ba3b
 
 name:           orbits
-version:        0.4
+version:        0.4.0.1
 synopsis:       Types and functions for Kepler orbits.
 category:       Physics
 homepage:       https://github.com/expipiplus1/orbits#readme
@@ -45,7 +45,8 @@
       Data.Metrology.Extra
   hs-source-dirs:
       src
-  default-extensions: DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
+  default-extensions:
+      DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
   ghc-options: -Wall -O2
   build-depends:
       ad >=4.3.2
@@ -64,7 +65,8 @@
       
   hs-source-dirs:
       test/doctest
-  default-extensions: DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
+  default-extensions:
+      DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
   ghc-options: -Wall -O2
   build-depends:
       base
@@ -86,7 +88,8 @@
       Paths_orbits
   hs-source-dirs:
       test
-  default-extensions: DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
+  default-extensions:
+      DataKinds FlexibleContexts GeneralizedNewtypeDeriving LambdaCase QuasiQuotes ScopedTypeVariables TemplateHaskell TypeApplications TypeFamilies TypeOperators ViewPatterns
   ghc-options: -Wall -O2 -threaded
   build-depends:
       QuickCheck
diff --git a/src/Physics/Orbit.hs b/src/Physics/Orbit.hs
--- a/src/Physics/Orbit.hs
+++ b/src/Physics/Orbit.hs
@@ -1,3 +1,4 @@
+{-# language CPP #-}
 {-# language QuasiQuotes #-}
 
 -- | Types and functions for dealing with Kepler orbits.
@@ -103,9 +104,16 @@
                                                 , Scalar
                                                 , auto
                                                 )
+#if MIN_VERSION_ad(4,5,0)
+import           Numeric.AD.Rank1.Halley        ( findZero
+                                                , findZeroNoEq
+                                                )
+#else
 import           Numeric.AD.Halley              ( findZero
                                                 , findZeroNoEq
                                                 )
+#endif
+
 import           Numeric.AD.Internal.Identity   ( Id(..) )
 import qualified Numeric.AD.Newton.Double      as Newton
 import           Physics.Orbit.Metrology