=try to reproduce figure 4 of Voigt+12