Numerical solutions to differential equations