急求大神！！！matlab中erf函数精度提高。

matlab代码如下

k=1079; z=3^0.5 ；z1=sqrt(-i*k)*(1+z);

mfun('erf',z1)

采用mfun('erf',x)，在matlab计算只有16位有效数字（和mathematica中得到的准确结果相比），那么如何才能在matlab中把计算精度提高到30位甚至更高的有效数字？现在很着急，做了很久都没成功，急求各位大神指点。

Answer 1

[kendin]

方法写给你了，照葫芦画瓢，用符号工具箱计算，最后再用vpa做符号转为数值的显示。

或者直接用mathematica算，普通的读写文件以及循环之类简单功能mathematica同样能实现。