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