用matlab求欧拉常数代码，谢谢各位

如何用matlab求欧拉常数？

1、首先我们根据欧拉常数的定义，写出其表达式，如下图所示。

2、从表达式我们看到，求和部分可以用symsum函数来求解

3、然后再用limit函数，求其n一﹥∞的极限

4、完整的代码如下

>>syms k n

>>S = symsum(1/k,k,1,n) - log(n)

>>vpa(limit(S,n,Inf),20)

5、也可以直接用下列命令来求解

>>-psi(1)

6、执行结果