标题:三角系数矩阵方程组的c程序
只看楼主
zhangzujin
Rank: 1
等 级:新手上路
帖 子:276
专家分:0
注 册:2005-5-9
 问题点数:0 回复次数:0 
三角系数矩阵方程组的c程序

#include<stdio.h> #define N 3 int a[N][N],b[N],choice; float x[N];

int qd() { int i,flag=1; for(i=0;i<N;i++) if(a[i][i]==0) { flag=0; break; } return flag; }

void main( ) { int i,j; float sum; printf("make a choice:ub(1) or lb(2):"); scanf("%d",&choice);

if(choice==1) { printf("input A:\n"); for(i=0;i<N;i++) for(j=i;j<N;j++) scanf("%d",&a[i][j]);

printf("input b:\n"); for(i=0;i<N;i++) scanf("%d",&b[i]);

if(qd( )) { x[N-1]=b[N-1]/a[N-1][N-1]; for(i=N-2;i>=0;i--) { sum=0; for(j=i+1;j<N;j++) sum+=a[i][j]*x[j]; x[i]=(b[i]-sum)/a[i][i]; } } else { printf("Error,the main stream must be positive or negative.\n"); return; } } if(choice==2) { printf("input A:\n"); for(i=0;i<N;i++) for(j=0;j<=i;j++) scanf("%d",&a[i][j]);

printf("input b:\n"); for(i=0;i<N;i++) scanf("%d",&b[i]);

if(qd( )) { x[0]=b[0]/a[0][0]; for(i=1;i<N;i++) { sum=0; for(j=0;j<i;j++) sum+=a[i][j]*x[j]; x[i]=(b[i]-sum)/a[i][i]; } } else { printf("Error,the main stream must be positive or negative.\n"); return; } }

printf("throuh calculating,the result is:\n"); for(i=0;i<N;i++) printf("x[%d]=%.2f\n",i+1,x[i]); }

搜索更多相关主题的帖子: 方程 系数 三角 矩阵 
2005-05-11 23:37



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-17932-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.278383 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved