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)