Newton-Raphson Method 사용
마디
질문 제목 : Newtddon-Raphson Method 사용한 소스인데 왜 안되는지 봐주세요.
질문 요약 :Newton-Raphson Method으로f(x)=exp(x)-1을 구하는 소스인데요.
초기값 설정에 따라출력되는값이 너무 차이가크다는점과
while(e0.0000001);이 부분 이해가 안되요.
질문 내용 :
#include stdio.h
#include math.h
int main(void)
{
double y, yd, x_1, e;
double x=2;//초기값설정 2로 해놨네요.
int i;
y = exp(x)-1;// 원래식
yd = exp(x);//미분한식
do //첫번째 한번은 돌아야 e값이 정의 되니까 do while사용.
{
x_1=x;
x=x-(y/yd); //x값은 점점 답에 가까워짐.
e=x_1-x;
}while(e0.0000001); //여기가 이해가 안되는데 e는 오차값인데 오차값이0.00000001보다 클때 반복하는게맞지 않나요? 그런데부등호는 이렇게 바꾸면 실행시켰을때 아무 작동도 안하더라구요.
printf(Solution by Newton-Rhapson method is %1.4lf\n, x); // 오차가 적은 근사치 x출력.
return 0;
}
제가 아는 부분은 주석 달았습니다.
단순히 생각해봤을 때 x값은 점점 답에 가까워고 그 오차값이 어느수준 이하로 떨어지면 멈추는데
이렇게 동작하는 경우 초기값이 달라도 같은 값이 출력되어야하지 않나요?
이 소스는 인터넷에서 찾은건데 초기값에 따라서 값 차이가 너무 크게 나네요..
그리고 두번째 질문은 while(e0.0000001); 이부분이에요. 주석에도 써놨는데
오차값이 0000001 이 값보다 클 때 동작해야된다고 생각하는데 부등호가 반대로 되어있네요.
그래서 부등호를 이쪽으로 바꿨더니 이경우는 아예 아무것도 안나오는..
exp(x)-1이 음수라 그런가 싶어서 exp(x)+1로도 해봤는데 역시나 안되더라구요..
질문 두가지 답변좀 해주세요.
-
민구
감사합니다^^