format long; % n = 40; n = 73; a = zeros(1,n); b = zeros(1,n); a(1) = 1; % a(2) = 1.01; % a(2) = 1.1; % a(2) = 2; a(2) = 1.5; for i = 3:n a(i) = a(i-1)+a(i-2); end b(n) = a(n); b(n-1) = a(n-1); for i = n-2:-1:1 b(i) = b(i+2)-b(i+1); end a0 = a(2)-a(1) b0 = b(2)-b(1)