m=[.0698; 0.0743; 0.1194]; T=[0.78; 0.80; 1.00]; dT=0.01;
x=m*(2*pi)^2;
y=T.^2; dy=2*T*dT;
A=dy.^(-1);
A(:,2)=x.*dy.^(-1);
b=y.*dy.^(-1);

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

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