积分请教

yqchenlolo

<P>syms x y T K; <br>f=12.5/pi/sqrt(1-400/T^2)*exp(-0.5*(x^2-40/T*x*y+y^2));<br>a= int(int(f,y,-inf,log(50/K)+0.07*T+0.4),x,-inf,log(50/K)+0.07*T+0.4);<br>这里a是个关于x,y的双重积分,T,K为变量<br>现在我用<br>b=subs(a,K,5000);<br>c=subs(b,T,40);<br>用cc=vpa(c)可以算出c的值为cc=1.2655795462299524135216594546451<br>也就是说T,K应该是a-cc=0的解<br>我用solve(subs(a-cc,K,5000),T)要解T,这时至少应该有一个解是cc<br>但是系统提示：<br>??? Error using ==&gt; solve<br>Error, (in solve) cannot solve expressions with, int(-2239906695008851/1125899906842624*Pi^(1/2)*2^(1/2)*exp(-1/2*x^2*(T-20)*(T+20)/T^2)*(erf(1/200*2^(1/2)*(-100*T*ln(2)-200*T*ln(5)+100*ln(5000)*T-7*T^2-40*T+2000*x)/T)-1)/((T^2-400)/T^2)^(1/2),x = -infinity .. -ln(100)+7/100*T+2/5), for, T</P>
<P>Error in ==&gt; D:\MATLAB6p5\toolbox\symbolic\@sym\solve.m<br>On line 49  ==&gt; [varargout{1:max(1,nargout)}] = solve(S{:});<br><br>提示的意思应该是这个表达式里面含有积分，无法解，<br>请问可有别的什么办法？</P>

[此贴子已经被作者于2006-4-6 13:48:20编辑过]

aspen

应该是你这个积分可以进行数值求解但是没有解析解