Create single_iter_lie.py

CFS/single_iter_lie.py

@@ -0,0 +1,8 @@
+def single_iter_lie(eta, h,vector_field,z):
+    # The list that will contain all the Lie derivatives is initiated. 
+    LT = sp.Matrix([h])
+
+    for i in eta:
+
+        LT = LT.jacobian(z)*vector_field[:,i]
+    return LT