see also JuliaDynamics/DynamicalSystems.jl#210
Would be nice to allow a 4th argument εy for the recurrences of y.
Only question is: what would be the recurrence type stored and displayed by the joint matrix? At the moment we store ε and we use it when we display the matrix. I guess we could make a new subtype called only for display purposes.