code/flowRatePulseA_f.m