copyright Peter Stewart, University of Glasgow 05-07-2019 files for Self-excited oscillations in a collapsible channel, with applications to retinal venous pulsation for ANZIAM Journal the files all run in MATLAB 2013a globscript_channel_trace_dimensionless.m is the global script to compute the growth rate of the mode. To do so, it first computes the static configuration using: static_twocomp_channel.m from which it then computes the growth rate using globfun_channel.m The global code uses a bisection method to compute the critical pc (isolating growth rate with zero real part) for a given pe and then takes a prescribed step in pe and restarts. This file returns a file eig_dimensionless_ten_1_rey_1.dat with 13 columns 1) dimensionless wall tension 2) Reynolds number 3) oscillation frequency 4) external pressure on compartment 1 5) external pressure on compartment 3 6) minimum height of membrane in compartment 1 7) corresponding x position to above 8) maximum height of membrane in compartment 1 9) corresponding x position to above 10) minimum height of membrane in compartment 2 11) corresponding x position to above 12) maximum height of membrane in compartment 2 13) corresponding x position to above and generates the corresponding static shape of the membrane in static_channel_dimensionless_run.dat with 3 columns 1) x positions 2) channel profile in compartment 1 3) channel profile in compartment 3