Computation of normal form coefficients of cycle bifurcations of maps by algorithmic differentiation

Computation of the normal form coefficients confirms that AD is as accurate as SD. Moreover, elapsed time in computations using AD grows linearly with the iteration number J, but more like Jd for d th derivatives with SD. For small J, SD is still faster than AD.