|
matprod[a is array, b is array] :=
|
|
{
|
|
c = makeArray[[length[a], length[b@0]], 0]
|
|
|
|
a_row = length[a]-1
|
|
a_col = length[a@0]-1
|
|
b_col = length[b]-1
|
|
|
|
for row = 0 to a_row
|
|
for col = 0 to b_col
|
|
for inc = 0 to a_col
|
|
c@row@col = c@row@col + (a@row@inc * b@inc@col)
|
|
|
|
return c
|
|
}
|