function de = determ(a, b, c, d)
% DETERM calculate determinant of 2x2 matrix [a b; c d]
   de = a*d - b*c;
end