A branch and prune algorithm for the approximation of non-linear AE-solution sets