Large N expansion of Wilson loops in the Gross–Witten–Wadia matrix model