This is the Mathematica notebook used to generate the above plot.
This table summarizes the data used to generate the above plot
| Test date | Maple | Mathematica | Sympy | # ODE’s |
| March 14, 2025 | 2024.2 (96.219%) | 14.2 (95.592%) | 1.13.3 (70.502%) | 19,466 |
| March 9, 2024 | 2024.0 (95.489%) | 14.0 (94.577%) | N/A | 15,472 |
| February 8, 2024 | 2023.2.1 (95.089%) | 14.0 (93.979%) | N/A | 13,784 |
| October 5, 2023 | 2023.1 (94.689%) | 13.3.1 (93.407%) | N/A | 10,997 |
| December 20, 2022 | 2022.2 (94.521%) | 13.2 (93.264%) | N/A | 10,258 |
| November 8, 2022 | 2022.2 (94.454%) | 13.1 (93.260%) | N/A | 10,044 |
These reports show the result of running Maple and Mathematica on my large collection of differential equations. Diagram illustrating the test system is below. Sympy was added on March 2025.
The command used for Maple is
dsolve(ode,y(x), singsol=all)
The command used for Mathematica is
DSolve[ode,y[x],x,IncludeSingularSolutions -> True]
The commands for sympy are also given in plain text for each problem in the report.
Each command was given 3 minutes of CPU time (not real time). If the command does not complete within this time it is counted as failed. in Mathematica TimeConstrained is used and in Maple timelimit is used. For sympy, timelimit is implemented using subprocess as shown here
The text books used are listed on the page of my own ODE solver at this link ODE solver.