Но ведь это в декартовых координатах. В полярных координатах у меня получилось громоздкое выражение, а хотелось бы увидеть что-то попроще.

. Так устроит? Как Вы записываете

в полярных координатах? Приведите Ваши выкладки.
В этих уравнениях я сделал замену

,

и в итоге получил громоздкую систему уравнений в показательно-полярных координатах. Но оказывается, что если бы я применил такую систему координат, в которой градиенты были бы ортонормальны, то система уравнений Коши-Римана не изменилась бы. Полярные координаты ортогональны, но не ортонормальны, поэтому хотелось бы нормализовать их. Есть у Вас готовое решение?