script for plotting gamma and Bk as functions of disturbance frequency and magnitude