#import rat

a =

<
   <1/1,  1/1,  1/1,   1/1>,
   <2/1,  4/1,  8/1,  16/1>,
   <3/1,  9/1, 27/1,  81/1>,
   <4/1, 16/1, 64/1, 256/1>>

b =

<
   <  4/1, -3/1,  4/3,  -1/4>,
   <-13/3, 19/4, -7/3,  11/24>,
   <  3/2, -2/1,  7/6,  -1/4>,
   < -1/6,  1/4, -1/6,   1/24>>

mmult = *rK7lD *rlD sum:-0.+ product*p

#cast %qLL

test = mmult(a,b)