уравнение шредингера с кубической нелинейностью

[erm]

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