[수치해석] PIVOTING 이용한 가우스 소거법

소개글

임의의 n*n 의 연립 방정식을 풀수 있는 전천후 코드 입니다.
일반 가우스 소거법이 아닌 pivoting와 scaling 알고리즘을 포함하고 있어 막강합니다.

컴파일 실행환경

// Gauss method

본문내용

#include<stdio.h>
#include<math.h>
#include "matrix.h"

#define SCALE_ON 1
#define PIVOT_ON 1

void error_code(int code);
void display_matrix(int n, double **A, double *B);
int pivot(int k,int n,double **A, double *B, double *DUMMY);
void scaling(int n,double **A,double *B,double *DUMMY);
FILE *fp;
void main()
{
double **A; // 계수 행열
double *B, *X, *DUMMY; // 벡터
double sum,number,m; // 각각 후위 대입과정의 임의 변수, 행열 요소값, multiplier
int i,j,k; // 행열 요소의 위치 첨자
int n,first=1; // 행열의 크기

fp=fopen("d:\gauss.txt","w");
printf("*** Gauss elimination method ***\n\n");
printf("*** Purpose to solve AX=B ***\n");
fprintf(fp,"*** Gauss elimination method ***\n\n");
fprintf(fp,"*** Purpose to solve AX=B ***\n");
// Make the matrix
while(n<=0)

참고 자료

없음