Integral transform method (Fourier or Laplace transform, etc) is more often effective to do the theoretical analysis for the stochastic processes. However, for the time-space coupled cases, e.g., Levy walk or nonlinear cases, integral transform method may fail to be so strong or even do not work again. Here we provide Hermite polynomial expansion approach, being complementary to integral transform method. Some statistical observables of general Levy walks are calculated by the Hermite polynomial expansion approach, and the comparisons are made when both the integral transform method and the newly introduced approach work well.