From 8299a9073587d12d865d6ce2a07ce7c611a5a22b Mon Sep 17 00:00:00 2001
From: chaitanyabatra <71209117+chaitanyabatra@users.noreply.github.com>
Date: Sun, 1 Oct 2023 12:51:50 +0530
Subject: [PATCH] Create energy_mass_equivalent.py

Fixes #9191

Added energy mass equivalence formulae for stationary as well as moving mass
---
 physics/energy_mass_equivalent.py | 55 +++++++++++++++++++++++++++++++
 1 file changed, 55 insertions(+)
 create mode 100644 physics/energy_mass_equivalent.py

diff --git a/physics/energy_mass_equivalent.py b/physics/energy_mass_equivalent.py
new file mode 100644
index 000000000000..386c95cb4427
--- /dev/null
+++ b/physics/energy_mass_equivalent.py
@@ -0,0 +1,55 @@
+import scipy.constants as const
+import math
+
+
+def energy_equivalent_for_stationary_mass(mass_kg):
+    """
+    Calculate the energy equivalent for a given mass using E=mc^2.
+
+    Parameters:
+        mass_kg (float): Mass in kilograms.(must be positive)
+
+    Returns:
+        float: Energy equivalent in joules.
+
+    >>>energy_equivalent_for_stationary_mass(4)                 
+    3.5950207149472704e+17
+
+    """
+    if mass_kg<0:
+        ValueError("mass of object cannot be negative")
+      
+    speed_of_light = const.speed_of_light  # Speed of light in meters/second
+    energy_joules = mass_kg * (speed_of_light ** 2)
+    return energy_joules
+
+
+def energy_equivalent_for_moving_mass(mass_kg, velocity_m_s):
+    """
+    Calculate the energy equivalent for a moving mass using relativistic energy-momentum relation(E^2= m^2*c^4 +p^2*c^2).
+
+    Parameters:
+        mass_kg (float): Mass in kilograms.(must be positive)
+        velocity_m_s (float): Velocity in meters per second.(can be negative as well as positive)
+
+    Returns:
+        float: Energy equivalent in joules.
+
+    >>>energy_equivalent_for_moving_mass(1,5675)       
+    1701322225868.2546
+    """
+    if mass_kg<0:
+        ValueError("mass of object cannot be negative")
+        
+    speed_of_light = const.speed_of_light  # Speed of light in meters/second
+
+    # Calculating momentum
+    momentum = mass_kg * velocity_m_s / math.sqrt(1 - (velocity_m_s**2 / speed_of_light**2))
+
+    # Calculating energy using the relativistic energy-momentum relation
+    energy_joules = math.sqrt((mass_kg * speed_of_light)**2 + (momentum * speed_of_light)**2)
+    return energy_joules
+
+if __name__ == "__main__":
+    import doctest
+    doctest.testmod()