> 
# Linearne rekurentne relacie 2.radu s konstantnymi koeficientami
> z[N]:=B1*z[N-1] +B0*z[N-2]+w[N];
> Y[1]:=Y1;Y[2]:=Y2;


> Digits:=10:
> m:=17:
> x:=array(1..m):
> A:=1:
> B:=1./3.:
> b1:=10./3.:b0:=-1.:
> a:=t->0.:
> x[1]:=A:x[2]:=B:
> for n from  3 by 1 to m do  x[n]:=b1*x[n-1] +b0*x[n-2]+a(n): od: print(x[m]);
# 
> with(plots):
> seqx:=1,x[1]:
>  for n from 2 to m do seqx:=seqx,n,x[n]:od:
> obr1:=plot([seqx],colour=black):
> display([obr1]);
#                    Teraz sa naucime ako sa riesia rekurentne relacie v MAPLE
> readlib(unassign):   # the purpose of these command is to eliminate the assignm
> ent 
> unassign('k');           
> b1:=10/3: b0:=-1:
> soly:=rsolve({y(k)=b1*y(k-1) + b0*y(k-2)+a(n),y(1)=A,y(2)=B },y);
# Zmente 10/3 na 10./3. a vykonajte znovu
>  evalf(subs(k=m,soly));         
> with(plots):
> seqy:=1,evalf(subs(k=1,soly)):
>  for i from 2 to m do
>  seqy:=seqy,i,evalf(subs(k=i,soly)):
> od:
> obr2:=plot([seqy],colour=red):
> display([obr2]);


> display([obr1,obr2]);


> seqx;
> seqy;
> with(plots):
> u:=evalf(subs(k=1,soly)):
> seqz:=1,(x[1]-u)/u:
>  for i from 2 to m do 
> u:=evalf(subs(k=i,soly)):
> seqz:=seqz,i,(x[i]-u)/u;
> od:
> obr5:=plot([seqz],colour=black):
> display([obr5]);


>