In High Energy Physics, double Higgs production plays a crucial role in assessing the Higgs selfcoupling (trilinear Higgs coupling), responsible for endowing elementary particles with mass and shaping the Higgs potential. Measuring the trilinear Higgs coupling at proton colliders necessitates high luminosity due to the rarity of processes involving it in the Standard Model. Nonetheless, Muon colliders offer distinct advantages over proton colliders, potentially mitigating some measurement challenges associated with the trilinear Higgs coupling. In my research, I have focused on investigating the production of two Higgs particles through the interaction of highenergy muon beams emitting collinear photons. Specifically, we employed the Effective Photon Approximation (EPA), LePDF, and MuonPDF to establish parton distribution functions (PDFs) and determine the total cross sections of these processes. This analysis was conducted within the framework of the Higgs Triplet Model (HTM). In the HTM, a noteworthy hierarchy is observed between the masses of singly charged Higgs boson and doubly charged Higgs boson, and this influence was thoroughly investigated in the context of this research at the 3 TeV Muon Collider. To perform these calculations, we had to manipulate large analytical expressions and generate numerical results. We utilized the BeoShock high-performance computer available at Wichita State University. In this presentation, I will discuss our results for various scenarios (arXiv:2312.12594). The development of advanced future detectors, computer algebraic methods and numerical analysis in particle physics research can have applications beyond the field itself, leading to technological innovations with real-world implications.