#!/usr/bin/env python

# author: Daniel Scheiermann
# email: daniel.scheiermann@stud.uni-hannover.de
# license: MIT
# Please feel free to use and modify this, but keep the above information.

"""
Script to check a calculation in a paper
and import :math:`a_0` and :math:`u` in kg
"""

from scipy.constants import hbar
from scipy.constants import pi
from scipy.constants import mu_0
from scipy.constants import m_u
import scipy.constants as sci_const

a_0: float = sci_const.physical_constants["Bohr radius"][0]
u_in_kg: float = sci_const.physical_constants["atomic mass constant"][0]

if __name__ == "__main__":
    a_0 = sci_const.physical_constants["Bohr radius"][0]
    mu_b = sci_const.physical_constants["Bohr magneton"][0]
    a = 100 * a_0
    mu = 6 * mu_b
    m_chromium = 52 * m_u
    g = 4 * pi * hbar ** 2 * a / m_chromium

    print(mu_0 * mu ** 2 / (3 * g))
