Sadovsky equation example

% Distance from the epicenter (meters)
R = 2000;

% TNT equivalent of the blast power (kg)
q = 20 * 1000;

% Pressure change in the blast wave (atm)
dP = 1.06 * q^(1/3) / R + 4.3 * q^(2/3) / R^2 + 14 * q / R^3;