가우스사이델 반복법 VBA
- 최초 등록일
- 2014.10.23
- 최종 저작일
- 2014.03
- 4페이지/ 한컴오피스
- 가격 1,000원
목차
없음
본문내용
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 3
void jacobi_iteration(int max_iter);
void gauss_seidel_iteration(int max_iter);
void main()
{
int max_iter=6;
jacobi_iteration(max_iter);
gauss_seidel_iteration(max_iter);
}
void gauss_seidel_iteration(int max_iter)
{
double A[N][N],b[N],x_old[N],x_new[N],x_exact[N];
int i,j,iter;
double temp, err;
A[0][0]=4.0; A[0][1]=3.0; A[0][2]=0.0;
A[1][0]=3.0; A[1][1]=4.0; A[1][2]=-1.0;
A[2][0]=0.0; A[2][1]=-1.0; A[2][2]=4.0;
b[0]=24.0; b[1]=30.0; b[2]=-24.0;
<중 략>
err=sqrt(err);
printf("iter=%2d\t x[0]=%.4f\t x[1]=%.4f\t x[2]=%.4f\t error=%.4e\n",iter,x_new[0],x_new[1],
x_new[2],err);
iter=iter+1;
for(i=0; i<N; i++)
x_old[i]=x_new[i];
}
printf("근사값은 \n");
for(i=0; i<N; i++)
printf("x[%d]=%.4f\t",i,x_new[i]);
printf("\n");
}
참고 자료
없음