function t = stab100(theta, mu)
t = 1 - 3*mu*(1 - exp(-i.*theta));