수치해석제2장

목차

없음

본문내용

연습 문제 2-4

9 . 다음 4차 다항식

는 구간 [-1, 0] 과 [0, 1]에 각각 하나씩의 실근을 가지고 있다. 이 근들은
이내의 정확도로 다음 각 방법을 이용하여 구하라.

(a)가위치법
프로그램
#include <stdio.h>
#include <math.h>

void main()
{
double fa[20],fb[20],fp[20],x;
double p[20],a[20],b[20];
double n[20],m[20];
int i,j;
int max=20;

for (j=1;j<=2;j++) {
n[1]=0; n[2]=0;
m[1]=-1; m[2]=1;
a[1]=n[j];
b[1]=m[j];
p[1]=0;
for ( i=1 ; i<=max ; i++) {
fa[i]= 230*pow(a[i],4) + 18*pow(a[i],3) + 9*pow(a[i],2)
- 221*a[i] - 9;
fb[i]= 230*pow(b[i],4) + 18*pow(b[i],3) + 9*pow(b[i],2)
- 221*b[i] - 9;
p[i+1]= b[i] -(b[i]-a[i]) * fb[i] / (fb[i]-fa[i]);
fp[i+1]= 230*pow(p[i+1],4) + 18*pow(p[i+1],3) +
9*pow(p[i+1],2) - 221*p[i+1] - 9;
printf(" P%d=%.10f
",i+1,p[i+1]);
if (fa[i]*fp[i+1]==0) goto label_1;
else if
(fabs(p[i+1]-p[i])<pow(10,-6)) goto label_1;
else if
(fa[i]*fp[i+1]<0) { a[i+1]=a[i];
b[i+1]=p[i+1]; }
else
a[i+1]=p[i+1];
b[i+1]=b[i];
}
label_1:
printf("
");
}
}

참고 자료

없음

압축파일 내 파일목록

num1-2.c
num2-11a.c
num2-11b.c
num2-11d.c
num2-662.c
num2-9-1.c
num2-9-2.c
num2-9-3.c
taylor1.c
num1-2.txt
...