function [out1, out2] = fibonacciFcn(in1, in2)
% documentation is here!!!
nArgIn = nargin;
if nArgIn == 0
   limit = 1e3;
   f = [0, 1];
elseif nArgIn == 1
   limit = in1;
   f = [0, 1];
elseif nArgIn == 2
   limit = in1;
   f = in2;
end
n = 1; % index for series generation
while f(n+1) < limit
   f(n+2) = f(n) + f(n+1);
   n = n + 1;
end
f(end) = [];
out1 = f;
if nargout == 2
   out2 = plot(f);
end
end