Buckling column

Model

using ReachabilityAnalysis

@taylorize function buckling_column!(dx, x, params, t)
    dx[1] = x[2]
    dx[2] = 2 * x[1] - x[1]^3 - 0.2 * x[2] + 0.1
    return dx
end

References