octave/fak.m

function f = fak(s)
  t = s - 1;
  if t <= 0 
    f = 1;
  else
    f = s * fak(t);
  endif 
endfunction