fragilityIDTFexample2.m

function Omegal=fragilityIDTFexample2(kp,kd,w,mu)
Omegal=(1/2).*(1+4.*w.^2).^(-3/4).*(1+9.*w.^2).^(-5/4).*csc((1/2).*mu.*pi).*(( ...
  -1).*kp+(-1).*(1+9.*w.^2).^(-1/4).*(w.*cos((1/2).*(mu.*pi+atan(3.*w)))+( ...
  1+4.*w.^2).^(1/4).*cos((1/2).*(mu.*pi+atan(3.*w))).*sin((1/2).*atan(2.* ...
  w))+(-1).*(1+4.*w.^2).^(1/4).*cos((1/2).*atan(2.*w)).*sin((1/2).*(mu.* ...
  pi+atan(3.*w)))).*(cos((1/2).*(mu.*pi+atan(3.*w))).*sin((1/2).*atan(3.* ...
  w))+(-1).*cos((1/2).*atan(3.*w)).*sin((1/2).*(mu.*pi+atan(3.*w)))).^(-1) ...
  ).*((1+4.*w.^2).^(3/4).*(2+9.*w.^2).*cos((1/2).*(mu.*pi+atan(3.*w)))+(( ...
  -1)+6.*w.^2).*cos((1/2).*(mu.*pi+(-1).*atan(2.*w)+atan(3.*w)))+(-3).*w.* ...
  (1+4.*w.^2).^(3/4).*sin((1/2).*(mu.*pi+atan(3.*w)))+5.*w.*sin((1/2).*( ...
  mu.*pi+(-1).*atan(2.*w)+atan(3.*w))))+(1/2).*w.^((-1)+(-1).*mu).*(1+4.* ...
  w.^2).^(-3/4).*(1+9.*w.^2).^(-5/4).*csc((1/2).*mu.*pi).*((w+(-6).*w.^3) ...
  .*cos((1/2).*(atan(2.*w)+(-1).*atan(3.*w)))+w.*(1+4.*w.^2).^(3/4).*(2.*( ...
  (-1)+mu)+9.*((-1)+2.*mu).*w.^2).*cos((1/2).*atan(3.*w))+(2.*mu+w.^2.*(5+ ...
  26.*mu+72.*mu.*w.^2)).*sin((1/2).*(atan(2.*w)+(-1).*atan(3.*w)))+3.* ...
  w.^2.*(1+4.*w.^2).^(3/4).*sin((1/2).*atan(3.*w))).*((-1).*kd+(-1).*w.^(( ...
  -1).*mu).*(1+9.*w.^2).^(-1/4).*sec((1/2).*(mu.*pi+atan(3.*w))).*(w.*cos( ...
  (1/2).*atan(3.*w))+(1+4.*w.^2).^(1/4).*cos((1/2).*atan(3.*w)).*sin((1/2) ...
  .*atan(2.*w))+(-1).*(1+4.*w.^2).^(1/4).*cos((1/2).*atan(2.*w)).*sin(( ...
  1/2).*atan(3.*w))).*((-1).*sin((1/2).*atan(3.*w))+cos((1/2).*atan(3.*w)) ...
  .*tan((1/2).*(mu.*pi+atan(3.*w)))).^(-1));
end