使用 Mathematica 7 调试 Mathematica 5 上的工作程序

Question

我目前正在阅读 Mathematica Guidebooks for Programming 并且我正在尝试制定这本书的第一个程序。基本上，当我运行以下程序时:

Plot3D[{Re[Exp[1/(x + I y)]]}, {x, -0.02, 0.022}, {y, -0.04, 0.042}, 
 PlotRange -> {-1, 8}, PlotPoints -> 120, Mesh -> False, 
 ColorFunction -> Function[{x1, x2, x3}, Hue[Arg[Exp[1/(x1 + I x2)]]]]]

要么我得到 1/0 错误和 e^\infinity 错误，或者，如果我降低 PlotPoints选项，例如 60，溢出错误。虽然我有一个工作输出，但它不是它应该的样子。色调似乎从左角扩散，而它应该从原点扩散(如原始 output 所示)

这是显然在 Mathematica 5 上运行的原始程序(Trott，Mathematica 编程指南):

Off[Plot3D::gval];

Plot3D[{Re[Exp[1/(x + I y)]], Hue[Arg[Exp[1/(x + I y)]]]},
       {x, -0.02, 0.022}, {y, -0.04, 0.042},
       PlotRange -> {-1, 8}, PlotPoints -> 120, Mesh -> False]

Off[Plot3D::gval];

然而， ColorFunction使用这种方式(第一个 Plot3D 参数)不起作用，所以我试图简单地适应它的新使用方式。

嗯，谢谢我猜!

Answer 1

如果您对 Mathematica 的默认设置感到满意，您可以使用旧版本的代码，只需删除 , Hue[Arg[Exp[1/(x + I y)]]]并且该功能工作正常。

您在新版本代码中遇到的问题似乎源于表达式 Exp[1/(x1 + I x2)] -- 有时这需要评估 1/0 .至少，如果我删掉 1/该程序执行(在 Mathematica 7 上)没有任何提示，但显然颜色错误。所以你可能需要重写你的颜色函数。