MultiPrecision<Pow2.N8> mu = 1.5; DormandPrinceAdaptiveODESolver<Pow2.N8> solver = new( v: (0, 1), (x, y) => (y, mu * (1 - x * x) * y - x), abstol: "1e-20", reltol: "1e-20", maxdepth: 12 ); for (int i = 0; i <= 4096; i++) { Console.WriteLine($"{i / 256d},{solver.X},{solver.Y}"); solver.Next(1d / 256); } Assert.IsTrue(solver.MaxError < 1e-14);