y=ziegelwanger2014_offaxis(p,x)
| p | off-axis time-of-arrival model parameters [SI-units] | 
| x | HRTF direction (azimuth,elevation) [rad] | 
toa=ziegelwanger2014_offaxis(p,x) calculates time-of-arrivals for given model parameters (p) and directions (x) with an off-axis time-of-arrival model.
H. Ziegelwanger and P. Majdak. Modeling the direction-continuous time-of-arrival in head-related transfer functions. J. Acoust. Soc. Am., 135:1278--1293, 2014.