A numerical technique based on the shifted Legendre polynomials for solving the time-fractional coupled KdV equations