This script computes and plot the otoacoustic emission for a 500 Hz sinusoid, This is generated using the cochlear model described in verhulst2012 In particular otoacoustic emissions are computed as the signal difference between the sound pressure at the middle ear with model non linearities and irregularities enabled, and the sound pressure at the middle ear in case of linear model.