ニュートン法で解を１つ求めて、解を使って３次式を２次式に変えて、解の公式で解を求めてみます。 初期値によっては収束しません。。 /* 3次方程式 */ #include <stdio.h> #include <math.h> #define EPS (1.0e-10) double cubic(double a, double b, double c, double d, double x0) { ＿double fx, ffx, x=x0; ＿int i; ＿for (i=0; i<1000; i++) { ＿＿fx=a*x*x*x+b*x*x+c*x+d; /* f(x) */ ＿＿ffx=a*x*x+b*x+c; /* f'(x) */ ＿＿if (fabs(ffx)<EPS) { ＿＿＿printf("0 div error\n"); ＿＿＿return 0; ＿＿} ＿＿if (fabs(fx)<EPS) break; ＿＿x=x-fx/ffx; /* x←x-f(x)/f'(x) */ ＿} ＿printf("x=%f\n",x); ＿return x; } void quad(double a, double b, double c) { ＿double dis; ＿if (fabs(a)<EPS) { /* a==0 */ ＿＿if (fabs(b)<EPS) { /* b==0 */ ＿＿＿printf("解なし\n"); ＿＿} else { ＿＿＿printf("x=%f\n",-b/c); ＿＿} ＿} else { ＿＿dis=b*b-4.0*a*c; /* b^2-4ac */ ＿＿if (fabs(dis)<EPS) { /* ≒0 */ ＿＿＿printf("x=%f\n", -b/(2*a)); ＿＿} else if (dis>0) { /* 実数解 */ ＿＿＿printf("x=%f\n""x=%f\n", (-b+sqrt(dis))/(2*a), (-b-sqrt(dis))/(2*a)); ＿＿} else { ＿＿＿printf("x=%f%+fi\n""x=%f%+fi\n", -b/(2*a), sqrt(-dis)/(2*a), -b/(2*a), -sqrt(-dis)/(2*a)); ＿＿} ＿} } int main(void) { ＿double a,b,c,d; /* ax^3+bx^2+cx+d=0 */ ＿double e,f,g; /* ax^2+ex+f=0 */ ＿double x,x0; ＿printf("a b c d ? "); scanf("%lf%lf%lf%lf",&a,&b,&c,&d); ＿printf("x0 ? "); scanf("%lf",&x0); //printf("%fx^3%+fx^2%+fx%+f=0\n",a,b,c,d); ＿x=cubic(a,b,c,d, x0); /* 解を１つ求める */ ＿/* 解を使って次数を下げる */ ＿e=x*a+b; ＿f=x*e+c; ＿g=x*f+d; ＿if (fabs(g)<EPS) { //＿printf("%fx^2%+fx%+f=0\n",a,e,f); ＿＿quad(a,e,f); /* 二次方程式を解く */ ＿} ＿return 0; } 実行例 a b c d ? 1 1 1 1 x0 ? 0 x=-1.000000 x=-0.000000+1.000000i x=-0.000000-1.000000i