с водородоподобным атомом понятно, а вот если мы захотим рассчитать волновую функцию в случае двух и более электронов, то нам нужно учитывать взаимодействия между этими электронами, и картина будет отличаться или электроны будут иметь такие же уровни энергии, квантовые числа и тд?
Самое первое, что вам необходимо понять (похоже, вы даже не в курсе), что если для одного электрона волновая функция - это функция в 3-мерном пространстве, то для двух электронов - это функция в 6-мерном пространстве. Координаты этого пространства -

одного электрона и

другого электрона.
Это отличие - гораздо сильнее, чем впечатлившее вас отличие волновой функции без спина от волновой функции со спином.
Для трёх электронов - соответственно, 9-мерное пространство. Для

электронов -

-мерное пространство. И это ещё без учёта того, что ядро тоже квантовое.
С этим связана и проблема получить хоть какое-то решение, не только аналитическое, но и численное.
Например, пусть мы численно решаем дифференциальное уравнение (уравнение Шрёдингера - дифференциальное), разбив диапазон координаты на 100 точек (это немного, но позволит грубо увидеть график). Тогда для 3-мерной функции мы получаем

точек - миллион. Это ещё можно решить. Для двух электронов мы получаем уже триллион точек - даже чтобы хранить их, понадобятся терабайты памяти. Для трёх электронов - квинтиллион. Экзабайты. А это всего лишь атом лития. И всего лишь память - время решения уравнения обычно растёт быстрее затрат памяти.