The Method of Lines Combined with Chebyshev Spectral Method respect to weighted residual (Collocation Points) for Space-Time fractional diffusion equation is considered, the direct way will be used for approximating Time fractional and the expiation of shifted first kind of Chebyshev polynomial will be used to approximate unknown functions, the structure of the systems and the matrices will be fund, the algorithm steps is illustrated, The tables and figures of the results of the implementation by using this method at different values of fractional order will be shown, with the helping of programs of matlab.