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