소개글
한양대학교 이관수 교수님 수치해석 연습문제 숙제입니다
다른 교수님들에 비해 빡쌔게 한 숙제입니다
목차
● 문 제 5-(c) Gauss 소거법
● 문 제 7-(d) Crout 법
● 문 제 8 TDMA 법
● 문 제 10 Gauss-Seidel 법
● 문 제 11-(d) SOR 법
● 문 제 20, 21 구조물 문제
●문 제 25 관회로망 문제
# 첨 부 : 함수 프로그램 (Gauss소거법, 피봇팅, 전진대입법, 후진대입법, Gauss-Seidel법, LU분해법)
본문내용
<Coding>
#include<stdio.h>
#include<math.h>
#define n 4
void elim(double A[n+1][n+1],double b[n]); //가우스소거법
void pivot (double A[n+1][n+1],double b[n+1], int k); //피봇팅
void back(double U[n+1][n+1],double b[n+1],double x[n+1]); //후진대입법
void main()
{
double A[n+1][n+1]={ //계수입력
{0,0,0,0,0},
{0,-1,1,2,3},
{0,3,-1,1,-1},
{0,-1,3,4,1},
{0,1,2,1,2},
};
double b[n+1]={0,11,3,7,7};
double x[n+1];
int i;
elim(A,b); //가우스소거법
back(A,b,x); //후진대입법
for(i=1;i<=n;i++)printf("x[%d]=%.5f\n",i,x[i]); //결과출력
}
<Result>
x_1 : 1.000000 x_2 : -1.000000
x_3 : 2.000000 x_4 : 3.000000
우리는 확산 마무리 오차를 줄이기 위한 부분 피봇팅을 통하여 소수점 여섯째 자리까지 오차없이 손으로 푼 것과 같은 해를 구할 수 있었다.
참고 자료
없음