程序在这个邮箱里,行列式的计算没有问题,好像问题出在了把行列式变换上,就是slove(1,n)这个函数的运算上,那位大神可以看看问题到底出在哪…… 不胜感激!!!
#include <stdio.h>
#include <math.h>
void input(int n);
void output(int n);
double slove(int i , int n);
double A[10][10],B[1][10],C[1][10],D[1][10];
int main()
{
int n,j;
double W;
printf("\n请输入方程组的未知数个数:");
scanf("%d",&n);
input(n);
output(n);
W=slove(0,n);
for(j=1;j<n+1;j++)
{
printf("第%d次",j);
C[0][j]=slove(j,n);
}
for(j=1;j<n+1;j++)
{
D[0][j]=C[0][j]/W;
}
printf("\n解得:");
for(j=1;j<n+1;j++)
{
printf("\n X%d=%f",j,D[0][j]);
}
return 0;
}
void input(int n)
{
int row,col;
for(row=0;row<n;row++)
{
printf("\n请输入第%d个方程的系数(a1 a2 ... an)(a1X1+a1X2+...+anXn=b1):",row+1);
for(col=0;col<n;col++)
scanf("%lf",&A[row][col]);
}
for(row=0;row<n;row++)
{
printf("\n请输入第%d个方程的值(bi)(a1X1+a1X2+...+anXn=b1):",row+1);
scanf("%lf",&B[0][row]);
}
}
void output(int n)
{
int row,col;
for(row=0;row<n;row++)
{
for(col=0;col<n;col++)
{
printf(" %f X%d +",A[row][col],col+1);
}
printf("\b = %f\n",B[0][row]);
}
}
double slove(int i,int n)
{
double temp,result=1.0;
int switchtime=0,flag=0;
int row,nextrow,col,stemp;
if(i!=0)
{
for(row=0;row<n+1;row++)
temp=B[0][row];
printf("temp=%f\n",temp);
A[row][i-1]=temp;
}
for(row=0;row<n;row++)
{
for(col=0;col<n;col++) /*将行列式的列进行变换!!!!!!!!!!*/
printf("\t%f",A[row][col]);
printf("\n");
}
for(row=0;row<n-1;row++)
{
nextrow=row+1;
if(A[row][row]==0)/* 开始处理第一列,如果行列式第一行第一个数为零,要交换行 */
{ while(A[nextrow][row]==0)
{
nextrow++; /* 如果行列式第二行第一个数为零,行增加继续寻找非零数值的行 */
if(nextrow==n)/* 如果遍历完行列式行列式第一列元素都为零,退出while循环 */
{ flag=1;
break;
}
}
if(flag==1) /* 退出while循环后回到for(row=0;row<n-1;row++)行加1?*/
continue; /* 从array[row][row]==0知列也相应加1,开始处理第二列 */
switchtime++; /* 每交换一次行,行列式符号变化1次,统计变化次数 */
for(col=0;col<n;col++) /* 交换非零行到行列式顶部 */
{
stemp=A[row][col];
A[row][col]=A[nextrow][col];
A[nextrow][col]=stemp;
}
}
for(nextrow=row+1;nextrow<n;nextrow++)
{ /* 类似高斯消去法,消第一行下各行第一列数值到零*/
temp=A[nextrow][row]/A[row][row];
for(col=0;col<n;col++)
A[nextrow][col]+=-temp*A[row][col];/* 化行列式为上三角行列式形式 */
}
}
for(row=0;row<n;row++)
result*=A[row][row];
if(switchtime%2)
return -result;
else
return result;
}
#include <stdio.h>
#include <math.h>
void input(int n);
void output(int n);
double slove(int i , int n);
double A[10][10],B[1][10],C[1][10],D[1][10];
int main()
{
int n,j;
double W;
printf("\n请输入方程组的未知数个数:");
scanf("%d",&n);
input(n);
output(n);
W=slove(0,n);
for(j=1;j<n+1;j++)
{
printf("第%d次",j);
C[0][j]=slove(j,n);
}
for(j=1;j<n+1;j++)
{
D[0][j]=C[0][j]/W;
}
printf("\n解得:");
for(j=1;j<n+1;j++)
{
printf("\n X%d=%f",j,D[0][j]);
}
return 0;
}
void input(int n)
{
int row,col;
for(row=0;row<n;row++)
{
printf("\n请输入第%d个方程的系数(a1 a2 ... an)(a1X1+a1X2+...+anXn=b1):",row+1);
for(col=0;col<n;col++)
scanf("%lf",&A[row][col]);
}
for(row=0;row<n;row++)
{
printf("\n请输入第%d个方程的值(bi)(a1X1+a1X2+...+anXn=b1):",row+1);
scanf("%lf",&B[0][row]);
}
}
void output(int n)
{
int row,col;
for(row=0;row<n;row++)
{
for(col=0;col<n;col++)
{
printf(" %f X%d +",A[row][col],col+1);
}
printf("\b = %f\n",B[0][row]);
}
}
double slove(int i,int n)
{
double temp,result=1.0;
int switchtime=0,flag=0;
int row,nextrow,col,stemp;
if(i!=0)
{
for(row=0;row<n+1;row++)
temp=B[0][row];
printf("temp=%f\n",temp);
A[row][i-1]=temp;
}
for(row=0;row<n;row++)
{
for(col=0;col<n;col++) /*将行列式的列进行变换!!!!!!!!!!*/
printf("\t%f",A[row][col]);
printf("\n");
}
for(row=0;row<n-1;row++)
{
nextrow=row+1;
if(A[row][row]==0)/* 开始处理第一列,如果行列式第一行第一个数为零,要交换行 */
{ while(A[nextrow][row]==0)
{
nextrow++; /* 如果行列式第二行第一个数为零,行增加继续寻找非零数值的行 */
if(nextrow==n)/* 如果遍历完行列式行列式第一列元素都为零,退出while循环 */
{ flag=1;
break;
}
}
if(flag==1) /* 退出while循环后回到for(row=0;row<n-1;row++)行加1?*/
continue; /* 从array[row][row]==0知列也相应加1,开始处理第二列 */
switchtime++; /* 每交换一次行,行列式符号变化1次,统计变化次数 */
for(col=0;col<n;col++) /* 交换非零行到行列式顶部 */
{
stemp=A[row][col];
A[row][col]=A[nextrow][col];
A[nextrow][col]=stemp;
}
}
for(nextrow=row+1;nextrow<n;nextrow++)
{ /* 类似高斯消去法,消第一行下各行第一列数值到零*/
temp=A[nextrow][row]/A[row][row];
for(col=0;col<n;col++)
A[nextrow][col]+=-temp*A[row][col];/* 化行列式为上三角行列式形式 */
}
}
for(row=0;row<n;row++)
result*=A[row][row];
if(switchtime%2)
return -result;
else
return result;
}