In this paper, a numerical method is developed to compute an approximate solution of high-order linear complex differential equations in elliptic domains. By using collocation points and Bessel polynomials, this method ...