High intensity focused ultrasound simulator

Simulates high intensity focused ultrasound beams and heating effects in layered media

Contains two modules. The first solves the KZK equation for the pressure distribution of an axisymmetric focused ultrasound transducer, produces salient plots, and dumps the computed data to the Matlab workspace. The second module uses the pressure distribution to compute heating and thermal dose according to a user-defined treatment protocol. Again plots are produced and data is available for postprocessing.