#include?
#include?
#define?N?10?//矩陣大小范圍
/*
*?使用已經(jīng)求出的x，向前計(jì)算x（供getx（）調(diào)用）
*?float?a[][]	矩陣U
*?float?x[]		方程組解
*?int	i		解的序號(hào)（數(shù)組X元素序號(hào)）
*?int	n		矩陣大小
*?return		公式中需要的和
*/
float?getmx（float?a[N][N]?float?x[N]?int?i?int?n）
{
float?mx?=?0;
int?r;
for（r=i+1;?r{
mx?+=?a[i][r]?*?x[r];
}
return?mx;
}
/*
*?使用已經(jīng)求出的y，向前計(jì)算y（供gety（）調(diào)用）
*?float?a[][]	矩陣L
*?float?y[]		數(shù)組Y
*?int	i		數(shù)組Y元素序號(hào)
*?int	n		矩陣大小
*?return		公式中需要的和
*/
float?getmy（float?a[N][N]?float?y[N]?int?i?int?n）
{
float?my?=?0;
int?r;
for（r=0;?r{
if（i?!=?r）?my?+=?a[i][r]?*?y[r];
}
return?my;
}
/*
*?解方程組，計(jì)算某x
*?float?a[][]	矩陣U
*?float?x[]		方程組解
*?int	i		解的序號(hào)
*?int	n		矩陣大小
*?return		方程組的第i個(gè)解（數(shù)組X的第i個(gè)元素值）
*/
float?getx（float?a[N][N]?float?b[N]?float?x[N]?int?i?int?n）
{
float?result;
if（i==n-1）?//計(jì)算最后一個(gè)x的值
result?=?（float）（b[i]/a[n-1][n-1]）;
else?//計(jì)算其他x值（對(duì)于公式中的求和部分，需要調(diào)用getmx（）函數(shù)）
result?=?（float）（（b[i]-getmx（axin））/a[i][i]）;
return?result;
}
/*
*?解數(shù)組Y，計(jì)算其中一元素值
*?float?a[][]	矩陣L
*?float?y[]		數(shù)組Y
*?int	i		數(shù)組Y元素序號(hào)
*?int	n		矩陣大小
*?return		數(shù)組Y的第i個(gè)元素值
*/
float?gety（float?a[N][N]?float?b[N]?float?y[N]?int?i?int?n）
{
float?result;
if（i==0）?//計(jì)算第一個(gè)y的值
result?=?float（b[i]/a[i][i]）;
else?//計(jì)算其他y值（對(duì)于公式中的求和部分，需要調(diào)用getmy（）函數(shù)）
result?=?float（（b[i]-getmy（ayin））/a[i][i]）;
return?result;
}
void?main（）
{
float?l[N][N]={0};?//定義L矩陣
float?u[N][N]={0};?//定義U矩陣
float?y[N]={0};?//定義數(shù)組Y
float?x[N]={0};?//定義數(shù)組X
float?a[N][N]={{211}{132}{122}};?//定義系數(shù)矩陣
float?b[N]={465};?//定義右端項(xiàng)
float?sum=0;
int?ijk;
int?n=3;
int?flag=1;