One of the most important uses of seismology science is subsurface imaging and seismic tomography that requires simulation. Simulation is accomplished by solving the wave equation. Due to the high processing cost, instead of solving the wave equation, using the WKB approximation and applying it to the wave equation, we obtain two simpler Eikonal and transport equations. In seismic exploration, the Eikonal solution is usually solved for high frequencies. This equation is a nonlinear equation for solving various methods that are divided into three categories of ray based, grid based and graph based methods.Fast marching method is a grid based method. In this method after the networking, the next receipt times are obtained based on the previous receipts. The fast marching method has a higher computational speed and more precision, and the method is more stable. In this method, taking into account a procedure (pseudo cone) at the top of which there is a source, and by obtaining time, at any given moment, one part of this procedure can be formed, thus, by obtaining the time of receipts in different moments, the shape of the cone is formed, and it is possible to obtain the motion of the wave front by obtaining this time of receipts in all points of the network.