funcao le_matriz (int l, int c): vetor(real, *, *)
    int i, j
    vetor(real, l, c) a
    para i em (0..l):
        para j em (0..c):
            leia a[i][j]
    retorne a

funcao mult_mat (vetor(real, *, *) a,
                 vetor(rael, *, *) b): vetor(rael, *, *)
    int i, j, k
    int l1 <- linhas(a)
    int c1 <- colunas(a)
    int l2 <- linhas(b)
    int c2 <- colunas(b)
    vetor(real, l1, c2) c
    para i em (0..l1-1):
        para j em (0..c2-1):
            para k em (0..c1-1):
                c[i,j] <- c[i,j] + a[i,k]*b[k,j]
    retorne c

funcao mostra_mat (vetor (real, *, *) a):
    int i, j
    int l <- linhas(a)
    int c <- colunas(a)
    para i em (0..l):
        para j em (0..c):
            mostre a[i,j]
        nova_linha


int m1, n1, m2, n2
vetor(real, *, *) a, b, c

leia m1, n1, m2, n2

se (n1 != m2):
    mostre "As matrizes não podem ser multiplicadas!"
    pare
    
a <- le_matriz(m1, n1)
b <- le_matriz(m2, n2)

c <- mult_mat (a, b)

mostra_mat (c)