The attached paper descibes the method. pages 65-66.

Thank you very much

you can use flipud in matlab, and you can get the opposite sequence.

for example:

clc

clear

cd('F:\PAPER CREAT\MK')

w=xlsread('.xls');

data=w;time=(1900:1990)';U=norminv(1-0.05/2,0,1);

N=size(data,1);n=2;nb=1;nB=2;ri=zeros(N,1);rbi=zeros(N,1);s=0;datac=flipud(data);

SFk=zeros(N,1);SBk=zeros(N,1);UFk=zeros(N,1);UBk=zeros(N,1);%data be limited in column

%in inverse order(RbI,UBK,SBK) & natural order(RI,UFK,SFK)

for n=2:N

ri(n)=sum(double(data(1:n-1)

Thank's for your help