print boundaries for used variables only
issue #4: computation of the derivations using combined numerical method is done in parallel (number of processes is automatically set based on number of cores/cpus/...)