function FePred = FePrediction(Fe,Hp,n)

disp('In FePrediction');

h = length(Fe)-n;

A = zeros(h,n);
x = Fe(end:-1:n+1)';
for j = 1:n
     A(:,j) = Fe(end-j:-1:n-j+1)';
end

alpha = A\x;

k = length(Fe);
Fe = [Fe zeros(1,Hp)];

for i = 1:Hp
     Fe(k+i) = Fe((k+i-1:-1:k+i-n))*alpha;
end

FePred = Fe(k+1:end);



end

