/*************************************************************************
Program:marks(math3-7.c)
Programmer:zhuhao'feng
TESTED/COMPILED-Windows xp/Borland Turbo c 2.0
*************************************************************************/
#include<stdio.h>
#include<math.h>
void main()
{
    float a[3][3] = {5,1,0,1,5,1,0,1,5}, b[3] = {17,14,7},x[3] = {0},y[3] = {0};
    float temp;
    int i,j,k;
    for(i = 0;i < 3;i ++)
    {
    temp = 0;
    for(j = 0,k = i + 1;j < i;j ++)
        temp += a[j][i]*a[i][j];
    a[i][i] = a[i][i] - temp;
    a[i][k] = a[i][k]/a[i][i];
    }
    y[0] = b[0]/a[0][0];
    for(i = 1;i < 3;i ++)
    {
    temp = 0;
    for(j = 0;j < i;j ++)
        temp += y[j]*a[i][j];
    y[i] = (b[i] - temp)/a[i][i];
    }
    printf("the changed array is:\n");
    for(i = 0;i < 3;i ++)
    {
    for(j = 0;j < 3;j ++)
        printf("%10.5f",a[i][j]);
    printf("|%10.5f\n",y[i]);
    }
    x[2] = y[2];
    for(i = 1;i >= 0;i --)
    {
    temp = 0;
    for(j = 2;j > i;j --)
        temp += x[j]*a[i][j];
    x[i] = y[i] - temp;
    }
    printf("the array answer is x[3]=[");
    for(i = 0;i < 3;i ++)
    printf("%10.5f",x[i]);
    printf("]\n");
}
我调试的时候发现数组b[0]的值在红线处自动变化。
我在windows和Linux下都调试出线问题
请高手解答

不是在你红线外出界哟，
for(i = 0;i < 3;i ++)
    {
    temp = 0;
    for(j = 0,k = i + 1;j < i;j ++)
        temp += a[j][i]*a[i][j];
    a[i][i] = a[i][i] - temp;
    a[i][k] = a[i][k]/a[i][i]; /* 是在这条，如i值是2，a[i][k]就是越限了。 */
   }

谢谢楼您提出来的问题！
当时我也想过这个问题，
追赶法解线性方程组的时，下面代码
    for(i = 0;i < 3;i ++)
    {
    temp = 0;
    for(j = 0,k = i + 1;j < i;j ++)
        temp += a[j][i]*a[i][j];
    a[i][i] = a[i][i] - temp;
    a[i][k] = a[i][k]/a[i][i];
    }
当i = 2越界后。a[i][i]就变为了a[2][2],是数组的最后一个变量了。
而a[i][k]变为a[2][3]没法赋值。
从而达到了LU分解的目的。
个人观点！

在TC上可以通过，在别的编译器上可能会提示错误。你的想法应适何在TC编译上。

我调试的时候发现数组b[0]的值在红线处自动随机变化。
这是为什么？
请高手解答

谈不上高手，讨论讨论。
原因就在你的a[i][k]越界时赋了值，b数组首地址是接在a数组后，a[2][3]刚好是b的首地址，你对a[2][3]赋值，就是对b[0]赋了值。b[0]被改动，就不是原值了。