  The paper presents the solutions for the zonal electrophoresis equations are
obtained by analytical and numerical methods. The method proposed by the
authors is used. This method allows to reduce the Cauchy problem for two
hyperbolic quasilinear PDE's to the Cauchy problem for ODE's. In some respect,
this method is analogous to the method of characteristics for two hyperbolic
equations. The method is effectively applicable in all cases when the explicit
expression for the Riemann-Green function of some linear second order PDE,
resulting from the use of the hodograph method for the original equations, is
known. One of the method advantages is the possibility of constructing a
multi-valued solutions. Compared with the previous authors paper, in which, in
particular, the shallow water equations are studied, here we investigate the
case when the Riemann-Green function can be represent as the sum of the terms
each of them is a product of two multipliers depended on different variables.
The numerical results for zonal electrophoresis equations are presented. For
computing the different initial data (periodic, wave packet, the Gaussian
distribution) are used.
