In this paper, in a new approach, we study the stability of the dynamical system (DS) of R1+ϵ in terms of two significant cosmological parameters, deceleration parameter and jerk parameter {q,j}. Other cosmographic parameters such as (l,s,…) have been obtained in terms of these two parameters. We have obtained critical points (qe,je), the best fitted current values of cosmographic parameters (q0,j0,l0,s0), best value for model parameter ϵ and best trajectory of dynamics of system in phase space by simultaneously solving the DS and best fitting the parameter by the SNIa data. By defining modified redshift in anisotropic cosmological model as 1+z(t,ˆp)=a(t0)a(t)(1−A(ˆn⋅ˆp)) (where A is a magnitude of anisotropy, ˆn is direction of privileged axis and ˆp is the direction of each SNe Ia sample to galactic coordinates), the luminosity distance has been obtained in terms of modified redshift using cosmography method. Using union 2 data, we have found the direction of privileged axis in the galactic coordinate. The results show that the magnitude of anisotropy is about |A|≃10−3 and the direction of privileged axis is (l,b)=(298−34+34,2−23+23). Also, our results are consistent with other studies in 1−σ confidence level.