Next: Continuation of limit cycles Up: Bifurcations Previous: Bifurcations   Contents


Branching point locator

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:

ì
ï
ï
í
ï
ï
î
f(u,a) + bp
=
0
fuT(u,a)p
=
0
pTfa(u,a)
=
0
pTp-1
=
0
(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).