function par=coefficients_PnP(xa,ra,beta,nu,Omega,zetaa,zetah,xih3,xia3,xia5,u)
xa = 0.2;
ra = 0.5;
beta = 0.2;
nu = 0.08;
Omega = 0.5;
zetaa = 0.01;
zetah = 0.01;
xih3 = 0;
xia3 = -1;
xia5 = 20;

M=[1 xa; xa ra^2];

C=[zetah+beta*u, 0; ...
    -nu*u, zetaa];

K=[Omega^2, beta*u^2;
    0, ra^2-nu*u^2];

A=[zeros(2,2), eye(2); -M\K, -M\C];
% [V,d]=eig(A,'vector');
% [autov,ind]=sort(abs(imag(d)));
% V=V(:,ind);
% par.T=[real(V(:,1)),imag(V(:,1)),real(V(:,3)),imag(V(:,3)),real(V(:,5)),imag(V(:,6))];
% W=par.T\A*par.T;
% W(abs(W)<1e-12)=0;
% par.W=W;
% par.invT=inv(par.T);
par.invM=inv(M);
par.A=A;

par.xih3=xih3;
par.xia3=xia3;
par.xia5=xia5;
par.ra=ra;
par.xa=xa;
par.Omega=Omega;
par.beta=beta;
par.nu=nu;


