n = 10;
historySize = 500;
outputSize = 15;
x=zeros(1,historySize);
for i=1 : historySize
    x(i) = mod(i,2);
end
while(true)
    entry = input('Enter 0 or 1: ');
    x = [x(2:historySize),entry];
    xNew = x;
    output = zeros(1,outputSize);
    for j = 1:outputSize
        A = lpc(xNew,n);
        result = 0;
        for i=1:n
            result = result - A(i+1)*xNew(historySize-i+1);
        end
        if (result > 0.5)
            output(j) = 1;
        else
            output(j) = 0;
        end
        xNew = [xNew(2:historySize),output(j)];
    end
    output
end