#include <stdio.h>
#include <stdlib.h>
#include <string.h>
int main(int argc, char const *argv[])
{
    double a[621][621];
    //for (int i = 0; i < 621; i++)
    //    for (int j = 0; j < 621; j++)
    //        a[i][j] = 0;


    FILE *fpread = fopen("f:\\666.txt", "r");
    if (fpread == NULL)
    {
        printf("ERROR::FAIL TO READ!\n");
        exit(-1);
    }
    int i = 0;
    for (int i = 0; i < 621; i++)
        for (int j = 0; j <621; j++)
        {
            fscanf(fpread, "%le", &a[i][j]);
        }
     for(int i = 0; i < 621; i++)
     {
         for(int j = 0; j < 621; j++)
         {
             printf("%le ", a[i][j]);
         }

     }
    printf("%le\n", a[0][1]);
    system("pause");
    return 0;
}

栈溢出问题怎么解决？这才621*621的矩阵，以后还要用更大的

那就用堆呗
（double a[621][621] 大约3M，而VC默认只有1M，虽然可以调整，但没必要，还是malloc比较好）