[问题] 帮忙求积分

[lijunbo1982]

各位高手，麻烦帮我看看这个求积分怎么做
syms ds;
p0=6*10^6;
v0=50*10^-3;
p2=16*10^6;
m=30;
n=1.4;
A=pi*0.125^2/4;
t=sqrt([2*ds*m*[p0-(3*ds*A/2*v0)^n]]/(p0*p2*A));
int(t,ds,0,0.53)

ans =

-131072000000000/5629734035232909*8042477193189870^(1/2)*9830400000000000000^(1/14)*629145600000000^(1/14)*(-1)^(13/14)/pi^6*(19663/944892805120000000000000000000*53000^(1/5)*209715200000000^(3/7)*29927310570716857910156250000000000000^(1/35)*pi^(67/10)*(-1)^(4/7)*hypergeom([-11/35, -1/10, 1/10, 3/10, 1/2, 7/10],[1/5, 2/5, 3/5, 24/35, 4/5],37778931862957161709568000000000000000000000000000000000000000000000000000/10572400258533/pi^7)*sec(1/5*pi)*csc(1/5*pi)*csc(1/10*pi)*csc(2/5*pi)-7/209715200000000*53^(4/5)*60514543933804184076288000000000000000000000000^(1/7)*209715200000000^(3/7)*(-1)^(4/7)*pi^(53/10)*hypergeom([-4/35, 1/10, 3/10, 1/2, 7/10, 9/10],[2/5, 3/5, 4/5, 31/35, 6/5],37778931862957161709568000000000000000000000000000000000000000000000000000/10572400258533/pi^7)*sec(3/10*pi)*sec(1/5*pi)*sec(2/5*pi)*csc(2/5*pi)+35/88*pi^5*csc(2/7*pi)*sec(27/70*pi)*sec(13/70*pi)*csc(3/35*pi)*csc(17/35*pi)*5^(1/2)*gamma(1/14)/gamma(4/7)+175/72*879609302220800000000000000000000000000000000000000000000000000^(1/7)*3^(29/35)*209715200000000^(3/7)*pi^(39/10)*(-1)^(4/7)*hypergeom([3/35, 3/10, 1/2, 7/10, 9/10, 11/10],[3/5, 4/5, 38/35, 6/5, 7/5],37778931862957161709568000000000000000000000000000000000000000000000000000/10572400258533/pi^7)*sec(1/10*pi)*csc(1/10*pi)*sec(1/5*pi)*sec(3/10*pi)+896000000000000/227529*737869762948382064640000000000000000000000000^(1/7)*81^(6/7)*pi^(5/2)*(-1)^(4/7)*hypergeom([2/7, 1/2, 7/10, 9/10, 11/10, 13/10],[4/5, 6/5, 9/7, 7/5, 8/5],37778931862957161709568000000000000000000000000000000000000000000000000000/10572400258533/pi^7)*csc(2/5*pi)*csc(1/10*pi)*sec(3/10*pi)*csc(3/10*pi)+183500800000000000000000000000/134138177*3022314549036572936765440000000000000000000000000000000000000000000^(1/7)*53^(3/5)*100^(2/5)*209715200000000^(3/7)*pi^(11/10)*(-1)^(4/7)*hypergeom([17/35, 7/10, 9/10, 11/10, 13/10, 3/2],[6/5, 7/5, 52/35, 8/5, 9/5],37778931862957161709568000000000000000000000000000000000000000000000000000/10572400258533/pi^7)*csc(1/5*pi)*csc(3/10*pi)*csc(1/10*pi)*csc(2/5*pi))
我想知道确切的数值，刚开始学matlab，水平有限，请各位高手帮忙！谢谢了

[tannengzhi0504]

clear;clc
syms ds;
p0=6*10^6;
v0=50*10^-3;
p2=16*10^6;
m=30;
n=1.4;
A=pi*0.125^2/4;
t=sqrt([2*ds*m*[p0-(3*ds*A/2*v0)^n]]/(p0*p2*A));
s=int(t,ds,0,0.53);
ss=vpa(s);
0.0044965893291238346267263438504854
不知道对不对，你参考一下

[justin_mao]

数值积分可以直接用quadl的呀···

p0=6*10^6;
v0=50*10^-3;
p2=16*10^6;
m=30;
n=1.4;
A=pi*0.125^2/4;
t=@(ds) sqrt([2.*ds.*m.*[p0-(3.*ds.*A./2.*v0).^n]]/(p0.*p2.*A));
quadl(t,0,0.53)

应该可以的