We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius' method and Moebius transformations in the vicinity of each of the singular points of the hypergeometric equation, which leads to a natural decomposition of the real axis into domains. In each domain, solutions to the hypergeometric equation are constructed via the wellconditioned ultraspherical spectral method. The solutions are matched at the domain boundaries to lead to a solution which is analytic on the whole compactified real line R∪∞, except for the singular points and cuts of the Riemann surface on which the solution is defined. The solution is further extended to the whole Riemann sphere by using the same approach for ellipses enclosing the singularities. The hypergeometric equation is solved on the ellipses with the boundary data from the real axis. This solution is continued as a harmonic function to the interior of the disk by solving the Laplace equation in polar coordinates with an optimal complexity Fourier-ultraspherical spectral method.

The inverse scattering approach for the defocusing Davey-Stewartson II equation is given by a system of D-bar equations. We present a numerical approach to semi-classical D-bar problems for real analytic rapidly decreasing potentials. We treat the D-bar problem as a complex linear second order integral equation which is solved with discrete Fourier transforms complemented by a regularization of the singular parts by explicit analytic computation. The resulting algebraic equation is solved either by fixed point iterations or GMRES. Several examples for small values of the semi-classical parameter in the system are discussed.

We present an efficient high-precision numerical approach for Davey–Stewartson (DS) II type equa- tions, treating initial data from the Schwartz class of smooth, rapidly decreasing functions. As with previous approaches, the presented code uses discrete Fourier transforms for the spatial dependence and Driscoll’s composite Runge–Kutta method for the time dependence. Since DS equations are non-local, nonlinear Schrödinger equations with a singular symbol for the non-locality, standard Fourier methods in practice only reach accuracy of the order of 10 −6 or less for typical examples. This was previously demonstrated for the defocusing integrable case by comparison with a numerical approach for DS II via inverse scattering. By applying a regularization to the singular symbol, originally developed for D-bar problems, the presented code is shown to reach machine precision. The code can treat integrable and non-integrable DS II equations. Moreover, it has the same numerical complexity as existing codes for DS II. Several examples for the integrable defocusing DS II equation are discussed as test cases. In an appendix by C. Kalla, a doubly periodic solution to the defocusing DS II equation is presented, providing a test for direct DS codes based on Fourier methods.