怎样对比较复杂的符号表达式进行二重积分

yiyeguzhou206

在求解非线性方程时候，用
f(5)=int(int(ty1*r1*cos(ae+theta1),y1,-b1*sqrt(1-(x1/a1)^2),b1*sqrt(1-(x1/a1)^2)),x1,-a1,a1)+int(int(ty2*r2*cos(ae+theta2),y2,-b2*sqrt(1-(x2/a2)^2),b2*sqrt(1-(x2/a2)^2)),x2,-a2,a2)求解二重积分，积分区域是椭圆，里面的x1,y1,x2,y2是积分项;a1,a2,b1,b2是常数；ty1,r1,ae,theta1,theta2是由自变量组合的中间变量,形式也比较复杂,并且其中还含有x1,y1,x2,y2等积分项。
由于积分项太复杂，int(int())求解二重积分已经不能使用（自己已经试过），NIT工具箱里面的quad2dggen（）也只能求解数值积分的问题。哪位知道什么别的办法？还请不吝赐教！

sogooda

这么复杂的积分，恐怕解析解不可行。

yiyeguzhou206

确实不能解出来啊。但是还必须要这个积分积出来的自变量之间的关系，因为这个积分只是方程组中的一个方程。

[ 本帖最后由 ChaChing 于 2009-8-22 21:55 编辑 ]