function t = stab222_BW(theta,mu)
t = stab221_BW(theta,mu) - 1.5*mu*mu*mu*(1-mu).*exp(-i.*theta).*(exp(i*theta)-1).*(1-exp(-i*theta)).^3;


