mu = 0;
sigma = 1;
syms u;
alpha = [0.1:0.1:0.9];
k=1;
F_1 = mu + sigma * sqrt(2*pi) * erfinv((-u+alpha-u*alpha)/(u+alpha-u*alpha));
mu_k = int((F_1)^k,u,0,1);

我想求F_1的k次方的积分，运行总是出现
??? Undefined function or method 'erfinv' for input arguments of type 'sym'.
错误提示，到底该怎么改啊？
虽然知道应该是比较简单的问题，但我还是没能解决，求帮助