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

#define max 1000                //最大繰り返し回数
#define eps 1.0e-5              //収束条件

double f(double x);
double df(double x);
void newton(void);

int main()
{
        newton();

        return 0;
}

/*ニュートン法の処理を行う関数*/
void newton(void)
{
        int count;
        double a,newa;

        count=0;

        printf("初期値を入力してください。\n");
        scanf("%lf",&a);

        for(;;) {
                count++;

                newa=a-f(a)/df(a);

                if(fabs(newa-a)<eps) break;

                a=newa;

                if(count==max) {
                        printf("収束しませんでした。\n");
                        exit(1);
                }
        }

        printf("解の値は %f\n収束するのに %d 回かかりました。\n",newa,count);

}

/*解く方程式定義する関数*/
double f(double x)
{
        return x*sin(x)+log(x);
}

/*解く方程式を一回微分したものを定義する関数*/
double df(double x)
{
        return sin(x)+x*cos(x)+1.0/x;
}