Cài đặt Giải hệ phương trình tuyến tính bằng phương pháp khử Gauss trên ngôn ngữ lập trình C
Code:
#include "math.h"
#include "conio.h"
#include "iostream.h"
/*Nhap ma tran he so*/
double**NhapMaTran(int n){
double**Temp;
*Temp = new double[n];
for(int i = 0; i<n; i++){
Temp[i] = new double[n+1];
for(int j = 0; j<=n; j++){
if(j!=n)
cout<<"a["<<i<<"]["<<j<<"] = ";
else
cout<<"b["<<i<<"] = ";
cin>>Temp[i][j];
}
}
return Temp;
}
/*Xuat ma tran*/
void XuatMaTran(double**A, int n) {
for(int i=0 ; i<n; i++){
cout<<endl;
for(int j=0 ; j<=n; j++)
cout<<A[i][j]<<"\t";
}
}
/*Xuat he phuong trinh tuyen tinh */
void XuatHPT(double**A, int n) {
for(int i=0 ; i<n; i++){
cout<<endl<<A[i][0]<<".x0";
for(int j=1 ; j<n; j++) {
if(A[i][j]<0)
cout<<"\t"<<A[i][j]<<".x"<<j;
else
cout<<"\t+"<<A[i][j]<<".x"<<j;
}
cout<<"\t = "<<A[i][n];
}
}
void Exit(){
asm {
mov ah,4ch
int 21h
}
}
/*Chuong Trinh Chinh*/
void main(){
clrscr();
int n;
cout<<"Nhap Cap Cua Ma Tran n = ";
cin>>n;
double**A = NhapMaTran(n);
//XuatMaTran(A,n);
XuatHPT(A,n);
A = NULL;
getch();
Exit();
}