demo2.m

prob = coco_prob();
prob = coco_add_func(prob, 'phi', @phi,... 
  [], 'zero', 'u0', [7.9829e-02   2.5133e+01  -2.0000e+00]);
prob = coco_add_func(prob, 'psi', @psi,... 
  [], 'inactive', 'mu1', 'uidx', [1:3]);
prob = coco_add_func(prob, 'iden1', @iden,... 
  [], 'active', 'mu2', 'uidx', [2]);
prob = coco_add_func(prob, 'iden2', @iden,...
  [], 'active', 'mu3', 'uidx', [3]);
%prob = coco_set_parival(prob, 'mu1', 0);
prob = coco_add_pars(prob, '' , 1:3, ...
  {'x' 'y' 'z'});
coco(prob, 'run1', [], 1, {'x' 'y' 'z' 'mu1'})