An asymptotic numerical method for continuation of spatial equilibria of special Cosserat rods