def plot_continuation_points(ax, lmbda_saddle, u_saddle, lmbda_hopf, u_hopf):
    """
    Plot the saddle and Hopf bifurcation points from the continuation results
    INPUT
    ax: matplotlib axes
    lmbda_saddle: array of the bifurcation parameter values corresponding to saddle points
    u_saddle: array of the infinite norms of the variable corresponding to saddle points
    lmbda_hopf: array of the bifurcation parameter values corresponding to Hopf points
    u_hopf: array of the infinite norms of the variable corresponding to Hopf points"""
    ax.plot(lmbda_saddle, u_saddle, "o", color="blue")
    ax.plot(lmbda_hopf, u_hopf, "o", color="red")

    return None