Оператор Лапласа

[shto-delat]

Вопрос таков.
Для повышений точности краевых условий разностной схемы требуется использовать уравнение теплопроводности в цилиндрических координатах на границе.

Известно, что в таком случае оператор Лапласа имеет вид:

$ \Delta f = {1 \over r} {\partial \over \partial r} \left( r {\partial f \over \partial r} \right) + {\partial^2f \over \partial z^2} + {1 \over r^2} {\partial^2 f \over \partial \varphi^2}$

Нужно применить этот оператор при $r = 0$.
Как это сделать?

[ewert]

Самое разумное -- использовать полуцелую сетку: $r_k={h\over2}+kh.$

[shto-delat]

ewert
Спасибо, выручили
Попробую так и сделать.