diff --git a/src/arap_precompute.cpp b/src/arap_precompute.cpp
new file mode 100644
index 0000000..09f8650
--- /dev/null
+++ b/src/arap_precompute.cpp
@@ -0,0 +1,14 @@
+#include "../include/arap_precompute.h"
+#include <igl/min_quad_with_fixed.h>
+#include <igl/arap_linear_block.h>
+#include <igl/cotmatrix.h>
+
+void arap_precompute(
+  const Eigen::MatrixXd & V,
+  const Eigen::MatrixXi & F,
+  const Eigen::VectorXi & b,
+  igl::min_quad_with_fixed_data<double> & data,
+  Eigen::SparseMatrix<double> & K)
+{
+// REPLACE WITH YOUR CODE
+}
diff --git a/src/arap_single_iteration.cpp b/src/arap_single_iteration.cpp
new file mode 100644
index 0000000..01d98c6
--- /dev/null
+++ b/src/arap_single_iteration.cpp
@@ -0,0 +1,12 @@
+#include "../include/arap_single_iteration.h"
+#include <igl/polar_svd3x3.h>
+#include <igl/min_quad_with_fixed.h>
+
+void arap_single_iteration(
+  const igl::min_quad_with_fixed_data<double> & data,
+  const Eigen::SparseMatrix<double> & K,
+  const Eigen::MatrixXd & bc,
+  Eigen::MatrixXd & U)
+{
+  // REPLACE WITH YOUR CODE
+}
diff --git a/src/biharmonic_precompute.cpp b/src/biharmonic_precompute.cpp
new file mode 100644
index 0000000..205bf80
--- /dev/null
+++ b/src/biharmonic_precompute.cpp
@@ -0,0 +1,15 @@
+#include "biharmonic_precompute.h"
+#include <igl/min_quad_with_fixed.h>
+#include <igl/cotmatrix.h>
+#include <igl/massmatrix.h>
+#include <igl/invert_diag.h>
+
+void biharmonic_precompute(
+  const Eigen::MatrixXd & V,
+  const Eigen::MatrixXi & F,
+  const Eigen::VectorXi & b,
+  igl::min_quad_with_fixed_data<double> & data)
+{
+  // REPLACE WITH YOUR CODE
+  data.n = V.rows();
+}
diff --git a/src/biharmonic_solve.cpp b/src/biharmonic_solve.cpp
new file mode 100644
index 0000000..e5ce727
--- /dev/null
+++ b/src/biharmonic_solve.cpp
@@ -0,0 +1,11 @@
+#include "biharmonic_solve.h"
+#include <igl/min_quad_with_fixed.h>
+
+void biharmonic_solve(
+  const igl::min_quad_with_fixed_data<double> & data,
+  const Eigen::MatrixXd & bc,
+  Eigen::MatrixXd & D)
+{
+  // REPLACE WITH YOUR CODE
+  D = Eigen::MatrixXd::Zero(data.n,3);
+}