[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
不知道对不对，你参考一下