fix expZ -> exp(Z)
actually use compute_eps_nonzero
all four eig* at each end
eps parameter to A, A_ode, adjoint_ode
use eps parameter, rewrite eps_nonzero
pdf with eig*