collocation_fun(fun, y, p, x, h)
This function lies in the core of the method. The solution is sought as a cubic C1 continuous spline with derivatives matching the ODE rhs at given nodes :None:None:`x`. Collocation conditions are formed from the equality of the spline derivatives and rhs of the ODE system in the middle points between nodes.
Such method is classified to Lobbato IIIA family in ODE literature. Refer to for the formula and some discussion.
Collocation residuals at the middle points of the mesh intervals.
Values of the cubic spline evaluated at the middle points of the mesh intervals.
RHS of the ODE system evaluated at the mesh nodes.
RHS of the ODE system evaluated at the middle points of the mesh intervals (and using :None:None:`y_middle`).
Evaluate collocation residuals.
Hover to see nodes names; edges to Self not shown, Caped at 50 nodes.
Using a canvas is more power efficient and can get hundred of nodes ; but does not allow hyperlinks; , arrows or text (beyond on hover)
SVG is more flexible but power hungry; and does not scale well to 50 + nodes.
All aboves nodes referred to, (or are referred from) current nodes; Edges from Self to other have been omitted (or all nodes would be connected to the central node "self" which is not useful). Nodes are colored by the library they belong to, and scaled with the number of references pointing them