Laboratory as well as Particle in cell (PIC) simulation experiments reveal the strong flow of energetic electrons co-moving with laser beam in laser plasma interaction. Equation governing the evolution of complex envelope in slowly varying envelope approximation is nonlinear parabolic equation. A Lagrangian for the problem is set up and assuming a trial Gaussian profile, we solve the reduced Lagrangian problem for beam width and curvature. Besides self-focusing and self-modulation of laser beam, we observe that stability properties of such plasma system are studied about equilibrium values using this variational approach. We obtained an eigen value equation, which is cubic in nature and investigated the criterion for stability using Hurwitz conditions for laser beam plasma system.