The code of the particular functions is not commented as it is fairly primitive. 

pshift.m -- Calculates sagnac phase shift with given modulation frequency fm, rotation rate in rad/s om, modulation index h and time t. Be sure to set up parameters of the gyroscope  c, L, r, n ,lam
spektr.m -- Calculates power of the output at the modulation frequency. Uses pshift.m. Takes fm -- modulation frequency,rot -- rotation rate, h -- modulation index
slope.m -- Calculates slope at particular frequency and index
slopeplot.m -- Calculates modulation index dependency
modsim.m -- Can be used to simulate output under various modulation frequency and modulation indices

Essentially, set up gyro parameters at pshift.m, then run either slopeplot.m for or modsim.m
