Vibration of a non-linear system with multiple degrees of freedom under random parametric excitations was evaluated by probabilistic method. The non-linear characteristic terms of system structure were quasi-linearized and excitation terms were remained as they were. An analytical method where the expectation values of square mean of error was minimized was used. The numerical results were compared with those obtained by Monte Carlo simulation. A linear congruential generator and Box-Muller method were used in Monte Carlo simulation. The comparison showed the results by probabilistic method agreed well with those by Monte Carlo simulation.