뉴튼방정식 질문입니다.
글쓴이: volition13 / 작성시간: 목, 2005/10/06 - 9:58오후
책에있는 예제대로 한번 만들어 봤습니다.
xn=1.5
xexac=1.36523
do n=0,100
e=(xexac-xn)/xexac*100
fxn=xn**3.+4*xn**2. +10
fxn1p=3.*xn**2+8*xn
print *,n,xn,fxn,e
xn1=xn-fxn/fxn1p
xn=xn1
end do
end
e가 줄어들지를 않습니다. 무엇이 잘못됬는지요..
Forums:
x^3 + 4x^2 + 10 = 0 을 푼 것인가요?그렇다면 xexa
x^3 + 4x^2 + 10 = 0 을 푼 것인가요?
그렇다면 xexac이 -4.495쯤 되야 하는데요.
기본적인 알고리즘은 맞게 코딩하셨습니다. 다만 Newton-Rhapson법의 종료
조건은 보통 일정 회수가 아니라 e값이나 f(x)가 주어진 오차 범위 이내에 들 때로
하지요.
[quote="doldori"]x^3 + 4x^2 + 10 = 0 을 푼
책에 있는 예제를 그대로 한건데요 문제는 예제에서 이 방정식의
정확해(exact solution)는 1.36523이라고 주었거든요..
정확해를 -4.5정도쯤으로 하면 e(오차)가 0.1대로 줄어드는데
책에서 주어진 정확해대로 한다면 400퍼센트가 넘게 납니다.
책이 틀린걸까요?
그런 것 같군요.
그런 것 같군요.
제목이 좀 .. ^^;
뉴튼의 중력을 수치적인 방법으로 푸시려는 건줄 알았어요..
제목이 "Newton-Rhapson Method에 대한 질문"였어야 할 것같다는...
사족을 달면서... ㅎㅎ 후다닥==33
댓글 달기