a = 0; b = 10; ContourPlot[{y + 3 ArcTan[y/(x - a)]}/20, {x, -b,
  b}, {y, -b, b}, Exclusions -> {x == a}, PlotPoints -> 50,
ContourStyle -> {Red}, Frame -> True, Axes -> True,
ContourShading -> None, PlotLegends -> {"\[Psi](x,y)"}]

6楼: Originally posted by xzczd at 2015-01-08 10:26:41
……跟Log没关系，这是极端条件下的数值误差。（楼主你很有才能啊，这么短时间内能踩出这么多雷。）嗯，有意思，WorkingPrecision选项对这种情况似乎无效。对于你这个问题，最简单的方法，就是别除那个300。除不除 ...

a = 0;
Limit[y + 3 ArcTan[y/(x - a)], x -> 0, Direction -> #] & /@ {-1, 1}
(* {y + (3 \[Pi] Sqrt[y^2])/(2 y), y - (3 \[Pi] Sqrt[y^2])/(2 y)} *)

a = 0; b = 10;
ContourPlot[{y + 3 Arg[(x - a) + I y]}/20, {x, -b, b}, {y, -b, b},
  ContourStyle -> Red, Frame -> True, Axes -> True,
  ContourShading -> None] /.
Graphics[a__] :>
  Legended[Graphics[a],
   Placed[LineLegend[{Red}, {"\[Psi](x,y)"}], After]]

2楼: Originally posted by xzczd at 2015-01-07 22:28:07
很遗憾，如果你坚持要用这个方程，那么这就是个不可能的任务，因为ContourPlot忠实地画出了你输入的东西：

a = 0;
Limit & /@ {-1, 1}
(* {y + (3 \ Sqrt)/(2 y), y - (3 \ Sqrt)/(2 y)} *)


如你所见 ...

2楼: Originally posted by xzczd at 2015-01-07 22:28:07
很遗憾，如果你坚持要用这个方程，那么这就是个不可能的任务，因为ContourPlot忠实地画出了你输入的东西：

a = 0;
Limit & /@ {-1, 1}
(* {y + (3 \ Sqrt)/(2 y), y - (3 \ Sqrt)/(2 y)} *)


如你所见 ...

z1 = -(0.5*Log[x^2 + y^2] - Arg[x + y I])/(2 Pi);
ContourPlot[{z1}/300, {x, -10, 10}, {y, -10, 10},
ContourStyle -> {Red}, Frame -> True, Axes -> False,
ContourShading -> None,
Exclusions -> {x == 0, y == 0}, PlotPoints -> 50]