function t = stab211(theta,mu)
t = stab111(theta,mu) + 1.5*mu*(1 - mu)*(1 - exp(-i.*theta)).*(1 - exp(i.*theta));
