This study explores a novel mathematical model for analyzing methimazole (MMI) therapy in patients with Graves' disease, a common autoimmune cause of hyperthyroidism. The model leverages a system of ordinary differential equations to capture the interplay between key factors: methimazole concentration, thyroid-stimulating hormone (TSH), free thyroxine (FT4) concentration, thyroid-stimulating hormone receptor antibody (TRAb) concentration, and functional thyroid gland size (V). Thirteen parameters account for the underlying physiological processes. Stability analysis were employedby Lyapunov stability theory and LaSalle's principle, to demonstrate the model's asymptotic stability. This implies the system converges towards a steady state under defined conditions. Additionally, a clinically relevant chart is constructed based on FT4 levels and treatment duration.For model validation, patient data is utilized within SimBiology simulations implemented in MATLAB. These simulations offer valuable insights into the dynamic response of hyperthyroidism to MMI treatment. This approach holds promise for future applications in optimizing treatment strategies and facilitating personalized medicine approaches.