这样编出来的程序是不是不稳定？

部分代码如下：
float **Mr;
Mr=matrix_float(1000,7200); //用于存放线性拟合结果
float **Ma;
Ma=matrix_float(1000,7200);//用于存放年平均值结果
float x[21];
for(int i=0;i<21;i++)
{
x[i]=i;
}
float y[21];
float a[2];//用于存放线性拟合系数
for(int i=0;i<_YDIM;i++)
{
for(int j=0;j<_XDIM;j++)
{
for(int k=0;k<21;k++)
{
y[i]=M1[k][i*_XDIM+k];
}
sqt1(x,y,21,a);
Mr[i][j]=a[1];
}
}
运行之后总是出现 ××处有未经处理的异常: 0xC0000005: 读取位置 0x000000d0 时发生访问冲突
求大神解啊。。。