平方根法.c

数值计算的源代码

#include<stdio.h>
#include<math.h>
void main()
{
	
float a[11][12],t;
	int n,i,j,k;

    printf("LL^Tx=b,L^Tx=y\n");
	printf("input unknown numbers<max=10>:");
    scanf("%d",&n);
	printf("input all quotieties:\n");

	for(i=1;i<=n;i++)
		for(j=1;j<=n+1;j++)
			scanf("%f",&a[i][j]);

printf("The equation:\n");
for(i=1;i<=n;i++)
{
for(j=1;j<=n+1;j++)
	printf("%f ",a[i][j]);
printf("\n");
}

	for(i=1;i<=n;i++)
	{t=0;
		for(k=1;k<=i-1;k++)t+=a[i][k]*a[i][k];
		a[i][i]=pow((a[i][i]-t),0.5);

		for(j=i+1;j<=n;j++)
		{t=0;
		for(k=1;k<=i-1;k++)t+=a[j][k]*a[i][k];

		a[j][i]=(a[j][i]-t)/a[i][i];
		}
	}
	

for(i=1;i<=n;i++)
for(j=i+1;j<=n;j++)
	a[i][j]=a[j][i];


printf("The change:\n");
for(i=1;i<=n;i++)
{
for(j=1;j<=n+1;j++)
	printf("%f ",a[i][j]);
printf("\n");
}


for(i=1;i<=n;i++)
{t=0;
for(k=1;k<=i-1;k++)t+=a[i][k]*a[k][n+1];
a[i][n+1]=(a[i][n+1]-t)/a[i][i];
}

printf("y:\n");
for(i=1;i<=n;i++)
printf("%f ",a[i][n+1]);

for(i=n;i>=1;i--)
{t=0;
for(k=i+1;k<=n;k++)t+=a[k][i]*a[k][n+1];
a[i][n+1]=(a[i][n+1]-t)/a[i][i];
}
printf("\nx:\n");
for(i=1;i<=n;i++)
printf("%f ",a[i][n+1]);


}