c - 构建错误导致代码块退出

Question

我正在编写 C 代码来解决数值问题，我定义了一个 100x100 矩阵并用值填充它，A 是矩阵，b 是解:A[i][j]=1/(i+j+1) ，b[i]=第 i 行中所有值的总和以下是我的代码:

#include <stdio.h>
#include <stdlib.h>
#define n 100
int main()
{
    double A[n][n];
    double b[n];
    int i,j;
    for (i=0;i<n;i++)
        for (j=0;j<n;j++)
        {
            b[i]=0;
            A[i][j]=0;
        }
    for (i=0;i<n;i++)
        for (j=0;j<n;j++)
        {
            A[i][j]=1/(i+j+1);
            b[i]+=A[i][j];
        }

    int c=10;

    for (i=0;i<n;i++)
        for (j=0;j<n;j++)
        {
             printf("%lf\t",A[i][j]);
             if (c==j) printf("\n");
             c=c*10;
        }


    return 0;
}

每当我单击构建按钮时，代码块就会终止，显示以下内容:

跨平台 IDE 停止工作

然后程序关闭。有人能帮我解决这个问题吗？!

Answer 1

您遇到的问题是所有A都将设置为0，只需使用1。请使用double而不是int

#include <stdio.h>                                                             
#include <stdlib.h>                                                            
#define n 100                                                                  
int main()                                                                     
{                                                                              
    double A[n][n];                                                            
    double b[n];                                                               
    int i, j;                                                                  
    for (i = 0; i < n; i++)                                                    
        for (j = 0; j < n; j++)                                                
        {                                                                      
            b[i] = 0;                                                          
            A[i][j] = 0;                                                       
        }                                                                      
    for (i = 0 ; i < n; i++)                                                   
        for (j = 0; j < n; j++)                                                
        {                                                                      
            A[i][j] = 1. / (i + j + 1);                                        
            b[i] += A[i][j];                                                   
        }                                                                      

    int c = 10;                                                                

    for (i = 0; i < n; i++)                                                    
        for (j = 0; j < n; j++)                                                
        {                                                                      
            printf("%lf\t", A[i][j]);                                          
            if (c == j) printf("\n");                                          
            c *= 10;                                                           
        }                                                                      


    return 0;                                                                  
}