using ApproxFun, Plots, ComplexPhasePortrait, ApproxFun,  SingularIntegralEquations,
        SpecialFunctions

using SingularIntegralEquations.HypergeometricFunctions
gr();

γ = Arc(0.,1., (0,π))

p = plot(γ; label="contour")
scatter!([0.],[0.]; label="singularity of a")
r = 0.5
for k = linspace(0.,π,10)
    plot!(Circle(exp(im*k),r); color=:green)
end
p