[求助]积分方程求解

zw123

我还有个问题请教<br>y2=(a1*(x2-x)*a+a2*(x-x1)*b)/(a1*(x2-x1)+a2*(x-x1));<br>y=sqrt(1+(diff(y2))^2);<br>现在需要对x求积分,积分限是x1~x2,<br>我使用j=int(y,x1,x2)<br>为什么不行,<br>应该怎么处理?<br>谢谢各位!<br>

[此贴子已经被aspen于2006-3-25 19:11:56编辑过]

yqchenlolo

试了一下可以的<BR>你具体的函数是怎么定义的？a1,a2,a,b都是常数么？<BR>检查一下是不是哪里不小心搞错了<BR>或者把你之前各个变量的具体定义贴出来看看

happy

系数太多了，是解出来，不过估计有几千项

freelive

对于^1/2的积分计算机很慢,而且产生虚数<BR>解决方法对 y 在(x1+x2)/2处8阶用taylor展开,在用int <BR>注意taylor只能对符号展开,要先将y定义为符号<BR>sym x;<BR>taylor( y(x),(x1+x2)/2,8);