% chaostest

clear
% *******************************************************
% *******************************************************
% *******************************************************
% ***  
% ***  Chaostest 
% ***  
% *******************************************************
% *******************************************************
% *******************************************************

load orbitdata
rho=intval(2e-7);
yy= [uuper1; uuhom1; uuhet12; uuper2; uuhom2; uuhet21];
hh=[hhper1; hhhom1; hhhet12; hhper2; hhhom2; hhhet21];
NN=length(hh);
index=[2:NN 1];

chaostest

% available data

% delta0
% delta1


% radius=rr
% rdurchrho=sup(rr/rho)
% maxeps=max(epsn)
% maxtilded_n = sup(maxeps/sqrt(1-delta3))
% faktor=Kinv*maxeps/Kinvd
% delta0strich=sup(Kinv*maxeps*sqrt((1+delta3)/(1-delta3))/(1-delta2))
% rrstrich=sup(delta0strich/(1-delta1))
% Linfnorm=Kinv*sqrt((1+delta3)/(1-delta3))/(1-delta2)
% sup(Linfnorm*maxeps)
