平方根法解方程组
需积分: 28 22 浏览量
2018-03-19
09:35:40
上传
评论 1
收藏 80KB DOCX 举报
[
6 2 1 − 1
2 4 1 0
1 1 4 − 1
− 1 0 −1 3
]
[
x
1
x
2
x
3
x
4
]
¿
[
−1
2
1
3
]
求解以上方程组
#include <math.h>
#include <stdio.h>
#dene N 4
void Cholesky(double a[][N],double b[],double L[][N],double d[],double x[],double
y[],int n)
{
int i,j,k;
double s;
for(i=0;i<n;i++)
L[i][i]=1;
d[0]=a[0][0];
for(i=1;i<n;i++)
{
for(j=0;j<i;j++)
{
s=0;
for(k=0;k<j;k++)
s=s+L[i][k]*d[k]*L[j][k];
L[i][j]=(a[i][j]-s)/d[j];
}
s=0;
for(k=0;k<i;k++)
s=s+L[i][k]*L[i][k]*d[k];
d[i]=a[i][i]-s;
}
y[0]=b[0];
for(i=1;i<n;i++)
{
s=0;
for(k=0;k<i;k++)
s=s+L[i][k]*y[k];
y[i]=b[i]-s;
}
x[n-1]=y[n-1]/d[n-1];
for(i=n-2;i>=0;i--)
{
s=0;
资源评论
