二重积分求值

uaidenp

计算二重积分的时候，结果就是一个式子，不是积分值，这是为什么呢？求教
, X4 R3 a* U5 B6 P6 O$ L
9 R) G: s  p: J代码如下：
- a3 Z1 F7 u6 ?3 O8 E$ g
$ @& Y$ a$ i8 l) F- F" f7 fclear all
0 y& [$ ^/ ~0 @) V) ?1 Jclose all
( U0 k0 o, p# c+ Zclc
) o" l6 U. R: J5 Z+ z: G* Msyms x y
' E1 ~9 V) J% P, Q6 q( \) Pomega=1;
9 |8 v4 o5 l' c+ o) nmu=0.001;
alpha=0.13;
: J) m2 c7 k8 [% c9 ?beta=0.0031;
u=0.001;v=-0.02;
tau=2.2;tau1=1.5;D=0.11;
; q( F( @$ a; a9 c9 ]tau01=1.42;tau02=1.5;tau03=1.56;tau04=1.9;
; B/ N- y, w9 m/ L7 q/ K' ?omega2=1-u*cos(omega*tau)-v*omega*sin(omega*tau01);
. |4 J1 p$ J  jeta=1+(1/mu)*(-(u/omega)*sin(omega*tau)+v*cos(omega*tau01));
K=D/(omega2.*(1+omega2.*(tau1.^2)));
% i9 O7 T( t  W7 Sf1=(1./(omega2.*K)).*exp(-(mu./(1536.*K)).*(15.*beta.*(x.^2+y.^2./omega2).^4-32.*alpha.*(x.^2+y.^2./omega2).^3+96.*(x.^2+y.^2./omega2).^2-768.*eta.*(x.^2+y.^2./omega2)));
2 K& f, R) ], K% Q! eF1=int(int(f1,x,0,10),y,-10,10)
: E( l- V8 @9 v; R: B


% _& Q0 ^( o6 r: F% H结果：
( b% s5 \" k2 B  Z+ ^  l3 pF1=
int(int((1029484437153*exp((76368780973521777*(x^2 + (2251799813685248*y^2)/2297649917499219)^3)/922337203685477580800 - (17623564840043487*(x^2 + (2251799813685248*y^2)/2297649917499219)^2)/9223372036854775808 - (546330510041348097*(x^2 + (2251799813685248*y^2)/2297649917499219)^4)/590295810358705651712000 - (3488543873076974433998624506509*x^2)/81129638414606681695789005144064 - (1162847957692324811332874835503*y^2)/27593854166072889506733609713664))/34359738368, x, 0, 10), y, -10, 10)
0 X' l. ~7 h2 ]3 E

勇往直前11

来学习一下

uaidenp

lupkpu9 发表于 2020-7-2 10:16
double(F1)
第一句加上这个

可以了，谢谢哦
0 t3 `; q* N% g: E- s

lupkpu9

double(F1)
. O7 M) B8 x; ~; |$ N第一句加上这个

llbnmo

积分结果是个关于y的函数，如何转化成具体数值呢