수치해석 연습문제풀이 3장

소개글

한양대학교 이관수 교수님 수치해석 연습문제 숙제입니다

다른 교수님들에 비해 빡쌔게 한 숙제입니다

목차

● 문 제 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

우리는 확산 마무리 오차를 줄이기 위한 부분 피봇팅을 통하여 소수점 여섯째 자리까지 오차없이 손으로 푼 것과 같은 해를 구할 수 있었다.

참고 자료

없음