The location of Hopf and limit points usually does not cause problems. However,
the location of branching points can give problems. The region of attraction of the
Newton type continuation method which is used, has the shape of a cone (see [1]).
In the localisation process we cannot assume to stay in this cone.
This difficulty can be avoided by introducing p Î Rn and b Î R and
considering the extended system:
| (41) |
We solve this system by Newton's method with initial data b = 0 and p: fuTp=mp where m is the real eigenvalue with smallest norm. A branching point (u,a) corresponds to a regular solution (u,a,0,p) of system (41) (see [2],p. 165). We note that the second order partial derivatives (Hessian) of f with respect to u and a are required.
The tangent vector at the singularity is also computed here. This is related to the processing of the branching point (computing the direction of the secondary branch).