function assert(equation)
%% assert('equation') would return if equation is evaluated to
%% true, and would raise an error if not.

  if evalin('caller', equation, '0') == 0
    error(sprintf('assertion %s fails', equation));
  end

