Моя электромагнитика ржавая, но у вашего вопроса есть простое аналитическое решение.
Электрический потенциал:
![enter image description here](https://i.stack.imgur.com/QCRZF.png)
и это то, на что это похожена плоскости z=0
![enter image description here](https://i.stack.imgur.com/i3yU1.jpg)
Теперь градиент равен
![enter image description here](https://i.stack.imgur.com/pQhdN.png)
и отмечается, что
![enter image description here](https://i.stack.imgur.com/mItZi.png)
Вы можете легко применить вышеупомянутое ко всем членам в уравнении градиента, чтобы получить решение в замкнутой форме, которое можно легко построить.
В MATLAB:
Вот пример, который показывает вам, как выполнить вышеупомянутое частичное дифференцирование в MATLAB.Затем вы можете опираться на это, чтобы получить полное решение.Я оставлю это на ваше усмотрение.
syms x y z x0 y0 z0
diff(1/sqrt((x-x0)^2+(y-y0)^2+(z-z0)^2),x)
ans =
-(x - x0)/((x - x0)^2 + (y - y0)^2 + (z - z0)^2)^(3/2)