pressure = kzk_transient(transducer_array, coordinate_grid, medium, time_struct,excitation_function,\
                         di_flag, fd_mode, dsig_fd, piston_pts, rho_max);

[time, pressure] = kzk_transient(transducer_array, coordinate_grid, medium,time_struct,\
                                 excitation_function, di_flag, fd_mode, dsig_fd, piston_pts, rho_max); 

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.
• time_struct A time samples struct created by set_time_samples.
• excitation_function An excitation function created by set_excitation_function.
• di_flag Diffraction effect is included if di_flag = 1.
• fd_mode Implicit backward finite difference (IBFD) method is used when fd_mode = 1 and Crank-Nicolson finite difference (CNFD) method is used when fd_mode = 0.
• dsig_fd Normalized axial step size.
• piston_pts The number of radial divisions within a piston radius.
• rho_max Normalized max radial units.

Output Parameters

• time A vector of points in time from tmin to tmax.
• pressure A 4-d array representing the transient pressure at each point in spacetime.