par.a = 75e3/2; % radius [m]
par.h = 1.004e3; % thickness [m]
par.E = 73.2e9; % Young's modulus [Pa]
par.nu = 0.155; % Poisson's ratio
par.rho = 2202; % density [kg/m^3]
%Calculate fundamental modes of the disk
[freqs, modes, shapes, x, y] = disk_frequencies(par, 10000, 1, 'shapes', 0.5e3);
%Now we extract the force profile from the COMSOL model
... 27 more lines ...
