n=input(' Index of Fibonacci number = ');             % p23_fib.m
f=zeros(n,1);
f(1)=1; 
f(2)=1;
 for ii = 3:1:n ;
   f(ii)= f(ii-1)+f(ii-2);
 end

 fprintf('Fibonacci(%d) = %d \n',n,f(n));
