aa=load('Wetzell.SLR.enu');          %lecture SLR
bb=load('Wetzell.GPS.enu');          %lecture GPS
cc=load('Wetzell.VLBI.enu');         %lecture VLBI
t1=aa(:,1)';                         % Le temps est dans la première colonne
L=aa(:,2)';                          % la coordonnées E dans la seconde
t2=bb(:,1)';                         % On va devoir se ramener à un même indice de
P=bb(:,2)';                          % temps et ramener les séries sur les mêmes 
t3=cc(:,1)';                         % temps pour comparer. Toutefois, seuls les 
R=cc(:,2)';                          % temps de début et de fin sont différents.
mi=max([min(t1),min(t2),min(t3)]);   % On cherche les temps du début le plus tard des 3
ma=min([max(t1),max(t2),max(t3)]);   % et le temps de fin le plus tot des 3
i1=find(t1>mi&t1<ma);                % On cherche dans chaque séries quels sont les
i2=find(t2>mi&t2<ma);                % points dans cet intervalle.
i3=find(t3>mi&t3<ma);
t1=t1(i1);                           % Et on limite chaque série à cet intervalle.
t2=t2(i2);
t3=t3(i3);
L=L(i1);
P=P(i2);
R=R(i3);
A=[[1 1 0];[1 0 1];[0 1 1]];        % Ca, c'est la matrice du système d'équations
b(1)=var(L-P);                      % Ca, c'est le membre de droite.
b(2)=var(L-R);
b(3)=var(P-R);
no=A\b'                             % Et on résoud le système.