Magic~笑颜常开 2020-08-10 16:58 采纳率: 0%
浏览 244
已采纳

使用Gauss-Jordan 消元法的计算

**学校的作业,让程序最后输出成如图这样,代码看了半天不知道怎么搞,自己的代码输出后是第三张图那样,一堆乱数,实在不明白,哪位高手来替小弟解下惑!
**

图片说明

图片说明

图片说明

图片说明

代码~~~~~~~

#include

#define N 100

void output(int n, int st, double a[100][N], double b[100])

{
int i, j;

printf("\n\nStep  %3d\n", st);
for (i = 1;i <= n;i++) {
    for (j = 1;j <= n;j++)
        printf("   %8.2f", a[i][j]);
    printf("   %8.2f\n", b[i]);
}

}
int main(void)
{
int i, j, k, n, q;
double p;
double a[N][N], b[N], x[N];

printf("namber of demension(n) =");        
scanf("%d", &n);
for (i = 1;i <= n;i++) {
    for (j = 1;j <= n;j++) {
        printf("a[%d][%d] =", i, j);
        scanf("%lf", &a[i][j]);
    }
    printf("b[%d] =", i);
    scanf("%lf", &b[i]);
}
printf("\n\nGauss eliminaition method\n\n");
printf("input data\n");
for (j = 1;j <= n;j++)
    printf("         x%d", j);
printf("       rhs\n");
for (i = 1;i <= n;i++) {
    for (j = 1;j <= n;j++)
        printf("   %8.2f", a[i][j]);
    printf("   %8.2f\n", b[i]);
}

for (i = 1;i <= n;i++) {
    p = a[i][i];
    for (j = i;j <= n;j++)
        a[i][j] = a[i][j] / p;
    b[i] = b[i] / p;
    for (k = 1;k <= n;k++) {
        if (k != i) {
            q = a[k][j];
            for (j = i;j <= n;j++)
                a[k][j] = a[k][j] - q * a[i][j];
            b[k] = b[k] - q * b[i];
        }
    }
    output(n, i, a, b);
}

printf("\n\n solution :\n\n"); 
for (i = 1; i <= n; i++)
    printf(" x[%d] =b[%d] %8.2f¥n", i, x[i], b[i]);
return 0;

}

  • 写回答

1条回答 默认 最新

  • 关注
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

悬赏问题

  • ¥15 file converter 转换格式失败 报错 Error marking filters as finished,如何解决?
  • ¥15 ubuntu系统下挂载磁盘上执行./提示权限不够
  • ¥15 Arcgis相交分析无法绘制一个或多个图形
  • ¥15 关于#r语言#的问题:差异分析前数据准备,报错Error in data[, sampleName1] : subscript out of bounds请问怎么解决呀以下是全部代码:
  • ¥15 seatunnel-web使用SQL组件时候后台报错,无法找到表格
  • ¥15 fpga自动售货机数码管(相关搜索:数字时钟)
  • ¥15 用前端向数据库插入数据,通过debug发现数据能走到后端,但是放行之后就会提示错误
  • ¥30 3天&7天&&15天&销量如何统计同一行
  • ¥30 帮我写一段可以读取LD2450数据并计算距离的Arduino代码
  • ¥15 飞机曲面部件如机翼,壁板等具体的孔位模型