2장-연습문제(8번, 14번) [문제 8] 어떤 오리피스 유량계의 유량계수()는 다음의 실험식을 만족한다. ... 따라서 이 식은 해를 여러 개 가진 함수이다. 3장-연습문제(28번,2-4) [문제 28] 다음 구조물의 무게에 걸리는 각각의 힘과 반력을 구하라. ... ; f2=fcn(x2new); printf("X1 = %f, X2 = %f\nf1 = %f, f2 = %f\n\n",x1,x2new,f1,f2); x1=x2new; x2new=x1-
수치해석 2장 과제 (a) Use the Newton-Raphson method and find the smallest positive root correct to at least ... method의 반복 횟수가 Secant method의 반복 횟수보다 적은 것을 확인할 수 있지만, Newton-Raphson method는 주어진 함수의 미분형을 구해야 한다는 점에서 각각 장단점이 ... 이런 결과가 나오게 된 이유는 주어진 식 의 근이 2개이기 때문이다. 위 식의 실제 그래프를 살펴보면 이와 같은 모양으로 생겼고, 두 개의 근이 있다는 것을 알 수 있다.
(c) 그래프로 본 결과와 컴퓨터로 해석한 결과값이 같음을 알 수 있다. (e) 그래프로 그려본 근과 해석적으로 컴퓨터를 이용해서 구한 근이 같은 것을 알 수 있다. 11. ... 그래프로 봤을 때 근은 1.5 와 2 사이에 있는 것을 알수있다. 정확한 근을 구하기 위에 C코딩을 하고 결과값을 보면 다음과 같다. ... (b) 근이 1.2 근처에 있는 것을 알 수 있다. 따라서 으로 놓고 계산한다. 역시 muller법이 빠르다는 것을 알 수 있다.
수치해석 2장 ►문제 1번 h는 1.4 m 와 1.5 m 사이의 값이다. h에 대한 함수를 만들어야 증가탐색법을 쓸 수 있으므로 에서 를 우변으로 넘겨서 로 만들었다. ... x2=0.9, x3=0, old, error, f1, f2, f3; for(j=1;j ... 0.02, x3, old=0, error, f1, f2, f3; for(i=1;i
수치계산 2장(방정식의 근) 연습문제 풀이-1조 문제 8. 의 가장 작은 양의 근인 를 Newton법으로 구하려고 한다. ... -pow(y1,2)-1); } //함수 g를 선언 double g(double x2, double y2) { return(pow(x2,2)+pow(y2,2)-4); } //함수 f를 ... =1.001; x2=10; do //secant법으로 계산 { f1=function(x1); f2=function(x2); x3=x2-((x2-x1)*f2)/(f2-f1); if(x3error_range
*f; >> v2f=v2. ... {(4-3)(4-1)(4-2)} 2+ {(x-4)(x-3)(x-2)} over {(1-4)(1-3)(1-2)} 0+ {(x-1)(x-3)(x-4)} over {(2-1)(2-3)( ... (3)/(2)*2+(3.4-1)*(3.4-3)*(3.4-4)/1/-1/-2*5 ans = 5.5600 -> 1차 거의 비슷하다. 15. 2 다음과 같은 데이터가 주어진다. x 1 2
D1은 , D2는 으로 계산한 결과이고 D가 extrapolation 결과이다. 19.11 수치미분을 위해서 gradient 명령을 쓰면 다음과 같다. ... 이는 해석적인 결과와 반올림을 제외하면 동일하다. 18.9 a)해석적 해석적 적분결과는 다음과 같다. b)dblquad 명령어를 사용하면 다음과 같다. ... 즉 오차의 한계가 0이기 때문에 참값이 구해진 것이다. 17.6 a)해석적인 방법 즉 참값은 이다. b)합성 사다리꼴 공식(n=2) 합성 사다리꼴 공식, n=2의 경우에는 3개의 점이
이 구간을 다시 2등분하면 새로운 x3가 얻어진다. { x}_{3 }= {0.25+0.5 } over {2 }=0.375 이 값과 앞서 구한 값에 대한 백분율 상대오차는 { ε}_{ ... 이분법으로 이 방정식을 풀어보면, x1=0과 x2=0.5일 때 함수값은 각각 f(0)=1, f(0.5)=-0.3934이다. f(0)·f(0.5)0이므로 근은 x3=0.25와 x2=0.5사이에 ... 방정식의 근이 유효숫자 3자리까지 만족하기 위한 허용오차 한계는 { ε}_{s }= 0.5× { 10}^{(2-3) }=0.05% 이다. 따라서, εa