function t = stab221_Fromm(theta,mu)
t = stab211_Fromm(theta,mu) + 1.5*mu*mu*(1 - mu).*(exp(i.*theta)-exp(-i.*theta)).*(1 - exp(-i.*theta)).^2;

