Page 67 -

FIGURE 2.5
Plot of multiple Hénon orbits having the same a = 0.25 but random starting points.

for n=1:1500
for m=1:40
x(1,m)=-0.99+2*rx(m);
y(1,m)=-0.99+2*ry(m);

x(n+1,m)=a*x(n,m)-b*(y(n,m)-(x(n,m))^2);
y(n+1,m)=b*x(n,m)+a*(y(n,m)-(x(n,m))^2);
end
end
plot(x,y,'r.')
axis([-1 1 -1 1])
axis square

2.9 Generation of Special Functions from Their Recursion
Relations*
In this section, we go back to more classical mathematics. We consider the
case of the special functions of mathematical physics. In this case, we need to

© 2001 by CRC Press LLC