{"references": ["Nag PK. Engineering thermodynamics. Tata McGraw Hill; 2013.", "Incropera FP, DeWitt DP, Bergman TL, Lavine AS, Incropera, F.P., Dewitt D.P., Bergman, T., Lavine A. Fundamentals of Heat and Mass Transfer. 8th ed. Hoboken, NJ: Wiley's; 2007. doi:10.1016/j.applthermaleng.2011.03.022.", "Huei YC. Benefits and introduction to python programming for freshmore students using inexpensive robots. Proc. IEEE Int. Conf. Teaching, Assess. Learn. Eng. Learn. Futur. Now, TALE 2014, 2015, p. 12\u20137. doi:10.1109/TALE.2014.7062611.", "Lin JWB. Why python is the next wave in earth sciences computing. Bull Am Meteorol Soc 2012;93:1823\u20134. doi:10.1175/BAMS-D-12-00148.1.", "Dumka P, Pawar PS, Sauda A, Shukla G, Mishra DR. Application of He's homotopy and perturbation method to solve heat transfer equations: A python approach. Adv Eng Softw 2022;170:103160. doi:10.1016/j.advengsoft.2022.103160.", "Dumka P, Deo A, Gajula K, Sharma V, Chauhan R, Mishra DR. Load and Load Duration Curves Using Python. Int J All Res Educ Sci Methods 2022;10:2127\u201334.", "Meurer A, Smith CP, Paprocki M, \u010cert\u00edk O, Kirpichev SB, Rocklin M, et al. SymPy: Symbolic computing in python. PeerJ Comput Sci 2017;2017:1\u201327. doi:10.7717/peerj-cs.103.", "Rocklin M, Terrel AR. Symbolic statistics with SymPy. Comput Sci Eng 2012;14:88\u201393. doi:10.1109/MCSE.2012.56.", "Rocklin M. Uncertainty Modeling with SymPy Stats. Proc 11th Python Sci Conf 2012:51\u20135. doi:10.25080/majora-54c7f2c8-009.", "Dumka P, Sharma S, Gautam H, Mishra DR. Finite Volume Modelling of an Axisymmetric Cylindrical Fin using Python. Res Appl Therm Eng 2021;4:1\u201311."]}

In this paper, an attempt has been made to develop a Python module for evaluating the first law of thermodynamics, which includes the process of work done and amount of heat gained or lost by the system and the amount of internal energy stored. The modules NumPy and Matplotlib were used to perform the stipulated task. In addition, the correctness of codes was checked against different numerical problems, and it has been observed that the program results match exactly with the results in the literature. As a result, the functions thus developed have shown high accuracy with the least effort and error in all the cases.