function e=stab_eq(alpha)

  % Lines that begin with % are comments.  You can explain
  % what you're doing to yourself!

  J = [0 1; -1 -alpha];
  e = eig(J);
