ニュートン法で方程式をとく

[kohaA]

ニュートン方で非線形方程式を解くプログラムを作ったつもりなのですが
まちがえているようで答えがおかしくなってしまいます
x^n+x^(n-1)+…x^2+x-(1/2)=0の式で
nを10、xを3、epsを0.0000000001でやったら3ちょっとくらいになるはずなのですがなりません…

どこがおかしいのでしょうか?

コード:

#include<stdio.h>
#include<math.h>

double f(double x , int n){
  double y=-1/2;
  while(n>=1){
    y += pow(x,n);
    n--;
  }
  return y;
}

double df(double x , int n){
  double y=0;
  while(n>=1){
    y += n*pow(x,n-1);
    n--;
      }
  return y;
}

int main(void)
{
  int i=0,n;
  double x,eps;
  printf("\n  x^n+x^(n-1)+...+x^2+x-(1/2)=0について\n");
  printf("  nの値:");
  scanf("%d",&n);
  printf("  初期値x:");
  scanf("%lf",&x);
  printf("  精度eps:");
  scanf("%lf",&eps);
  
  while(fabs(f(x,n))>eps)
    {
      x = x-f(x,n)/df(x,n);
      i++;
      printf("\n%3d回目    x=%lf    f(x)=%lf\n",i,x,f(x,n));
    } 
}

[みけCAT]

x^n+x^(n-1)+…x^2+x-(1/2)=0の式で
nを10、xを3、epsを0.0000000001でやったら3ちょっとくらいになるはずなのですがなりません…

どこがおかしいのでしょうか?

• この予想がおかしいです。3のn乗(nは1以上10以下の整数)を全部足したら明らかに-0.5を超え、方程式は成立しません。
• 関数fの最初

  コード:

    double y=-1/2;

  は、-1を2で割った商の0をyに代入しています。

  コード:

    double y=-1.0/2.0;

  としてください。

[kohaA]

ありがとうございます
小数点をつけて書くことの重要さが分かっていませんでした・・・
3ちょいは記憶違いで0.3ちょいでした^^;
毎回ありがとうございます！！！