We develop an explicit solution of problem describing collinear four-waves mixing in medium with cubic nonlinear response. This solution is carried out for set of Schrödinger equations using plane wave approximation for the case of phase matching for interacting waves. This solution allows to do full analysis of four-wave interaction modes in dependence of the problem parameters.We have shown, in particular, an existence of bistable mode for energy conversion from pump waves to signal wave under certain conditions. In general case, there are greater than 10 various modes of four-wave interaction. Knowledge about these modes is very important for spectroscopic experiment results understanding using four-waves mixing because its result depends on them in strong way. Analytical solutions and derived modes can explain complicated regime of four-wave interaction which may be appeared at high intensity of interacting waves.