|
nitial guess를 0으로 주고 iteration 했을 때는 근을 찾지 못하고 종료하는 것을 보았다. 이로써 적절한 initial guess를 얻고 대입하는 것이 modified Secant Method를 사용하는데 있어서 또 하나의 과제라고 하겠다. 반대로 이 문제의 경우 friction factor의 방
|
|
|
Secant Method
뉴튼법을 수행하는데 생길수 있는 문제점은 도함수의 계산이다 다항식이나 다른 많은 함수에서는 문제되지 않는다고 하더라도 어떤 함수에서는 도함수를 계산하는 것이 매우 어려울 수 있다 그러한 경우 도함수는 후진 유한제차
|
|
|
Secant method at p0=-3 , p1=-2
for the value of x with x^3 + 3x^2 - 1 = 0
*/
#include<stdio.h>
#include<math.h>
#define TOL 0.0001
float f(float x);
float absol(float x);
main()
{
int n;
float p0, p1, p2;
n = 1;
p0 = -3;
p1 = -2;
p2 = 0;
while(absol(p1-p0)>=TOL)
{
p2 = p1 - f(p1)*(p1
|
|
|
method at %f\n",p2);
while(absol(p2-p1)>=TOL)
{
p1 = p2;
p2 = p1 - f(p1)/fp(p1);
n++;
printf("n= %d Pn= %f f(Pn)= %f \n",n,p2,f(p2));
}
printf("Solution is %f\n\n",p2);
return 0;
}
(2)Result
(3)Compare to 2.2-9 and 2.3-9
( Bisection method )
( Secant method )
I used the Bisection method for excer
|
|
|
Secant법을 추천하는 바입니다.
참고 문헌
Numerical Methods for Engineers(공학도를 위한 수치해석) Ⅰ. 수치해법이 필요한 이유
Ⅱ. 정확한 해를 구하는 과정
1. 구간법
a. 이분법
b. 선형보간법
2. 개방법
a. 고정점 반복법
b
|
|
메뉴펼치기
