function x=bivarchaos(x)
% input vector of initial two values for x
for i=2:100
    x=[x 1-1.4*x(i)+0.3*x(i-1)];
end
plot(x,'.')