A fractional‐order Jacobi Tau method for a class of time‐fractional PDEs with variable coefficients