We describe an algorithm for the numerical solution of an integral equation of the form $$ -\frac1{i}ıt_{-1}^1eft[(y-x)n|y-x|-h(x,y)\right]\frac{u(y)\,dy}{qrt{1-y^2}}=f(x),\quad-1