A=[0 1; -1 -2]; B=[0;1]; C=[1 2]; D=0;

w=[0.1:0.05:1, 1.5:0.5:10, 15:5:100 150:50:1000];
m=[]; I=eye(size(A));
for i=w,
   m=[m; C*inv(sqrt(-1)*i*I-A)*B+D];
end
length(w)   % 顯示頻率點個數(shù)

w1=0.1:.1:1000; m1=[];
for i=w1,
   m1=[m1; C*inv(sqrt(-1)*i*I-A)*B+D];
end
length(m1)

w2=0.1:1:1000; m2=[];
for i=w2,
   m2=[m2; C*inv(sqrt(-1)*i*I-A)*B+D];
end
length(m2)

mysum = 0; i=1;
while (i<=100),
   mysum=mysum+i; i=i+1;
end