m=[.0698; 0.0743; 0.1194]; T=[0.78; 0.80; 1.00]; dT=0.01; dm=0.0001;

y=m; dy=dm;
x=T.^2; dx=2*T*dT; sigma=sqrt(dx.^2+dy^2);

A=sigma.^(-1);
A(:,2)=x.*sigma.^(-1);
b=y.*sigma.^(-1);

[Q,R]=QRdec(A); c=QRback(Q,R,b);

err=(R'*R)^(-1);
k=c(2)*(2*pi)^2
dk=sqrt(err(2,2))*(2*pi)^2
