An advanced numerical method of parallel shooting is proposed for solving steady state boundary problems in membrane electrochemistry. The advanced algorithm is described. The possibilities of the method are discussed. The possible applications to solving some boundary tasks of ordinary and singular perturbated differential equations in many-layers membrane systems are indicated. The transport task of two ions in three-layers systems at over-limiting current is solved.