请高手帮我解一个数值积分，谢谢：）

[angel4296]

function Y1=chao(x)
Y1 =27671378696711990071590590550091/3749001307303900038450491949056*x^(19081/10000)*besselk(16719/5000,1/562949953421312*7177179097769788095511012573184^(1/2)*x^(1/2))*(-pi^(1/2)*erf(10*500/707*x)+pi^(1/2))

把上面的保存成chao.m的文件，然后对chao求数值积分，在0，1区间，大家知道怎么解吗？我用quad，总是提示我：

??? Error using ==> mpower
Matrix must be square.

Error in ==> chao at 2
Y1 =27671378696711990071590590550091/3749001307303900038450491949056*x^(19081/10000)*besselk(16719/5000,1/562949953421312*7177179097769788095511012573184^(1/2)*x^(1/2))*(-pi^(1/2)*erf(10*500/707*x)+p
Error in ==> quad at 63
y = f(x, varargin{:});