Можно вместо Python сделать на маткаде или С++.
Вид задачи:
1 краевая задача для уравнения Лагранжа.
Форма и расположение области:
Шар радиуса r=1 с центром в (0,0,0).
Точное решение для задания правых частей уравнений и краевая условий задачи и оценки погрешности получаемых приближенных решений:
u(x1,x2,x3)=(x1^2+x2^2+(x3-2)^2)^(-1/2)+1