restart; with(LinearAlgebra): interface(rtablesize=infinity): with(linalg): ScalarMultiply(H,xi);