Modify RTACalc for increasing the number of layer

I am using the script RTA calc for modelling multilayer behavior.

The script has been designed for 2 layers, but I would like to raise it to 3.

%%%

layers = {'fused silica' 'Pt' 'SiO2' 'Ag' 'air'};

film1thicknesses = [500]; % nm Pt

film2thicknesses = [1]; % nm SiO2

film3thicknesses = [30]; % nm Ag

%%%%

bilayerDispersionStudy(lam0, lam1, dlam, date, suffix, layers, film1thicknesses, film2thicknesses, film3thicknesses, angle, polarization, single, control);

clearvars date suffix lam0 lam1 dlam layers film1thicknesses film2thicknesses polarization control single;

However, the matrix that has been originally set is for two layers, so I get a series of mistakes