source: mm.c procedure: mm format: rose loop: 0 known(ambn > 0) known(an > 0) known(bm > 0) shift_to(1,1,4) print