function t = stab210_BW(theta,mu)
t = stab110(theta,mu) - 1.5*mu*(1 - mu)*((1 - exp(-i.*theta)) - (exp(-i.*theta) - exp(-2*i.*theta)));

