diff --git a/examples/cavity/binomial_code.jl b/examples/cavity/binomial_code.jl
new file mode 100644
index 00000000..177ed09b
--- /dev/null
+++ b/examples/cavity/binomial_code.jl
@@ -0,0 +1,13 @@
+using QuantumCollocation
+using NamedTrajectories
+
+transmon_levels = 3
+cavity_levels = 5
+
+system = MultiModeSystem(transmon_levels, cavity_levels)
+
+g0 = multimode_state("g0", transmon_levels, cavity_levels)
+g1 = multimode_state("g1", transmon_levels, cavity_levels)
+g4 = multimode_state("g4", transmon_levels, cavity_levels)
+
+ψ_init = [g0, g1]