A =[16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1];

B = [1     1
     0     1];
     
@(A,B)cell2mat(arrayfun(@(n)n*B,A,'un',0));

ans(A,B)