Математическое моделирование двумерной задачи

Нужно написать программу в Matlab для разностной схемы. сигма r и сигма z принять пока за какие-то матрицы. кси это логарифм от r, z это ось y сетки. Нужно найти U по сетке. сигма r это электропроводимость по горизонтали. сигма z это электропроводимость по вертикали. электропроводимость это 1 делить на удельное электрическое сопротивление. формулы расчета электрического сопротивления могу скинуть, если решите вместе с формулами и данными полностью написать программу.