%pt = floor(length(branch1.point)/2);
pt = 1;
w0 = first_hopf.omega;
second_hopf.omega = w0;

while abs(second_hopf.omega - w0) < 0.1 && pt < length(branch1.point)
    pt = pt + 1;
    hopf = p_tohopf(branch1.point(pt));
    method=df_mthod('hopf');
    [hopf,success]=p_correc(hopf,1,[],method.point);
    second_hopf = hopf;
    hopf.stability=p_stabil(hopf,method.stability);
end
