Для вычисления производных координатных функций поверхности методом обобщенного дифференцирования получен прямоугольный конечный элемент оператора частной производной. Приведены результаты тестирования алгоритма численного дифференцирования для двух модельных поверхностей, заданных графиками аналитических функций двух переменных. Для сравнения тестовые задачи решены методом конечных разностей. Результаты иллюстрируют повышенную точность обобщенного численного дифференцирования.