수다닷컴

질문 제목 : newton rahpson 알고리즘에 관하여..질문 요약 :맞는지가...질문 내용 : 책에 뉴턴 랩슨 알고리즘이 있는데 ..

/*
𝑓(𝑥)= 3𝑥^3 −2𝑥−1=0 의 근을 𝑥0=2
*/#include stdio.h
#include math.h
#include stdlib.hint a[] ={1, 3, 0, -2, -1};
int b[5];
int c[5];double differential(int x0) // 미분 함수
{
int k;
int n = 4;c[0] = b[0];for(k=1; k=n;k++){
c[k] = b[k]+ x0 * c[k-1];
printf(%d\n,c[k-1]);
}
printf(========differential===========\n);
return c[3];
}double synthetic_func(double x0){ // 조립제법 함수

int k;
double n=5;

//double x0 = 2;

b[0] = a[0];for(k=1; k=n; k++){
b[k] = a[k]+ x0 * b[k-1]; // b= 1 5 10 20 38 75
printf(%d\n,b[k-1]);
c[k] = b[k];
}
printf(=======synthetic_func============\n);
return b[5];
}void func(){// newton rahpson
double x0 = 2; // x0는 2
double km = 10; // 반복 횟수 10double e1=0, e2=0;
double f0, f1;

double x1, k;f0 = synthetic_func(x0);for(k=1; k=km; k++){
f1 = differential(x0); // 미분 함수 호출
x1 = x0-f0/f1;
//f0 = synthetic_func(x1);

if(fabs(x1-x0) e1)
{
printf(%lf\n, x0);
break;
}else if(fabs(f0) e2)
{
printf(%lf, %lf \n, x1,f0);
break;
}
else
x0=x1;printf(%.10lf\n, x1);
}

}int main()
{

func();
//int b = synthetic_func(2);return 0;
}
(순서도에 따라 출력한 결과.. 틀린것 같은 느낌이...)

이렇게 짜 봤는데.. 맞는 지 모르겠습니다.. ㅋ 오차범위를 생각하지말고 짰는데..

알고리즘에 파란색 부분이 표시되 있는데 책 밑에 순서도에는 저 부분이 없어서 어떤걸 해야맞는건지 모르겠는데

일단 순서도에 따라 저 부분을 없애면 숫자가 출력이 되고 저 부분을 넣으면 쓰레기 값이 나옵니다.. -0.idx00 머이런식으루..

맞는것 일까요..?? 지금 풀어 보는 중인데 .. 혹시나 다른 답이 있을까 해서... 올려봅니다.. 글구 문법같은거

이부분은 어떻하면 더 유용하거나 그런 스킬들 있으시면.. 알려주세요.. 그럼 감사~하겠습니다 ㅠ ~

신청하기