added Fig 3 code for Brooks et al
added complete fig 3. working on fig 4