Analytic study on the generalized ($$3+1$$3+1)-dimensional nonlinear Schrödinger equation with variable coefficients in the inhomogeneous optical fiber