n원 1차 연립 방정식 해 구하기

수치해석 시간 때 레폿으로 했던거 같군요;; ㅠㅠ

n원 1차 연립방정식 해 구하는 프로그램 입니다.

가우스 소거법을 이용해 역진대입법으로 풀었습니다.

소스코드 입니다.

#include <stdio.h> #include <stdlib.h> double** Get_Memory_2(int b); double* Get_Memory_1(int b); int main(void) { int Unknown_Variable; //모르는 변수의 개수 double** add_matrix; //첨가행렬의 원소를 담을 포인터 double* answer; double Temp_Num; //행과 행 나눈 것을 곱해서 빼야 하는데 그때 그 수를 저장 해둘 변수 while(1) { system("cls"); printf("몇 원 1차 방정식을 입력하실 건가요??(원 입력) : "); scanf("%d",&Unknown_Variable); add_matrix = Get_Memory_2(Unknown_Variable); //계수와 상수를 넣을 메모리 할당 받음 answer = Get_Memory_1(Unknown_Variable); //답을 넣을 메모리 할당 printf("Ax1 + Bx2 + ‥ + Cx3 = d 라면 A, B, C, d를 입력해주세요~! \n"); printf("모르는 변수가 3개이고 식이 3개라면 입력 개수는 12개! \n\n"); printf("==== 입력 시작 ====\n"); for(int i=0; i<Unknown_Variable; i++) { for(int k=0; k<=Unknown_Variable; k++) { scanf("%le",&add_matrix[i][k]); } } /*가우스 소거법 하는 루프문*/ for(int k=0; k<Unknown_Variable-1; k++) //첫번째 행을 가지고 끝내면 두번째 행을 가지고...하는 것 { for(int i=k; i<Unknown_Variable-1; i++) //다음 식을 빼는 루프문 { Temp_Num = add_matrix[i+1][k]/add_matrix[k][k]; for(int j=0; j<=Unknown_Variable; j++) //식과 식을 빼는 루프문 { add_matrix[i+1][j] = add_matrix[i+1][j] - ((Temp_Num)*add_matrix[k][j]); } } } printf("\n==== 가우스 소거법 후 ====\n"); /*가우스 소거법을 한 후 출력*/ for(int i=0; i<Unknown_Variable; i++) { for(int j=0; j<=Unknown_Variable; j++) { printf("%.3f ",add_matrix[i][j]); } printf("\n"); } /*역진대입법 들어감 마지막 변수부터 구함*/ for(int k=0; k<Unknown_Variable; k++) { double Temp=0; for(int i=0; i<Unknown_Variable; i++) { if(i!=Unknown_Variable-k-1) Temp += add_matrix[Unknown_Variable-k-1][i]*answer[i]; } add_matrix[Unknown_Variable-k-1][Unknown_Variable] = add_matrix[Unknown_Variable-k-1][Unknown_Variable] - Temp; /*해가 무수히 많을 경우와 해가 없을 경우 체크*/ if(add_matrix[Unknown_Variable-k-1][Unknown_Variable] == 0 && add_matrix[Unknown_Variable-k-1][Unknown_Variable-k-1] == 0) { printf("\n==== 에러 ====\n"); printf("해가 무수히 많다 \n"); break; } else if(add_matrix[Unknown_Variable-k-1][Unknown_Variable-k-1] == 0 && add_matrix[Unknown_Variable-k-1][Unknown_Variable] != 0) { printf("\n==== 에러 ====\n"); printf("해가 없다 \n"); break; } answer[Unknown_Variable-k-1] = add_matrix[Unknown_Variable-k-1][Unknown_Variable] / add_matrix[Unknown_Variable-k-1][Unknown_Variable-k-1]; } /*각 변수 출력*/ printf("\n==== 각 변수 ====\n"); for(int i=0; i<Unknown_Variable; i++) printf("x%d = %.3f \n",i+1,answer[i]); char a; printf("다시 계산 하려면 Enter키를 눌러주세요 ^^ \n"); fflush(stdin); a = fgetc(stdin); } return 0; } double** Get_Memory_2(int b) { double** a = (double **)malloc(sizeof(double) * b); for(int i=0; i<b; i++) { a[i] = (double *)malloc(sizeof(double) * (b+1)); } return a; } double* Get_Memory_1(int b) { double* a = (double *)malloc(sizeof(double) * b); return a; }

실행결과 입니다.

풀어본 문제

A + B = 0

3B + C = 0

4A + 2B + 3C + D = 6

6B + 2C + 3D + E = -15

4A + 6C + 3E = -22

5개의 변수를 구하여라

ㅇ

빨간색으로 네모친 부분이 사용자가 입력할 부분입니다.

필요하신 분은 사용하세요 ^^;

역진대입법.exe