clc;
clear all;
a=[1 2 3 4 5 6 7 8 9 10];
l=linspace(0,0,10);
l(1)=0.98;
for i=2:10
l(i)=1-2*l(i-1)*l(i-1);
end
[lsort,lindex]=sort(l);
for i=1:10
k=lindex(i);
b(i)=a(k);
end
b
[Vsort,Vindex]=sort(lindex);
for i=1:10
k=Vindex(i);
c(i)=b(k);
end
c