In this paper, we have obtained the numerical solutions of complex differential equations with variable coefficients by using the Legendre Polynomials and we have performed it on two test problems. Then, we applied with different technical of error analysis to the test problems. When we compared exact solutions and numerical solutions of tables and graphs, we realized that our method is reliable, practical, and functional.