In this paper, the orthogonal polynomial approximation on triangle, proposed by Dubiner, is studied. Some approximation results are established in certain non-uniformly Jacobiweighted Sobolev space, which play important role in numerical analysis of spectral and triangle spectral element methods for differential equations on complex geometries. As an example, a model problem is considered.