测量平差编程(10)

}

l=k*n+k;

a[k][k]=1.0/a[k][k]; for(j=0;j<=n-1;j++)

if(j!=k) {

u=k*n+j;a[k][j]=a[k][j]*a[k][k]; }

if(i!=k) if(j!=k) { } if(i!=k) { }

u=i*n+k;a[i][k]=-a[i][k]*a[k][k]; u=i*n+j;

a[i][j]=a[i][j]-a[i][k]*a[k][j];

for(i=0;i<=n-1;i++) for(j=0;j<=n-1;j++)

for(i=0;i<=n-1;i++)

}

for(k=n-1;k>=0;k--) {

if(js[k]!=k)

for(j=0;j<=n-1;j++) { }

for(i=0;i<=n-1;i++)

u=k*n+j;v=js[k]*n+j;

p=a[k][j];a[k][j]=a[js[k]][j];a[js[k]][j]=p;

if(is[k]!=k)