reccurHom.mws

Решение однородного рекуррентного уравнения

с постоянными коэффициентами

>    restart;

Дано уравнение и начальные члены

>    a:=-2:b:=3:x0:=-4:x1:=5:

>    eq:=x(n+2)=a*x(n+1)+b*x(n);

eq := x(n+2) = -2*x(n+1)+3*x(n)

Характеристическое уравнение

>    h:=z^2-a*z-b;

h := z^2+2*z-3

>    lambda:=solve(h,z);

lambda := 1, -3

Решение рекуррентного уравнения

>    X:=c1*(lambda[1])^n+c2*(lambda[2])^n;

X := c1+c2*(-3)^n

Находим константы из начальных условий

>    u1:=subs(n=0,X)=x0;

u1 := c1+c2 = -4

>    u2:=subs(n=1,X)=x1;

u2 := c1-3*c2 = 5

>    solve({u1,u2},{c1,c2});

{c2 = -9/4, c1 = -7/4}