clear all
aa=load('../Données/data_air_PIB.txt');
d=aa(:,1);
x=aa(:,2);
lat=aa(:,3);
lon=aa(:,4);
dir=abs(mod(angle(lat-lat(1)+i*(lon-lon(1)))*180/pi-130,180));
[r,t]=corrsign(x,dir);
A=[0*dir'+1 ; dir'];
y=A'\x;
mod=A'*y;
res=x-mod;
C=inv(A*A');
s=sqrt(sum(res.^2)/(length(x)-2));
r=C(2,1)/sqrt(C(1,1)*C(2,2))