1. f:= ((1+8x)/(7+x))
Orb:=proc(f,x,x0,K1,K2) local x1,i,L:
x1:=1:
for i from 1 to 1000 do
 x1:=subs(x=x1,f):   
L := [op(L), x1]; 
od:
end:
L:=[x1]:




2.  f := k*x(n - 1)*(1 - x(n - 1))
K := 3.5:  
Orb:=proc(f,x,1,2,2.5, 3.5, 1000) local x1,i,L:  
x1:=1:  
for i from 1 to 1000 do   
x1 := subs(x = x1, f); 
L := [op(L), x1]; 
end 


3. xn= (x(n-1)+8x(n-2))/((x(n-1)+9x(n-2)), x0= 0.5, x1=0.7


f:= (x(n-1)+8x(n-2))/((x(n-1)+9x(n-2)):  
Orb:=proc(f,x,0.5,0.7,1,1000) local x1,i,L: 
 x1:=1:  
for i from 1 to K1 do   
x1:=subs(x=x1,f):    
x2:=subs(x=x2,f)   
od:  
L := [x1, x2]