Observe um código m-file do software Matlab.

clear all

clc

format long

T=0.01;

x1=[0];

x2=[1];

y=[0];

u=[0];

tempo=[0];

for k=2:1:500

u(k)=1;

x1(k)=x1(k-1)+T*x2(k-1);

x2(k)=(-6*T)*x1(k-1)+(-5*T+1)*x2(k-1)+T*u(k-1);

y(k)=x1(k);

tempo(k)=(k-1)*T;

end

plot(tempo,y,'r')

hold on

plot(tempo,x1,'b')

plot(tempo,x2,'o')

O valor do terceiro elemento de y é: