Is there a way to get taylor series of expressions, containing itensors? so something like that would work: load("itensor.lisp"); taylor(rho * V([], [i]), [rho, V([], [i])], [rr, vv([], [i])], 2);