Verification of Taylor's theorem

Multivariate function calculus is an important part of mathematical analysis courses, and most conclusions can be found and generalized in univariate calculus. However, the biggest difficulty in teaching multivariate calculus lies in its abstraction, such as Taylors theorem, multiple integral regions drawing, and integral variable transformation. At the same time, ordinary differential equations are also one of the basic courses of the profession, and dynamic systems based on ordinary differential equations have extensive applications in mathematical models of continuity problems and optimal control problems. Software such as Mathematica, Python, Matlab, etc. can solve similar problems. Therefore, this article will use the visualization and computational capabilities of Mathematica to validate important definitions and conclusions in multivariate calculus, and compare the differences among the three software in solving approximate numerical solutions of dynamic systems of ordinary differential equations from different perspectives.