Figure 2.24:

Bifurcation diagrams for the pitchfork bifurcation.

Code for Figure 2.24

Text of the GNU GPL.

main.m


1
2
3
4
5
6
7
8
npts = 51;
x = linspace(0, 1, npts)';
y = linspace(0, 1, npts)';
xy = y.^2;
ynull = zeros(size(x));

table = [-x, x, ynull, xy, y, -y, -xy];
save "pitchforkb.dat" table