  In this document, we examine exact and efficient numerical approaches to the
MIT Bag Model, a theoretical framework used to describe the properties of bound
quarks in Hadrons. We present the exact and Boundary Value Problem (BVP)
numerical approaches. Both methods are effective in calculating the
eigen-functions and energy levels. Notably, the precision of the BVP approach
matches up to 10 decimal places when compared to the exact approach.
