The formation of soliton in optical fiber governed by nonlocal nonlinear Schrödinger (NLS) equation with fourth order dispersion is studied. The model of nonlocal NLS equation with fourth order dispersion is solved using the new extended auxiliary method and the solutions are obtained. The solutions are in the form of hyperbolic and trigonometric functions which are based on Jacobi elliptic function m. Shape changing soliton in optical fiber for nonlocal fourth order dispersive NLS equation is discussed by suitably choosing the values of kerr and quintic nonlinearities and by varying fourth order dispersion term. The effect of fourth order dispersion on soliton in fiber for different conditions of kerr and quintic nonlinearity is also discussed. In addition, the phase portraits of the system have been investigated and the stability of wave in optical fiber for nonlocal NLS equation is discussed using fourth order Runge-Kutta algorithm. This paper addresses a significant gap in the current literature by examining the impact of fourth order dispersion on the nonlocal NLS equation in optical fiber.