calculate subcortical activation

Adds an example to examples/inverse using the spm faces dataset. Calculates the activation for a subcortical volume. 
t is monotonic, so `set_xlim(t[0, -1])` should work