G exp(-alpha*(x^2+y^2)) 高斯函数
u_i c0*beta*x*G "x 向速度初始条件"
v_i c0*beta*y*G "y 向速度初始条件"
p_i rho0*c0^2*G "压力初始条件"
r sqrt((x-Ma*t)^2+y^2) "距离变量（动坐标系）"
p_a 1/(2*alpha)*integrate((cos(L*t)-beta/(2*alpha)*L*sin(L*t))*L*besselj(0,L*r)*exp(-L^2/(4*alpha)),L,0,10,1e-3) "解析压力表达式"