依据的拟声波公式如下： du/dt=1/rho*dp/dx , dw/dt=1/rho*dq/dz , dp/dt=rho*vpx^2*du/dx+rho*vp*vpn*dw/dz ,dq/dt=rho*vp*vpn*du/dx+rho*vp^2*dw/dz ,vpx^2=vp^2*(1+2*epsilon);vpn^2=vp^2*(1+2*delta);还是直接上程序：