example1.log

./example1.out 

 "mfb.txt" No.00 -- focused beam --
medium density                       [kg/m^3] :            1000
speed of sound in the medium            [m/s] :            1500
frequency                                [Hz] :          100000
sound pressure amplitude                 [Pa] :     10+      0I
radius of focused transducer              [m] :           0.045
curvature radius of the transducer        [m] :           0.072
x-component of translation vector         [m] :               0
y-component of translation vector         [m] :               0
z-component of translation vector         [m] :               0
rotation parameter theta                [rad] :               0
rotation parameter psi                  [rad] :               0
sampling number for Gauss-Legendre quadrature :              64
-- additional information --
wavelength                                [m] :           0.015
energy passing through a plane orthogonal to the acoustic axis [W]
                                              : 1.4965495613794e-06

sound pressure and particle velocity at r=( 0.0005,-0.003, 0.001)
p  = 5.86589876401821e-01 -6.85792615191279e-01 I
v_x= 1.18761996799894e-08 +6.06129630464936e-09 I
v_y=-7.12571980799364e-08 -3.63677778278960e-08 I
v_z= 3.62391275553684e-07 -3.99114592538701e-07 I
Real value at r=( 0.0005,-0.003, 0.001), t=0.01 
p  = 5.86589876402261e-01
v_x= 1.18761996799855e-08
v_y=-7.12571980799130e-08
v_z= 3.62391275553941e-07

Directional derivative of sound pressure at r=( 0.0005,-0.003, 0.001)
n=( 0.107084, 0.275436, 0.955336),|n|=1
dpdn= 2.45456827863654e+02 +2.05994566683486e+02 I
Directional derivative using central difference. h=1.5e-09
dpdn= 2.45456827419055e+02 +2.05994566330562e+02 I