功能:求解Whittaker函数

格式: WhittakerW(a,b,z)
a : 矩阵变量或者数值
b : 矩阵变量或者数值
z : 矩阵变量或者数值

说明:
1. 当a,b,z任意两个为矩阵变量时,其维度必须一致。执行结果返回同大小矩阵变量.
2. 当参数不合理时,将返回非数字。
3. 原理
$$
\mathbf{WhittakerW}(a,b,z) = \dfrac{\exp(-z/2)z^a}{\Gamma(1/2-a+b)}\int_0^\infty t^{-a-1/2-b}(1+t/2)^{a-1/2+b}\exp(-t)dt
$$

原理: {Whittaker函数<https://mathworld.wolfram.com/WhittakerFunction.html>}

例子:
a = WhittakerW(0.4,1.2,0.58)//回车得到如下结果
a =
[ 2.23235838694637 ]