Vieta–Fibonacci operational matrices for spectral solutions of variable-order fractional integro-differential equations