%Comparing Maclaurin series to the function
% Order starts at zero
function retval = eval_series(x_in, order)
  retval = 1;

  for i = 1:(order)
    retval = retval + x_in^i;
  endfor
endfunction