For convenience, the following code should work out of the box:
% computation of primary branch (EP)
bd1 = coco('1', 'ode', 'isol', 'ep', ... % branch, TBF, IPT, BT
@brus, [1;0], {'A' 'B'}, [1; 0], ... % F, x0, PNames, p0
'B', [0 3]); % 'continuation parameter, computational domain
% computation emerging periodic orbits (PO)
HB = coco_bd_labs(bd1, 'HB');
bd2 = coco('2', 'ode', 'HB', 'po', ... % prob.struct., branch, TBF, IPT, BT
'1', HB(1), ... % primary branch name, solution label
{'B'}, {[0 3]}); % 'continuation parameter, computational domain
Having cont.NAdapt=1 removes the preliminary termination of the computation.
Due to backward compatibility, mesh adaptation is disabled by default and must be enabled by the use of the coco_set command. Error estimation for trajectory segments and termination when the error estimate exceeds a threshold is enabled by default. The latter may be disabled by the use of the coco_set command. As a consequence, continuation of solutions to boundary-value problems, e.g., periodic orbits from a Hopf bifurcation, may result in premature termination at an MXCL point, corresponding to a solution for which the error estimate exceeds the threshold.