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