Numerical solution of multiterm variable‐order fractional differential equations via shifted Legendre polynomials