This study presents computational results for proton radiative capture by triton (3H) using the Argonne V18 (AV18) potential model enhanced with three-body forces. We develop a comprehensive computational framework combining the AV18 nucleon-nucleon potential with Urbana IX three-nucleon forces to calculate the astrophysical S-factor, scattering length, and effective range for the 3H(p, γ)4He reaction. Our methodology involves solving the Schrödinger equation numerically using variational Monte Carlo techniques with explicit treatment of three-body correlations. The results show improved agreement with experimental data compared to two-body potential models, particularly in the low-energy region relevant for astrophysical applications. We provide detailed comparisons with previous theoretical studies and experimental measurements, demonstrating the importance of three-body forces in accurately describing this radiative capture process.