F=[0.001:0.001:1];
z=F.^(-2/3);
figure(1); clf;

plot(F,z,'r-','linewidth',2)
axis([0 1 0 10])
hold on 
y=F.^(1/3);
plot(F,5*y,'k--','linewidth',1.5)
set(gca, 'fontsize',12,'fontweight','bold');
xlabel('Force amplitude (a.u.)','fontsize',13,'fontweight','bold');
ylabel('Gain Hopf bifurcation','fontsize',13,'fontweight','bold');
text(0.4,4,'Output driven Hopf bifurcation: F^{(1/3)}','fontsize',13,'fontweight','bold', 'rotation',10);
text(0.4,2,'Gain','fontsize',13,'fontweight','bold', 'rotation',-12);
text(1.2,5,'1.0','fontsize',12,'fontweight','bold');

print -djpeg Fig5-24