A=[1 1
1 0]
l1=(1+sqrt(5))/2
l2=(1-sqrt(5))/2
v1=[l1;1]
v2=[l2;1]
A*v1-l1*v1
A*v2-l2*v2
x=[1;0];
c=[v1 v2]\x;c1=c(1);
c1*l1^100*v1
Feig100=ans(1)

for i=1:100,
  x=A*x;
%disp(['F( ',num2str(i),') = ',num2str(x(1))]);
end
F100=x(1)
disp(['relative error (Feig100-F100)/Feig100 is  ',...
          num2str((Feig100-F100)/Feig100)]);
