求助！关于多元函数积分的问题

yousipu

最近在编一个程序，一部分代码如下：
syms x u d r n
F=1/((1+x)*sqrt(2*pi*d^2))*exp(-(log(1+x)+d^2/2)^2/(2*d^2));   %对数正态分布的分布函数
f=1/(x*sqrt(2*pi*d^2))*exp(-1/2*((log(x)-u)/d)^2);    %对数正态分布的密度函数
frn=(vpa('n!')/(vpa('(r-1)!')*vpa('(n-r)!')))*F^(r-1)*(1-F)^(n-r)*f;  %次序统计量密度函数
Ern=int(x*frn,x,0,1)

其实Ern就是求次序统计量的期望，但是估计式子太复杂，int积分不出来，matlab提示：
Warning: Explicit integral could not be found.

然后我考虑用数值积分函数quadl(),将上面int函数一行换成
[I,num]=quadl(frn,0,1)
但也不对，matlab提示：
??? Error using ==> fcnchk
If FUN is a MATLAB object, it must have an feval method.
Error in ==> quadl at 54
f = fcnchk(funfcn);
Error in ==> tmp3 at 11
[I,num]=quadl(frn,0,1)

是我对quadl使用错误吗？请问quadl能对多元函数数值积分吗？这个问题该如何解决？恳请高人指点，我在网上查了半天没弄明白，谢谢谢谢谢！！！
我QQ 93660036，在线等回复