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