Интегро-интерполяционный метод в цилиндрических координатах

[Ketsyki]

Здравствуйте. Подскажите, пожалуйста, правильные ли у меня формулы для интегро-интерполяционного метода для цилиндрической системы координат?

Модель имеет вид:

$\frac{1}{r^n} ( \frac{d}{d r} (r^n k(r) \frac{du}{dr}) - q(r) u ) = f(r)$

Граничные условия:

$u|_{r=0}$- ограничена, $-k \frac{du}{dr}|_{r=R} = \chi u|_{r=R} - \nu$

И формулы для самого метода:

для $i=1,\dots , N-1$

$$ -(r_{i+1/2} k_{i+1/2} \frac{v_{i+1} - v_i}{h_{i+1}} - r_{i-1/2}k_{i-1/2} \frac{v_i - v_{i-1}}{h_i} - \hbar_i r_i q_i v_i)=\hbar_i r_i f_i $$

для $i=N$

$$ -(-R(\chi v_i - \nu) - r_{i-1/2}k_{i-1/2} \frac{v_i - v_{i-1}}{h_i} - \hbar_i r_i q_i v_i)=\hbar_i r_i f_i $$