请教该积分为何不能求解

simonwoo

syms fs ls  r m0 q k x L s0
[email=f=@(x)sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2]f=@(x)sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2[/email]);
Q=quad(f,0,L)
结果为：??? Undefined function or method 'isfinite' for input arguments of type 'sym'.
Error in ==> quad at 75
if ~isfinite(y(1))

或改为：
syms fs ls  r m0 q k x L s0
[email=f=@sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2]f=sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2[/email]);
Q=int(f,x,0,L)
结果为：Warning: Explicit integral could not be found.
> In sym.int at 58

Q =

int((1+((m0+q/k^2)*(k*sinh(k*x)-k*tanh(1/2*k*L)*cosh(k*x))-q*x+1/2*q*L)^2/s0^2)^(1/2),x = 0 .. L)

请教高人，这是为什么啊？怎么不能求解呢？
谢谢！

[ 本帖最后由 simonwoo 于 2009-6-10 16:39 编辑 ]

feifeifool

结果说明没有解析解，可以采用数值积分的办法求解。