a=[12345420212345432123484321];
subplot（121）;stem（a‘filled‘）;
xlabel（‘序列個數n‘）;
ylabel（‘幅值‘）;
title（‘原序列‘）;
axis（[0?25?0?25]）;
b=input（‘請輸入閾值b:\n‘）;
switch?b
????case?b==0.05
????????x=1;
????otherwise
????????x=2;
end
while（1）
????a=sort（a）;
????i=1:length（a）;
????n=length（a）;
????T=[0.988?0.889?0.780?0.698?0.637?0.683?0.635?0.597?0.679?0.642?0.651?0.641?0.616?0.595?0.577?0.561?0.547?0.535?0.524?0.514?0.505?0.497?0.489;0.341?0.765?0.642?0.560?0.507?0.554?0.512?0.477?0.576?0.546?0.521?0.546?0.525?0.507?0.490?0.475?0.462?0.450?0.440?0.430?0.421?0.413?0.406];
????if?n<=7
????????r101=（a（n）-a（n-1））/（a（n）-a（1））;
????????r102=（a（1）-a（2））/（a（1）-a（n））;
????????g=T（xn-2）;
????????if?r101>g
????????????fprintf（‘剔除的數是?%f\n‘a（n））;
????????????a（n）=[];
????????else?fprintf（‘沒有粗大誤差\n‘）;
????????end
????????if?r102>g
????????????fprintf（‘剔除的數是?%f\n‘a（1））;
????????????a（1）=[];
????????else?break;
????????end
????end
????if?（n<=8）&&（n<=10）
?