function [Rxx]=autom(x)
% [Rxx]=autom(x)

N=length(x);
Rxx=zeros(1,N);
for m=1: N+1
for n=1: N-m+1
Rxx(m)=Rxx(m)+x(n)*x(n+m-1);
end;
end;

x=[8.121132
8.580819
8.274361
6.2823853
5.209783
2.6048915
4.4436383
4.4436383
12.411542
5.209783
0.45968673
2.4516625
-2.9113493
-0.15322891
-0.91937345
-1.3790601
-0.15322891
-0.7661445
-1.3790601
-0.45968673
0.15322891
1.2258313
-2.7581203
3.5242648
-1.9919758
-0.30645782
-1.9919758
-1.0726024
-2.1452048
0.91937345
-0.45968673
2.1452048
-0.30645782
0.61291564
-1.2258313
0.30645782
-3.8307228
-3.371036
-0.30645782
1.532289
0.61291564
-0.61291564
-3.6774938
-0.45968673
-3.064578
0.45968673
-1.8387469
1.8387469
2.1452048
-1.8387469
-0.61291564
-0.15322891
-0.45968673
-6.742072
-2.7581203
-2.9113493
2.6048915
10.113108
1.3790601
0.30645782
-1.3790601
-6.2823853
-5.66947
-2.7581203
-1.532289
3.5242648
-1.0726024
-0.91937345
-0.30645782
-1.9919758
-3.5242648
1.2258313
1.8387469
-1.685518
2.9113493
0.7661445
1.0726024
3.217807
3.064578
1.532289
1.3790601
-0.7661445
1.0726024
3.371036
4.7500963
1.9919758
7.6614456
5.5162406
4.1371803
-0.15322891
0.45968673
0.15322891
1.3790601
];