maxima - 无法画出电场偶极子

Question

我正在尝试使用 Maxima 来模拟电场偶极子。 这是代码:

P : 1$
eps0 : 8.854$
Ex : 3*P*y*x/(4*π*eps0*(x^2 + y^2)^2.5)$
Ey : P*(3*y^2/(x^2 + y^2) + 1)/(x^2 + y^2)^1.5/(4*π*eps0)$
ew : sqrt(Ex^2 + Ey^2)$

contour_plot (ew,[x, -4, 4], [y, -4,4]);

但我只得到直线。
可能有什么问题，因为 Wolfram Mathematica 中的类似代码可以正常工作。

Answer 1

contour_plot 并不像 Mathematica 中的相应函数那么智能。你可以帮忙一点。

P : 1$
eps0 : 8.854$
Ex : 3*P*y*x/(4*%pi*eps0*(x^2 + y^2)^2.5)$
Ey : P*(3*y^2/(x^2 + y^2) + 1)/(x^2 + y^2)^1.5/(4*%pi*eps0)$
ew : sqrt(Ex^2 + Ey^2)$

BIG: subst([x=1/4, y=1/4], ew)$
cap_log(e, c):= if e > c then log(c) else log(e)$

set_plot_option ([gnuplot_preamble, "set cntrparam levels 8"])$
contour_plot ('cap_log(ew, BIG), [x, -4, 4], [y, -4, 4])$