function t = stab222_Fromm(theta,mu)
t = stab221_Fromm(theta,mu) - 0.75*mu*mu*mu*(1-mu).*(exp(i*theta)-exp(-i.*theta)).*(1-exp(-i*theta)).^3;


