Пытаюсь получить из уравнения дивергенции в декартовых координатах

уравнение дивергенции в цилиндрических координатах

Получаю в формуле (D2) лишне слагаемое

.
Использую для этого две формулы


и формулы


которые получил из формул


из условия

Подставляя формулы (1a), (1b), и (2a), (2b) в уравнение (D1) получаю

Использование формул (1a), (1b), и (3a), (3b) позволило мне получить уравнение градиента в цилиндрических координатах.
Поэтому подозреваю ошибку в формулах (2a), (2b).
И где же у меня ошибка?
Заранее спасибо.