В рамках одноэлектронной задачи для решения системы разностных уравнений относительно коэффициентов разложения волновой функции по атомным орбиталям (амплитуд атомных функций) предложен метод, явно не опирающийся на технику функций Грина. Метод содержит в себе одновременно идею кластерных и зонных расчетов и дает существенные преимущества при анализе локализованных состояний, что показано на примере описания как одиночных примесей, так и их комплексов в ионных кристаллах со структурой флюорита. Рассмотрен случай периодического распределения дефектов.