function f=Fib(a);
r(1)=0;
r(2)=1;
i=3;
while(i<=a+1)
  r(i)=r(i-1)+r(i-2);
  i=i+1;
end
f=r(a+1);
