clc
clear all

modi_considerati=15;


freq=[9.170855E-05
9.220868E-05
9.384325E-05
6.455951E-04
7.482965E-04
8.344773E-04
7.329947E+01
2.072568E+02
2.173950E+02
2.269527E+02
2.331142E+02
3.516379E+02
3.628181E+02
3.875460E+02
4.057810E+02
];

freq=freq(1:modi_considerati);

lambda=2*pi*freq.^2;


num_rigid_modes=6;

PHI_struct = load('Phi_finale.mat');                        % permette di caricare il .mat in una struttura
Phi = PHI_struct.PHI_finale;                                % estrae la matrice vera e propria da utilizzare

Calcolo_Fattori_di_partecipazione(Phi, lambda, num_rigid_modes);

 