PeriodicOrbit.jl
abstract type PeriodicOrbit end
include("SimpleShooting.jl")
include("PeriodicOrbitFD.jl")
function Flow(Fper, dFper, x0, M::Int, options::NewtonPar)
options.verbose = false
xnew = similar(x0)
for ii=1:M
xnew, hist, flag = newton(
u->Fper(u, x0),
u->dFper(u, x0),
x0,
options)
@assert flag == true "Newton method for computing the Flow did not converge at i = $ii, hist = $hist, norm(x0) = $(norm(x0, Inf64))"
x0 .= xnew
end
return xnew
end