function [a  b] = M(n)     % M Memoized version of the Fibonacci number
  if n < 3                 % computation
     b=1; a=1;
  else
     [ b  c ] = M(n-1);
     a = b + c;
end
