disp('水准网间接平差')
disp('已知高程')
Ha=5.015
Hb=6.016
disp('观测高差（m）')
L=[1.359;2.009;0.363;1.012;0.657;-0.357]
disp('系数矩阵')
B=[1,0;0,1;1,0;0,1;-1,1;-1,0]
l=[0;0;4;3;7;2]
disp('C:')
C=1*ones(1,6);
S=[1.1,1.7,2.3,2.7,2.4,4.0]
P=C./S
P=diag(P)
disp('参数的解')
x=inv(B'*P*B)*B'*P*L
disp('wucha')
V=B*x-l
L1=L+V/1000
disp('jingdu')
n=6;
t=2;
delta=sqrt(V'*P*V/(n-t))