解方程组（方程组中含有二重积分，并且积分的上下限中含有变量）

问题来源：http://www.ilovematlab.cn/thread-57006-1-3.html

问题描述：http://hi.csdn.net/attachment/200912/20/37536_1261296414KS8i.png

主要问题就是方程组中含有二重积分，并且积分的上下限中含有变量，所以我现在是无从下手，望高手指点。

Forcal代码：

(::p,q,m,C1,C2,C3,C4,k,g,Txing,T2)= { p=0.020,q=0.219,m=10369.6,C1=800,C2=2,C3=6,C4=8,k=3,g=4,Txing=8,T2=12 }; t_T2(u:a:p,q,m)= a=exp[-(p+q)*u],m*p*(p+q)^2*a/(p+q*a)^2; t2_T2(t::h_t_T2,T2,p,q,g)= if[t==T2,return(0)], XSLSF::fpqg[h_t_T2,t,T2,1e-6]/[g*(p+q)]; f(t2,T1,y1,y2:a,b:h_t_T2,h_t2_T2,p,q,m,C1,C2,C3,C4,k,g,Txing,T2)= //函数定义 { a=exp[-(p+q)*Txing], b=m*p*(p+q)^2*a/(p+q*a)^2, y1=-C2*XSLSF::fpqg[h_t2_T2,t2,T2,1e-6]+C4*b*[1-k*(p+q)] +C3*{k*(p+q)*XSLSF::fpqg[h_t_T2,T1,Txing,1e-6]+b*k*(p+q)*t2-k*(p+q)*Txing}, y2=C2*T1*T1/[2*g*(p+q)]-k*(p+q)*C3*(t2-T1)-C4*[1-k*(p+q)] }; main(:x,i,t2,T1:h_t_T2,h_t2_T2)= { h_t_T2=HFor("t_T2"),h_t2_T2=HFor("t2_T2"), x=new[rtoi(real_s),rtoi(4),rtoi(EndType),1.0,1.0], i=XSLSF::netn[HFor("f"),1e-6,0.1,0.1,x,100], x.XSLSF::getra[0,&t2,&T1], printff{"/r/n实际迭代次数={1,i}, t2={2,r}, T1={3,r}/r/n",i,t2,T1}, delete[x] };

结果：

实际迭代次数=4, t2=5.6891908532692748, T1=3.4017560817759569

当然，选择其他的初始参数时，还会有其他解。