A numerical solution to the Gel'fand-Levitan-Marchenko equation (Q1126626)

A numerical method is developed to solve the inverse scattering problem associated with two Gel'fand-Levitan-Marchenko integral equations. This inverse problem is represented by two coupled partial differential equations \[ {\partial A_1(x,y)\over\partial y}+ {\partial A_1(x,y)\over\partial x}= q(x) A_2(x, y),\tag{1} \] \[ {\partial A_2(x,y)\over\partial y}- {\partial A_2(x, y)\over\partial x}= q^*(x) A_1(x,y),\tag{2} \] where \[ q^*(x)= -2A_2(x,x),\tag{3} \] \(0\leq x\), \(y\leq L\). An approximate solution of the system (1)--(2) is found in the form of the interpolated polynomials \[ A^N_1(\alpha, \beta)= \sum^N_{i=0} \sum^N_{j=0} \gamma_{ij}\varphi_i(\alpha)\varphi_j(\beta),\tag{4} \] \[ A^N_2(\alpha, \beta)= \sum^N_{i= 0} \sum^N_{j= 0} \zeta_{ij}\varphi_i(\alpha)\varphi_j(\beta),\tag{5} \] where \(\varphi_k(t)={1\over N(N+ 1)L_N(t_k)} {(t^2- 1)L_N'\over(t- t_k)}\), \(k= 0,1,\dots, N\), \(L_N(t)\), \(-1\leq t\leq 1\), denote the Legendre polynomial of order \(N\), with the Legendre-Gauss-Lobatto nodes, \(\alpha= 2x/L- 1\), \(\beta= 2y/L- 1\). From equation (3) the expression \[ q^N(\alpha)= -2 \sum^N_{i=1} \sum^N_{j=0} \zeta^*_{ij}\varphi_i(\alpha) \varphi_j(\beta)\tag{6} \] occurs. Substituting (4), (5) and (6) in equations (1) and (2) the authors obtain a system of algebraic equations for the determination of the coefficients \(\zeta_{ij}\) and \(\gamma_{ij}\). Numerical examples are given.