function [diff] = spacing(theta_max)
global hglob
global theta_min

diff = quad8('inteval', theta_min, theta_max, 1.e-8) - hglob;