The most important aspect in the design of a systolic arrays is the mapping of the algorithm to the processor array. However, not all algorithms can be systolized. Only highly regular algorithms with the structure of nested loops are suitable for systolic implementation. Therefore, one has to choose an appropriate mathematical method that can be systolized. In this paper we analyze the problem of choosing the most suitable mathematicam method for systolic implementation. We illustrate this problem on the example of matrix multiplication.