This paper presents a new numerical method to solve time fractional Fokker-Planck equation. The space dimension is discretized to the Gauss-Lobatto points, then we apply pseudo-spectral successive integration matrix for this dimension. This approach shows that with less number of points, we can approximate the solution with more accuracy. The numerical results of the examples are displayed.

We compute point sets on the triangle that have low Lebesgue constant, with sixfold symmetries and Gauss-Legendre-Lobatto distribution on the sides, up to interpolation degree 18. Such points have the best Lebesgue constants among the families of symmetric points used so far in the framework of triangular spectral elements.

We introduce a new and eecient Chebyshev-Legendre Galerkin method for elliptic problems. The new method is based on a Legendre-Galerkin formulation, but only the Chebyshev-Gauss-Lobatto points are used in the computation. Hence, it enjoys advantages of both the Legendre-Galerkin and Chebyshev-Galerkin methods.