|
()
{
double x1=0;
double x2=1;
double x3;
double x3o=0;
double f1=-1;
double f2=-0.12319;
while(1)
{
x3=x2-f2*(x2-x1)/(f2-f1);
printf("%lf\n",x3);
double f3=func(x3);
if (fabs((x3-x3o)/x3)*100 < 0.0005 )
break ;
if (fabs((x3-x3o)/x3)*100 > 0.0005)
{
x3o=x3;
x1=x2;
x2=x3;
f1=func(x1);
f2
|
|
메뉴펼치기
