diff --git a/CMakeLists.txt b/CMakeLists.txt
index 4b7d4f4a7d20637a144d52637ddf107abbb130fd..73b480fdd48f43154eb6363fdc94717f657fb7b5 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -13,6 +13,7 @@ set(CMAKE_CXX_STANDARD 11)
 # Library setup
 add_library(usgsephem SHARED
             src/TestObject.cpp
+            src/eal.cpp
 )
 set_target_properties(usgsephem PROPERTIES
     VERSION ${PROJECT_VERSION}
diff --git a/include/eal.h b/include/eal.h
new file mode 100644
index 0000000000000000000000000000000000000000..8ea965672de4b51726ef99109d25830294771de3
--- /dev/null
+++ b/include/eal.h
@@ -0,0 +1,15 @@
+#ifndef EAL_H
+#define EAL_H
+
+#include <vector>
+#include <string.h>
+
+using namespace std;
+
+vector<double> getPosition(vector<vector<double>> coords, vector<double> coeffs, char *interp, double time);
+vector<double> getVelocity(vector<vector<double>> coords, vector<double> coeffs, char *interp, double time, bool interpolation);
+
+vector<double> getPosition(vector<double> coeffs, char *interp, double time);
+vector<double> getVelocity(vector<double> coeffs, char *interp, double time, bool interpolation);
+
+#endif // EAL_H
diff --git a/src/eal.cpp b/src/eal.cpp
new file mode 100644
index 0000000000000000000000000000000000000000..b523c36ce7f82c289f2a58ace1262f9a856c8218
--- /dev/null
+++ b/src/eal.cpp
@@ -0,0 +1,49 @@
+#include "eal.h"
+
+// Positional Functions
+
+// Position Data Functions
+vector<double> getPosition(vector<vector<double>> coords, vector<double> coeffs, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+vector<double> getVelocity(vector<vector<double>> coords, vector<double> coeffs, char *interp, double time, bool interpolation) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+// Postion Function Functions
+vector<double> getPosition(vector<double> coeffs, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+vector<double> getVelocity(vector<double> coeffs, char *interp, double time, bool interpolation) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+// Rotation Data Functions
+
+// Rotation Data Functions
+vector<double> getRotation(char *from, char *to, vector<vector<double>> rotations, vector<double> times, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+vector<double> getAngularVelocity(char *from, char *to, vector<vector<double>> rotations, vector<double> times, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+// Rotation Function Functions
+vector<double> getRotation(char *from, char *to, vector<double> coefficients, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}
+
+vector<double> getAngularVelocity(char *from, char *to, vector<double> coefficients, char *interp, double time) {
+  vector<double> coordinate = {0.0, 0.0, 0.0};
+  return coordinate;
+}