clear all
n=17;
y = zeros(n,1);
y(1)=3.0001;  % y(ii) = 3^ii*x(ii)
y(2)=3.0001;
for ii = 3:n
     y(ii) = 13*y(ii-1) - 12*y(ii-2);
     ystar_ii = 3.0001;  % the true solution
     RelErr = abs(ystar_ii-y(ii))/ystar_ii;
     [ii  y(ii) ystar_ii  RelErr ]
%    fprintf('x(%2.0f)=%20.8d, xast(%2.0f)=%20.8d,’, ...
%    ’RelErr(%2.0f)=%14.4d\n’, ii,x(ii),ii,xn,ii,RelErr);
end