From 604009db9a20290810634238e1a4d7942afb09ff Mon Sep 17 00:00:00 2001
From: pulipakaa24 <pulipakaaditya@gmail.com>
Date: Wed, 19 Nov 2025 22:00:13 -0600
Subject: [PATCH] find equilibrium height

---
 equilibrium.py | 13 +++++++++++++
 1 file changed, 13 insertions(+)
 create mode 100644 equilibrium.py

diff --git a/equilibrium.py b/equilibrium.py
new file mode 100644
index 0000000..759837c
--- /dev/null
+++ b/equilibrium.py
@@ -0,0 +1,13 @@
+import sympy as sp
+
+g = 9.81
+m = 6
+z = sp.symbols('z')
+term1 = (6.167266375e8*((0.4446789110e5)**2))/((0.2466906550e10*z + 0.6886569338e8)**2)
+term2 = 0.3042813963e19*z*((0.4446789110e5)**2)/((0.2466906550e10*z + 0.6886569338e8)**3)
+equation = sp.Eq(4*(term1 - term2), m*g)
+solution = sp.solve(equation, z)
+print("Equilibrium position (z):")
+for sol in solution:
+    if sol.is_real and sol > 0:
+        print(sol.evalf())
\ No newline at end of file