alpha = 0.7
p(1) = 2
p(2) = 1
p(3) = .5
for i = 1:100
   p(i+4) = sin(p(i+3) + p(i+2) + p(i+1) - p(i)) + alpha*p(i)*(1-p(i))
end

plot(p,'o')