I have the system natural frequency and zeta, need to solve for stiffness and damping matrices

I am scaling down a car to 1/10th mass and 1/5th dimensions. I have already conducted an analysis of the full scale car and have determined the natural frequencies, modeshapes and zeta. I now need to replicate these values for the scaled down version of my car.

I have the eigenvalues from the eig function for the full scale:

 -23.8950 + 80.2990i
 -23.8950 - 80.2990i
 -24.3328 + 78.5293i
 -24.3328 - 78.5293i
 -24.1746 + 79.0562i
 -24.1746 - 79.0562i
 -24.4534 + 79.4662i
 -24.4534 - 79.4662i
  -1.2450 + 7.8696i
  -1.2450 - 7.8696i
  -1.6017 + 8.5027i
  -1.6017 - 8.5027i
  -1.6518 + 9.4437i
  -1.6518 - 9.4437i

I would like to sub these into the eigenvalue function to solve for kf,kt,kr and cs. My matlab skills in this area are somewhat lacking, so any help would be greatly appreciated..

if true
%Known Rig Dimensions
%======================================================================
n=7;
lf=0.34; lr=0.37;tf=0.14;tb=0.14;
% tf=distance from COG to front axle; tb=distance from COG to rear axle
L=[1 1 1 1
  lf lf -lr -lr
  tf -tf -tb tb];
Ms=[180 7.8 8.85]; %Sprung mass Kg pitch moment of inertia(kg-m^2), roll moment of inertia (kg-m^2)
Ms=diag(Ms);
Mu=4.2*ones(1,4); 
Mu=diag(Mu); % Unsprung Mass
%Unknown Stiffness and Damping Values
%=====================================================================
Bs=[cs cs cs cs]; %Damping values N.s/m
Bs=diag(Bs);
Ks=[kf kf kr kr]; %Spring Stiffness
Ks=diag(Ks);
Kt = [kt kt kt kt]; % Tyre Stiffness
Kt=diag(Kt); % N/m
Km=[Ks+Kt -Ks*L';-L*Ks L*Ks*L']; % without anti-roll bar
Bm=[Bs -Bs*L';-L*Bs L*Bs*L'];
Mm=[Mu zeros(4,3); zeros(3,4) Ms]; 
%Attempt to solve for cs,kf,kr,kt
%=====================================================================
EigA=[ -23.8950 + 80.2990i
-23.8950 - 80.2990i
-24.3328 + 78.5293i
-24.3328 - 78.5293i
-24.1746 + 79.0562i
-24.1746 - 79.0562i
-24.4534 + 79.4662i
-24.4534 - 79.4662i
-1.2450 + 7.8696i
-1.2450 - 7.8696i
-1.6017 + 8.5027i
-1.6017 - 8.5027i
-1.6518 + 9.4437i
-1.6518 - 9.4437i];
A=[zeros(n) eye(n);
  -Mm\Km -Mm\Bm];
Solve(eig(A)==EigA,kr,kf,cs,kt)
end