|
초기치 벡터 x를 입력
printf(\"\\nx(0)를 입력해 주세요 : \\n\");
for(i=0;i<n;i++)
scanf(\"%lf\", &x[i]);
Jacobi(n, Max, esp);
}
void Jacobi(int n, int Max, double esp)
{
int i, j, k;
double sum, err, next_x[50];
for(k=1;k<=Max;k++){
printf(\"\\nk = %d \\n\", k);
for(i=0;i<n;i++){
sum=0.;
fo
|
|
|
Jacobi Method
#include <stdio.h>
#include <math.h>
#define x 1
#define y 2
#define del_x 0.05
#define del_y 0.05
#define NM 21
#define JM 41
#define err_max 0.01
void main()
{
int i,j;
int n=1;
double err;
double T[NM][JM];
double old_T[NM][JM];
double fabs(double);
//Matrix 초기화
for(
|
|
|
Jacobi법으로 해를 구하여보자.
--
이고 으로 가정
이므로
이다.
x
계산
비고
첫단계
두번째단계
세번째단계
네번째단계
x
계산
비고
다섯번째단계
여섯번째단계
일곱번째단계
여덟번째단계
x
계산
비고
아홉번째단계
열번째단계
열한번째단계
|
|
|
for(j=0;j<n;j++){
if(i!=j)
sum=sum+A[i][j]*x[j];
}
next_x[i]=(A[i][n]-sum)/A[i][i];
printf(\"%lf \\t\", next_x[i]);
}
// 2차_norm 으로 오차를 계산.
norm = 0.;
sum = 0.;
for(i=0;i<n;i++){
for(j=0;j<n;j++)
sum = sum + pow(A[i][j],2);
}
norm = pow(sum,0.5);
if(norm <= esp) // 허용오차보다
|
|
|
저지른 첫 번째 잘못인 중용을 지키지 못한 것과 두 번째의 잘못인 거짓말은 하나의 맥락에서 설명될 수 있을 것이다. 우선 거짓말을 한 현수는 습관화된 덕을 가지지 못하고 있다. 항상 정직한 성향을 가지고 있거나 그러한 습관을 가지고 있
|
|
메뉴펼치기
