Computing Néron–Tate heights of points on hyperelliptic Jacobians