C追赶法解三对角方程组
巡山小妖精
807次浏览
2020年07月29日 00:32
最佳经验
本文由作者推荐
埋怨读音-梦周
int main(void)
{
double c[100];
double b[100];
double a[100];
double l[100];
double u[100];
double f[100];
double x[100];
double y[100];
int n;
char ch;
printf("-------本程序为追赶法解三对角方程组计算程序-------
");
do
{
printf("请输入您要计算的三对角方程组是几阶方阵(默认最大100阶方阵):
");
scanf("%d",&n);
printf("请依次输入主对角线元素，中间以空格或回车隔开:
");
for(int i=0; i
scanf("%lf",&b[i]);
}
printf("请依次输入上次对角线元素，中间以空格或回车隔开:
");
for(i=0; i
scanf("%lf",&c[i]);
}
printf("请依次输入下次对角线元素，中间以空格或回车隔开:
");
for(i=0; i
scanf("%lf",&a[i]);
}t
printf("请依次输入方程组右边矩阵元素，中间以空格或回车隔开:
");
for(i=0; i
scanf("%lf",&f[i]);
}
for(i=0; i
if(i==0)
{
l[i] = b[i];ttt
}
else
{
u[i-1] = c[i-1]/(float)l[i-1];
l[i] = b[i] - a[i-1]*u[i-1];
}
}
for(int j=0; j
if(j==0)
{
y[j] = f[j]/(float)l[j];
}
else
{
y[j] = (f[j]-a[j-1]*y[j-1])/(float)l[j];
}
}
for(int m=n-1; m>=0; m--)
{
if(m==n-1)
{
x[m] = y[m];
}
else
{
x[m] = y[m]-u[m]*x[m+1];
}
}
for(j=0; j
printf("x%d的值为",j+1);
printf("%lf
",x[j]);
}
printf("是否继续下一个三对角方程组的计算？ 输入‘Y’继续
");
}
while((ch=getchar())=='Y'||(ch=getchar())=='y');
return 0;
}