A quadrature tau method for fractional differential equations with variable coefficients