Modeling Near Vertical Incidence Sounding (NVIS) High Frequency (HF) signal propagation in the ionosphere is important. Because, ionosondes which are special types of radars probing the ionosphere with certain HF frequencies (between 2 and 30 MHz), work mostly in NVIS mode (where elevation angle is between 89 and 90 degrees). In this work, we are going to propose a new method for NVIS wave propagation in the ionosphere by discretizing the NVIS wave propagation path into mediums in which the refractive index changes linearly, where we solve the ray propagation in each medium analytically using calculus of variations and use Snell's Law at medium changes. The main advantage of the proposed solution is the reduced computational complexity and time. This algorithm can be used to simulate and compare the behavior of vertical ionosondes together with other ray tracing algorithms. (C) 2018 Institute of Seismology, China Earthquake Administration, etc. Production and hosting by Elsevier B.V. on behalf of KeAi Communications Co., Ltd.