diff --git a/src/shape_functions.jl b/src/shape_functions.jl
new file mode 100644
index 0000000..f035904
--- /dev/null
+++ b/src/shape_functions.jl
@@ -0,0 +1,19 @@
+function C2D4(ξ::Float64, η::Float64)
+     values::Array{Float64, 1} = [(1 - ξ) * (1 - η),
+                                  ξ * (1 - η),
+                                  ξ * η,
+                                  (1 - ξ) * η]
+    return values
+end
+
+function C3D8(ξ::Float64, η::Float64, μ::Float64)
+    values::Array{Float64, 1} = 1 / 8 * [(1 - ξ) * (1 - η) * (1 - μ),
+                                         (1 + ξ) * (1 - η) * (1 - μ),
+                                         (1 + ξ) * (1 + η) * (1 - μ),
+                                         (1 - ξ) * (1 + η) * (1 - μ),
+                                         (1 - ξ) * (1 - η) * (1 + μ),
+                                         (1 + ξ) * (1 - η) * (1 + μ),
+                                         (1 + ξ) * (1 + η) * (1 + μ),
+                                         (1 - ξ) * (1 + η) * (1 + μ)]
+    return values
+end