Gaussian Process-Homotopy Analysis Method for Tuberculosis Transmission With Stiff Dynamics and Time-Varying Contact Rate. NIJOR 2026, 3 (2), 185-194.