Abstract: Scoliosis deformity is characterized as a lateral deviation of the spine in three-dimensions accompanied by axial rotation of the vertebrae. In this paper by surveying on the mathematical model of spine movement, we simulate the growth of scoliosis disease of a patient by using finite element method for vertebrae T1-T6 of spine. For this purpose by using Newmark method, we first deal the numerical solution of the given mathematical model. We will also explain the growth process of this disease by using geometric drawing in Matlab environment. Moreover, numerical simulations based on the mentioned process are presented.