plot을 그려보긴했는데 맞는지 모르겠네요 저도 독학으로만 끼적대본 수준이라

clear

clc

n=1:2000;

y(2)=0.1;

y(1)=0.2; 


x(1)=0.4;

x(2)=0.3;

x(3)=0.5;


for i=3:max(n)+3;

    y(i)=1.6*y(i-1)-0.64*y(i-2)+0.04*x(i-2);

    x(i)=(y(i)-1.6*y(i-1)+0.64*y(i-2))/0.04;

end;


subplot(211)

plot(n,y(4:2003))


subplot(212)

plot(n,x(4:2003))