a=115.5434
b=152.1584
c=181.1240
A=43.7495*pi/200
B=62.9515*pi/200
C=93.2910*pi/200

T1=0.5*a*b*sin(C)
T2=0.5*a*c*sin(B)
T3=0.5*b*c*sin(A)

omega=200/pi

sigmaside2=power(0.01,2)
sigmaangle2=power(0.002,2)
K=diag([sigmaside2/8 sigmaside2/3 sigmaside2/5 sigmaangle2/2 sigmaangle2/4 sigmaangle2/3])
G=[T1/a T1/b 0 0 0 T1/(tan(C)*omega);T2/a 0 T2/c 0  T2/(tan(B)*omega) 0;
0 T3/b T3/c T3/(tan(A)*omega) 0 0]

VarT=G*K*transpose(G)
VarT(1,2)/sqrt(VarT(1,1)*VarT(2,2))

p=[1,1.4399,0.8553]
w=p/sum(p)
w*VarT*transpose(w)

ones=[1; 1; 1]
T=[T1; T2; T3]
inv(transpose(ones)*inv(VarT)*ones)*transpose(ones)*inv(VarT)*T
inv(transpose(ones)*inv(VarT)*ones)
