diff --git a/system.jl b/system.jl
index c557325..d52f6c2 100644
--- a/system.jl
+++ b/system.jl
@@ -52,12 +52,3 @@ Z_or_N(s::system, p::Bool)::Int = p ? s.Z : s.N
 
 "Create an empty array for the size of the mesh"
 zero_array(s::system) = zeros(1 + s.divs)
-
-"Create linear interpolations for the meson fields"
-function fields_interp(s::system)
-    S_interp = linear_interpolation(rs(s), s.Φ0)
-    V_interp = linear_interpolation(rs(s), s.W0)
-    R_interp = linear_interpolation(rs(s), s.B0)
-    A_interp = linear_interpolation(rs(s), s.A0)
-    return (S_interp, V_interp, R_interp, A_interp)
-end