ニュートンの法則をつかって三次方程式を解こうと思ったのですが....
受付中
回答 1
投稿
- 評価
- クリップ 0
- VIEW 851
退会済みユーザー
前提・実現したいこと
Cを最近勉強しはじめた学生です
ニュートンの法則を使って三次方程式を解くコードを書こうと思い、コンパイルまではいけたのですが、自分の思った処理どおりに動きません涙
どなたか助けてください涙
該当のソースコード
include<stdio.h>
include<math.h>
define ER 0.001
double app( double x, double a, double b, double c, double d);
int main(){
int i;
double tmp;
double Sol[3];
double L=-10;
double J=10;
double a =1;
double b =-3;
double c =0;
double d =3;
tmp=L;
for(i=0 ; i<3 ; i++){
while( tmp <=J ){
Sol[i]=app(tmp,a,b,c,d);
while(fabs(Sol[i]-tmp) > ER){
tmp=Sol[i];
Sol[i]=app(tmp,a,b,c,d);
}
tmp=tmp+0.1;
if(i >0&&Sol[i] != Sol[i-1])break;
}
}
for(i=0 ; i<3 ; i++){
if(Sol[i] >=L && Sol[i] <=J)printf("%lf\n",Sol[i]);
}
return 0;
}
double app( double x, double a, double b, double c, double d){
return(x-( a*x*x*x + b*x*x + c*x + d)/(3*a*x*x + 2*b*x + c));
}
-
気になる質問をクリップする
クリップした質問は、後からいつでもマイページで確認できます。
またクリップした質問に回答があった際、通知やメールを受け取ることができます。
クリップを取り消します
-
良い質問の評価を上げる
以下のような質問は評価を上げましょう
- 質問内容が明確
- 自分も答えを知りたい
- 質問者以外のユーザにも役立つ
評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。
質問の評価を上げたことを取り消します
-
評価を下げられる数の上限に達しました
評価を下げることができません
- 1日5回まで評価を下げられます
- 1日に1ユーザに対して2回まで評価を下げられます
質問の評価を下げる
teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。
- プログラミングに関係のない質問
- やってほしいことだけを記載した丸投げの質問
- 問題・課題が含まれていない質問
- 意図的に内容が抹消された質問
- 過去に投稿した質問と同じ内容の質問
- 広告と受け取られるような投稿
評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。
質問の評価を下げたことを取り消します
この機能は開放されていません
評価を下げる条件を満たしてません
質問の評価を下げる機能の利用条件
この機能を利用するためには、以下の事項を行う必要があります。
- 質問回答など一定の行動
-
メールアドレスの認証
メールアドレスの認証
-
質問評価に関するヘルプページの閲覧
質問評価に関するヘルプページの閲覧
+3
面倒かとは思いますが、一気に作り上げようとせずに、
「この段階までの計算ではこういう値が出てくればOK」のような確認をしつつ、
段階的に作っていった方が良いと思います。
今の所やるべきは、どの段階でおかしい計算になっているのかを見つけることだと思います。
投稿
-
回答の評価を上げる
以下のような回答は評価を上げましょう
- 正しい回答
- わかりやすい回答
- ためになる回答
評価が高い回答ほどページの上位に表示されます。
-
回答の評価を下げる
下記のような回答は推奨されていません。
- 間違っている回答
- 質問の回答になっていない投稿
- スパムや攻撃的な表現を用いた投稿
評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。
15分調べてもわからないことは、teratailで質問しよう！
- ただいまの回答率 88.12%
- 質問をまとめることで、思考を整理して素早く解決
- テンプレート機能で、簡単に質問をまとめられる
質問への追記・修正、ベストアンサー選択の依頼
granfa_yuzo
2017/02/12 23:20
自分の思った処理どおりに動きません とは具体的にどのように違いがあるのでしょうか？
退会済みユーザー
2017/02/12 23:23
L~J の区間にある実数解を出したいのですが、出力後何も表示されないままなんです