An approach to design fault-tolerant hexagonal systolic array (SA) for multiplication of rectangular matrices is described. The approach comprises three steps. First, redundancies are introduced at the computational level by deriving three equivalent algorithms but with disjoint index spaces. Second, we perform the accommodation of index spaces to the projection direction to obtain hexagonal SA with optimal number of processing elements (PE) for a given problem size. Finally, we perform mapping of the accommodated index spaces into fault-tolerant systolic array using valid transformation matrix. As a result, we obtain SA with optimal number of PEs which performs fault-tolerant matrix multiplication. In the case of square matrices of order N x N, this array comprises N 2 + 2N PEs with active computation time tc -----5N -4 time units. Fault tolerance is achieved through triplicated computation of the same problem instance and majority voting. We have proposed two hardware solutions for the voting process: one when voting is performed at the end of the computation, i.e., at the output of the SA, and the other where voting is performed after each computational step. With the proposed method, any single transient or permanent fault can be detected and corrected. Experimental results show that with the proposed schemes a lot of multiple error patterns can be tolerated, also.