Wastewater data, including Population of the Municipality Served by the Wastewater Treatment Plants (WWTPs), Capacity, Number of WWTPs, Amount of Daily Wastewater Discharged per Person and Wastewater Treated in WWTPs (y(t)) obtained from Turkish Statistical Institute (TUIK) for Kayseri province for 2003-2020 were mathematically modelled and analysed with Artificial Neural Network (ANN) and Differential equations. Firstly, the data were augmented through ARIMA since TUIK data numbers are insufficient to ANN training and later the data were normalized. The augmented normalized data was trained with ANN twice, thus the effect of other variables on the y(t) variable was shown and mathematical ANN activation functions in the form of a tangent hyperbolic function was proposed for this variable. Then, arbitrary parameters used for a linear system consisting of differential equations representing the five variables mentioned above were estimated using normalized original data and thus the ODE (ordinary differential equation) model was proposed. Two ANN models and ODE model were evaluated on normalized real TUIK data and the performances of these three models were compared. Among these mathematical models, the model that gave the minimum MSE (mean squared error) has been determined as the ODE model. Finally, future predictions were made for the y(t) variable with the ODE model.