To complete the topic of this thread professionally, I will answer the question that no one asked me:
What is the well-known work on which your 3 formulas were based and derived?
Its answer is simple to those who are familiar with LTspice simulator.
They are derived from the work of John Chan in modeling a non-linear coil.
To those who don’t have LTspice, I attached a PDF of the LTspice help page which is related to John Chan model.
His two formulas, one for the positive BH branch and the other for the negative one, are the Bup(H) and Bdn(H) respectively.
Obviously, if someone will be able to prove that any of my 3 formulas is not valid, he simply proves me that John Chan model is also not valid.
Regards,
Kerim