Journal of Neutrosophic and Fuzzy Systems JNFS 2771-6449 2771-6430 10.54216/JNFS https://www.americaspg.com/journals/show/1614 2021 2021 Department of fundamental sciences, Higher Institute for Applied Sciences and Technology, Aleppo, Syria George A. Toma Department of mathematics, Faculty of Science, Aleppo University, Aleppo, Syria Fahed Farhood Department of mathematics, Faculty of Science, Aleppo University, Aleppo, Syria Taqi A. Alkhatib In this paper, an analytical method (Homotopy perturbation method HPM) is used for solving the initial value problem represented by a neutrosophic nonlinear Van der Pol oscillator equation (N-VDP) arising in applied dynamics using the thick function. We find the solutions of the (N-VDP) equation by HPM and then compare the numerical results with fourth order Runge-Kutta method (RK4). The results showed that the HPM lead to accurate and efficient results. Furthermore, these results of the HPM scheme and RK4 are implemented in Matlab. 2023 2023 14 20 10.54216/JNFS.060102 https://www.americaspg.com/articleinfo/24/show/1614