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

