

  %
  clear all
  tp=1;
  r=1;
  NU=linspace(0,10,1001);
  for ia=1:length(NU)
      nu=NU(ia);
  
  A1=[r-nu,nu;nu,-r-nu];
  A2=[-r-nu,nu;nu,r-nu];
  AE1=expm(A1*tp/2);AE2=expm(A2*tp/2);
  Phi=AE2*AE1;
  mu=max(abs(eig(Phi)));
  lambda(ia)=log(mu)/tp;
  end

  figure(45)
  plot(NU,lambda,'r','linewidth',2)
  xlabel('$\nu$','Interpreter','latex')
  ylabel('$\lambda$','Interpreter','latex')
  grid on
  set(gca,'Fontsize',24)



