diff --git a/include/units/Vector2D.hpp b/include/units/Vector2D.hpp
index 865d9fe..a144176 100644
--- a/include/units/Vector2D.hpp
+++ b/include/units/Vector2D.hpp
@@ -193,19 +193,6 @@ template <isQuantity T> class Vector2D {
          */
         constexpr T magnitude() const { return sqrt(square(this->x) + square(this->y)); }
 
-        /**
-         * @brief difference between two vectors
-         *
-         * This function calculates the difference between two vectors
-         * a.vectorTo(b) = {b.x - a.x, b.y - a.y}
-         *
-         * @param other the other vector
-         * @return Vector2D<T>
-         */
-        constexpr Vector2D<T> vectorTo(const Vector2D<T>& other) const {
-            return Vector2D<T>(other.x - this->x, other.y - this->y);
-        }
-
         /**
          * @brief the angle between two vectors
          *
diff --git a/include/units/Vector3D.hpp b/include/units/Vector3D.hpp
index 6cf211e..1caacee 100644
--- a/include/units/Vector3D.hpp
+++ b/include/units/Vector3D.hpp
@@ -220,21 +220,6 @@ template <isQuantity T> class Vector3D {
          */
         constexpr T magnitude() const { return sqrt(square(this->x) + square(this->y) + square(this->z)); }
 
-        /**
-         * @brief difference between two vectors
-         *
-         * TODO: figure out if this is even necessary, you could just use the - operator overload
-         *
-         * This function calculates the difference between two vectors
-         * a.vectorTo(b) = {b.x - a.x, b.y - a.y, b.z - a.z}
-         *
-         * @param other the other vector
-         * @return Vector3D<T>
-         */
-        constexpr Vector3D<T> vectorTo(const Vector3D<T>& other) const {
-            return Vector2D<T>(other.x - this->x, other.y - this->y, other.z - this->z);
-        }
-
         /**
          * @brief the angle between two vectors
          *