공학수치해석

소개글

공학수치해석에 대한 보고서 자료입니다.

목차

1. Problem

2. Idea

3. Programing with C

4. Result

5. Discussion

6. Reference

본문내용

.9e \n", i+1, p, fp);
// 오차한계보다 작을 때까지 지속
if ((fabs(p - p1) < tolerance)) {
printf("\nApproximate solution P = %15.9f \n", p);
printf("with Func(P) = %15.9f\n", fp);
printf("Number of iterations = %3d\n", i+1);
printf("Tolerance = %15.9e\n", tolerance);
ok = false;
} else {
i++;
p0 = p1; f0 = f1; // 새로운 p0, f0, p1, f1 값
p1 = p; f1 = fp;
}
}
if (ok) {
printf("\nIteration number %3d", num);
printf(" gave approximation %15.9f\n", p);
printf("Func(P) = %15.8f not within tolerance : %15.9e\n", fp, tolerance);
}
}
// 구하려는 해의 함수
double func(double x)
{
return (-1/(sqrt(x))+(1/k)*(log(RE*sqrt(x)))+(14-5.6/k));
}
4. Result
a. double p0 = 0.0001; // 임의값1
double p1 = 0.0002; // 임의값2
b. double p0 = 0.000001; // 임의값1
double p1 = 0.000002; // 임의값2
5. Discussion
그래프에서 알 수 있듯이 새로 정의한 함수는 x값이 양수일 때만 정의된다.
따라서 초기 임의값 두 개를 잘 못정하게 되면 아래 그래프처럼 새로운 p값에서 함수가 정의되지 않아 더 이상 진행하지 못하고 오류가 뜨게 된다.
따라서 초기 값을 지정할 때에는 0.0001 0.0002 와 같이 매우 작은 숫자로 하였다.
그 결과 0.005121894 라는 해를 얻을 수 있었다.
6. Reference
http://blog.naver.com/lovedh0920?Redirect=Log&logNo=150003022763
http://www.zappygo.com
http://yscec.yonsei.ac.kr/-Basics of CFD coding.pptx
Gerald and Wheatley, Applied Numerical Analysis, 7e, Pearson Education, Inc.

키워드

공학수치해석 ,   프로그램,   MATLAB