|       SUBROUTINE SMXM(A,N,B,M,C,L) | 
 | ** | 
 | **    C = A * B | 
 | **    A ET B MATRICES A(N,M)  B(M,L) ==> C(N,L) | 
 | ** | 
 | *>A      PREMIERE MATRICE | 
 | *>N      PREMIERE DIMENSION DE A ET DE C | 
 | *>B      DEUXIEME MATRICE | 
 | *>M      DEUXIEME DIMENSION DE A ET PERMIERE DE B | 
 | *<C      MATRICE PRODUIT DE A ET DE B | 
 | *>L      DEUXIEME DIMENSION DE B ET DE C | 
 | *A R. SANCHEZ ( EARLY WINTER 1987 ) | 
 | *V M.F. ROBEAU | 
 | *M AM BAUDRON - AVRIL 94 | 
 | *: ERREUR DANS L'APPEL A L'UTILITAIRE SGEMM | 
 | *: APPEL A L'UTILITAIRE SGEMM DE LA LIBRAIRIE BLAS SUR HP | 
 | *M AM BAUDRON - NOVEMBRE 1991 | 
 | *: ERREUR ( SOMME SUR LES TERMES PAS FAITE ) | 
 | *: APPEL A L'UTILITAIRE SGEMM DE LA LIBRAIRIE BLAS SUR RISC | 
 | *M AM BAUDRON - MAI 1993 | 
 | *: CHANGEMENT DES %IF LOCAL SUN MIPS SUITE A INTRODUCTION VERSION IBM | 
 |       DIMENSION A(N,M),B(M,L),C(N,L) | 
 |       DO 20 I=1,N | 
 |       DO 20 J=1,L | 
 |          R=0. | 
 |          DO 10 K=1,M | 
 |             R=R+A(I,K)*B(K,J) | 
 |    10    CONTINUE | 
 |          C(I,J)=R | 
 |    20 CONTINUE | 
 |       RETURN | 
 |       END |