diff --git a/scripts/plot_cubic.py b/scripts/plot_cubic.py
new file mode 100644
index 0000000..5752d9e
--- /dev/null
+++ b/scripts/plot_cubic.py
@@ -0,0 +1,11 @@
+#!/usr/bin/env python3
+
+import numpy as np
+import matplotlib.pyplot as plt
+
+def cost(x, y):
+    return 2*x**2 + 2*y**2 + 1
+
+fig = plt.figure()
+ax = plt.axes(projection='3d')
+