수치해석 secant method

소개글

수치해석 secant method c++로 풀어보는 문제

목차

1. OS : Window Xp

2. Source code

3. Output

4. Graph

5. Discussion

본문내용

2. Source code
#include <stdio.h>
#include <math.h>

double f(double x)
{
double a;
a = sqrt(16 - (x+1)*(x+1)) + 2; // 1사분면의 값으로 계산할 경우
//a = -sqrt(16 - (x+1)*(x+1)) + 2; // 4사분면의 값으로 계산할 경우
return a;
}
int main() {
int n;
double x1, x2, x3;

x1 = 0.5;
x2 = 3.0;

for(n = 1; n<10; n++)
{
x3 = x2 - ( f(x2)*(x1-x2) ) / (f(x1)-f(x2) );
printf(x1 : %f , x2 : %f, x3 : %f\n, x1, x2, x3);

x1 = x2;
x2 = x3;

if(x3 > 3)
{
printf(%d번째에서 발산하다, n);
return 0;
}
}
return 0;
}

5. Discussion
Xi-1 = 0.5일 때 1사분면에서 양수값 그리고 4사분면에서 음수값을 갖는다. 이 때 1사분면의 값을 택하여 계산하면 x축 값이 3을 넘게 되어 아무런 y값도 가질 수 없게 된다. 반면 4사분면의 값을 택하여 계산하면 위의 결과값과 같이 2.464102에 수렴하게 된다.

참고 자료

없음