pressure = kzk_cw(transducer_array, coordinate_grid, medium, f0, nharm,dsig_fd, piston_pts, rho_max,\
                  exflag); 

Arguments

• transducer_array A transducer_array.
• coordinate_grid A coordinate grid struct like the ones created by set_coordinate_grid.
• medium A medium struct like the ones created by set_medium.
• f0 Frequency of the transducer in Hz.
• nharm The number of harmonics to use in the calculation.
• dsig_fd Normalized axial step size.
• piston_pts The number of radial divisions within a piston radius.
• rho_max Normalized max radial units.
• exflag For circular transducers only. If 0, the untransformed coordinate system is used in the calculation. If 1, the coordinate system is transformed before the pressure calculation.

Output Parameters

• time A vector of points in time from tmin to tmax.
• pressure A 4-d array representing the complex pressure at each point in space for each harmonic such that pressure(x, y, z, n) is the pressure at point (x, y, z) for the nth harmonic.