diff --git a/Raphael/DWBA_ZR.py b/Raphael/DWBA_ZR.py
new file mode 100755
index 0000000..92a7417
--- /dev/null
+++ b/Raphael/DWBA_ZR.py
@@ -0,0 +1,15 @@
+#!/usr/bin/env python3
+
+from solveSE import WS, Coulomb, SO, WSSurface, SolvingSE
+
+
+boundState = SolvingSE(16, 8, 1, 0, -4.14)
+boundState.SetRange(0, 0.1, 300)
+boundState.PrintInput()
+
+boundState.ClearPotential()
+boundState.AddPotential(WS(-40, 1.10, 0.65))
+
+boundState.SetLJ(0, 0.5)
+
+print(boundState.GetPotentialValue(0))
\ No newline at end of file