An integrable hierarchies connected with linear stationary Schrödinger equation with energy dependent potentials (in general case) are considered. Galilei-like and scaling invariance transformations are constructed. A symmetry method is applied to construct invariant solutions.