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


