lyapunov-stability 0.1.0

Lyapunov stability analysis for agent systems — convergence, divergence, basin estimation

lyapunov-stability

Lyapunov stability analysis for agent systems. Is your fleet converging, diverging, or oscillating?

Uses real mathematics — numerical differentiation, Hessian analysis, positive definiteness checking, and basin of attraction estimation. Zero dependencies (no numpy).

Usage

from lyapunov_stability import LyapunovAnalyzer

analyzer = LyapunovAnalyzer()

# Analyze V = x^2 + y^2 (stable paraboloid)
result = analyzer.analyze(
    V=lambda x, y: x**2 + y**2,
    dVdx=lambda x, y: 2*x,
    dVdy=lambda x, y: 2*y,
    equilibrium=(0.0, 0.0),
)

print(result.stability)           # STABLE
print(result.is_positive_definite) # True
print(result.basin_radius)        # 4.99+ (large basin)

Supports 2D systems. Analytical or numerical gradients. Central difference approximation.

pip install lyapunov-stability