i try to use InvariantErrorPlot
command for a generic system of ODE but i find only example with pre-building equations from Mathematica Packages. This is one of them:
Needs["DifferentialEquations`NDSolveProblems`"];
Needs["DifferentialEquations`NDSolveUtilities`"];
system = GetNDSolveProblem["PerturbedKepler"]
vars = system["DependentVariables"]
H = system["Invariants"]
time = {T, 0, 290};
step = 1/25;
Met = {"SymplecticPartitionedRungeKutta", "DifferenceOrder" -> 10,
"PositionVariables" -> {Subscript[Y, 1][T], Subscript[Y, 2][T]}};
sol = NDSolve[system, time, Method -> Met, StartingStepSize -> step,
MaxSteps -> Infinity];
ParametricPlot[
Evaluate[{Subscript[Y, 1][T], Subscript[Y, 2][T]} /. sol],
Evaluate[time]]
InvariantErrorPlot[H, vars, T, sol, PlotStyle -> {Blue, Red}]
Now the question is simple: can I modify it to use InvariantErrorPlot
for my equations? And if so, how?
Thanks a lot.