In this paper, we provide a mathematical and statistical methodology using heteroscedastic estimation to achieve the aim of building a more precise mathematical model for complex financial data. Considering a general regression model with explanatory variables (the expected value model form) and the error term (including heteroscedasticity), the optimal expected value and heteroscedastic model forms are investigated by linear, nonlinear, curvilinear, and composition function forms, using the minimum mean-squared error criterion to show the precision of the methodology. After combining the two optimal models, the fitted values of the financial data are more precise than the linear regression model in the literature and also show the fitted model forms in the example of Taiwan stock price index futures that has three cases: (1) before COVID-19, (2) during COVID-19, and (3) the entire observation time period. The fitted mathematical models can apparently show how COVID-19 affects the return rates of Taiwan stock price index futures. Furthermore, the fitted heteroscedastic models also show how COVID-19 influences the fluctuations of the return rates of Taiwan stock price index futures. This methodology will contribute to the probability of building algorithms for computing and predicting financial data based on mathematical model form outcomes and assist model comparisons after adding new data to a database.