Identification of working fluids and development of their mathematical models should always precede construction of a proper model of the analysed thermodynamic system. This paper presents method of development of a mathematical model of working fluids in a gas turbine system and its implementation in Python programming environment. Among the thermodynamic parameters of the quantitative analysis of systems, the following were selected: specific volume, specific isobaric and isochoric heat capacity and their ratio, specific enthalpy and specific entropy. The development of the model began with implementation of dependencies describing the semi-ideal gas. The model was then extended to the real gas model using correction factors reflecting the impact of pressure. The real gas equations of state were chosen, namely due to Redlich–Kwong, Peng–Robinson, Soave– Redlich–Kwong, and Lee–Kesler. All the correction functions were derived analytically from the mentioned equations of real gas behaviour. The philosophy of construction of computational algorithms was presented and relevant calculation and numerical algorithms were discussed. Created software allowed to obtain results which were analysed and partially validated.