G exp(-alpha*(x^2+y^2)) 高斯函数
p_i rho0*c0^2*G "压力初始条件"
u_i c0*beta*x*G "x 向速度初始条件"
v_i c0*beta*y*G "y 向速度初始条件"
dp_i "-c0^2*rho0*(d(u_i, x) + d(v_i, y))" "压力时间导数初始条件"
r sqrt(x^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) "解析压力表达式"