clear all;
%y(t+2) + 0.25y(t) = 5, t>0
%y(1) = 5, y(2) = 4,5
y(1) = 5;
y(2) = 4.5;
for m=3:30;
    y(m)=-0.25*y(m-2)+5;
end
t=[1:30];
plot(t,y,'*')
xlabel ('tempo')
ylabel ('y')
title ('Trajetória obtida no Matlab (*)')
pause
t=[1:30];
ayte(1)=5;
for i = 1:29
yte(i)=(0.5)^(i)*(cos((pi/2)*(i))+sin((pi/2)*(i))) + 4;
ayte(i+1)=yte(i);
end
plot(t,ayte,'o')
xlabel ('tempo')
ylabel ('y')
title ('Trajetória teórica (o)')
pause
plot(t,y,'*',t,ayte,'o')
xlabel ('tempo')
ylabel ('y')
title ('Trajetória obtida no Matlab (*) e teórica (o)')