绘制复函数

Question

我试图绘制以下函数，但没有结果。这是一个非常复杂的函数，我真的不知道这是原因还是一个愚蠢的错误，因为我不太熟悉Wolfram。任何关于如何用其他语言解决这个问题的建议都是非常受欢迎的。你知道吗

f[y_, n_] := n*exp (-n^2/y)/(sqrt (Pi)*y^((n + 3)/2))*(HermiteH[n, (y + 2 n)/(2*sqrt (y))] - 
    sqrt (y)*HermiteH[n - 1, (y + 2 n)/(2*sqrt (y))])

Plot[f[y, 500], {y, 0, 6000}]

谢谢

Answer 1

当我在第12版中尝试Gilles Philippe Paillé的例子时，我得到了

General munfl 4.697873421421826967046763728`16.*^-9856 is too small to represent as a
normalized machine number, precision may be lost.

一次又一次的重复，接着是一个看起来像是空白的情节

试着用这个来解决问题

f[y_, n_] := n*Exp[-n^2/y]/(Sqrt[Pi]*y^((n + 3)/2))*(HermiteH[n, (y + 2 )/
  (2*Sqrt[y])] - Sqrt[y]*HermiteH[n - 1, (y + 2 n)/(2*Sqrt[y])]);
ListPlot[Table[f[y, 500], {y, 10, 6000,10}],PlotRange->All,Joined->True]

它仍然会给我关于数字大小的警告，但随后会显示绘图。（我从10开始计算，以避免被0除）

Answer 2

考虑到你问题的标签，你混合了Python和Mathematica。以下代码适用于Mathematica 10.4（参见Bill对版本12的回答）：

f[y_, n_] := n*Exp[-n^2/y]/(Sqrt[Pi]*y^((n + 3)/2))*(HermiteH[n, (y + 2 n)/(2*Sqrt[y])] -
    Sqrt[y]*HermiteH[n - 1, (y + 2 n)/(2*Sqrt[y])])
Plot[f[y, 500], {y, 0, 6000}, PlotRange -> All]