## Determine Hill's equation stability regions
A = linspace(0, 1, 10);
E = linspace(0, .5, 12);

[AA, EE] = meshgrid(A, E);
Z = hillstab(A, E);

mesh(AA, EE, Z);

plot([Z, 2 * ones(length(ans), 1), -2 * ones(length(ans), 1)]);

## ;;; Local Variables: ***
## ;;; mode: octave ***
## ;;; End: ***
