function [s,p] = addmult(x,y)
% addmult  Compute sum and product of two matrices
s = x+y;
p = x*y;
