Chaps11.mw

Maple Text Commands for Chapter 11 -  

LINEAR DISCRETE DYNAMICAL SYSTEMS 

 

Example 1: Solving recurrence relations x[n+1] = `+`(1, 3/100)*x[n], x[0] = 10000. 

> restart:b:=evalf(rsolve({x(n+1)=(1+(3/(100)))*x(n),x(0)=10000},x)):
n:=5:evalf(b,7);
 

11592.74 

Example 2(i): Second-order linear difference equation x[n+2] = 4*x[n+1]-4*x[n], x[0] = 1, x[1] = 2. 

> n:='n':
factor(rsolve({x(n+2)=x(n+1)+6*x(n),x(0)=1,x(1)=2},x));
 

1/5*(-2)^n+4/5*3^n 

Example 4: Leslie matrices. 

> with(LinearAlgebra):
 

> L:=Matrix([[0,3,1],[0.3,0,0],[0,0.5,0]]);
 

L := Matrix(%id = 147796120) 

> X0:=<<1000,2000,3000>>;
 

X0 := Vector[column](%id = 151079112) 

> X10:=L^(10).X0;
 

X10 := Vector[column](%id = 151378876) 

> Eigenvectors(L);
 

Vector[column](%id = 147067464), Matrix(%id = 145940868)
Vector[column](%id = 147067464), Matrix(%id = 145940868)
Vector[column](%id = 147067464), Matrix(%id = 145940868)
Vector[column](%id = 147067464), Matrix(%id = 145940868)
Vector[column](%id = 147067464), Matrix(%id = 145940868)
Vector[column](%id = 147067464), Matrix(%id = 145940868)
 

End of Chapter 11 Commands