前提・実現したいこと
二分法を用いて3次方程式ax³+b²x+cx+d=0で実数解を求めるプログラムを作成したいです
(実数解が複数ある場合は考慮せず、指定範囲で避けるようにしています)
発生している問題・エラーメッセージ
下記のようなコードを作ったのですが正しい解が出力されない状況です。
該当のソースコード
c
1#include <stdio.h> 2#include <math.h> 3double func(double x,double a,double b,double c,double d); 4 5int main(void) 6{ 7 double a,b,c,d,x,i,j,k,Fi,Fj,Fk; 8 9 printf("abcd\n"); 10 scanf("%lf %lf %lf %lf",&a,&b,&c,&d); 11 12 printf("区間\n"); 13 scanf("%lf,%lf",&i,&j); 14 15 Fi=func(i,a,b,c,d); 16 Fj=func(i,a,b,c,d); 17 18 k=(i+j)/2; 19 20 Fk=func(i,a,b,c,d); 21 22 if(Fi*Fk!=0){ 23 while(j-i>pow(10,-6)) 24 { 25 Fi=func(i,a,b,c,d); 26 27 k=(i+j)/2; 28 29 Fk=func(i,a,b,c,d); 30 31 if(Fk!=0) 32 { 33 if(Fi*Fk<0) 34 { 35 j=k; 36 }else 37 { 38 i=k; 39 } 40 41 } 42 43 } 44 45 printf("解は%f\n",k); 46 47 }else 48 { 49 printf("解なし\n"); 50 } 51 return 0; 52 } 53 54 double func(double x,double a,double b,double c,double d) 55 { 56 return a*pow(x,3)+b*pow(x,2)+c*x+d; 57 } 58
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。