a=2;
b=9;
syms x
I0=besseli(0,a*x/b);
F(x)=(x/b)*exp(-(x^2+a^2)/(2*b));
FUN=x*F(x);
mean=quad(FUN,0,100)
我得到这个错误:
Error using fcnchk (line 107)
If `FUN` is a MATLAB object, it must have an feval method.
Error in quad (line 57)
f = fcnchk(funfcn);