arr5[i, 2] = -arr1[i] / m * f; }

//3.b.计算误差方程式Vy中dZs的系数 for (int i = 1; i < 2 * n; i += 2) {

arr5[i, 2] = -arr1[i] / m * f; }

//4.a.计算误差方程式Vx中dΦ的系数 for (int i = 0; i < 2 * n; i += 2) {

arr5[i, 3] = -f * (1 + arr1[i] * arr1[i] / f * f); }

//4.a.计算误差方程式Vy中dΦ的系数 for (int i = 1; i < 2 * n; i += 2) {

arr5[i, 3] = -arr1[i - 1] * arr1[i] / f; }

//5.a.计算误差方程式Vx中dΨ的系数 for (int i = 0; i < 2 * n; i += 2) {

arr5[i, 4] = -arr1[i] * arr1[i + 1] / f; }

//5.b.计算误差方程式Vy中dΨ的系数 for (int i = 1; i < 2 * n; i += 2) {

arr5[i, 4] = -f * (1 + arr1[i] * arr1[i] / f * f); }

//6.a.计算误差方程式Vx中dk的系数 for (int i = 0; i < 2 * n; i += 2) {

arr5[i, 5] = arr1[i + 1]; }

//6.b.计算误差方程式Vy中dk的系数 for (int i = 1; i < 2 * n; i += 2) {

arr5[i, 5] = -arr1[i - 1]; }

//定义外方位元素组成的数组

double[] arr6 = new double[6];//--------------------外方位元素改正数矩阵（X） //定义常数项元素组成的数组

double[] arr7 = new double[2 * n];//-----------------常数矩阵（L） //计算lx的值

for (int i = 0; i < 2 * n; i += 2)