質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.83%

  • C#

    6012questions

    C#はマルチパラダイムプログラミング言語の1つで、命令形・宣言型・関数型・ジェネリック型・コンポーネント指向・オブジェクティブ指向のプログラミング開発すべてに対応しています。

二分法 C言語のプログラム

解決済

回答 4

投稿

  • 評価
  • クリップ 1
  • VIEW 819

tennis1202

score 11

二分法のC言語のプログラムについて、エラーが出てうまくいきません。
どこが間違えているか、どこを直したらうまく出力できるか
教えていただけると幸いです。
以下がプログラムです。

include <stdio.h>

include <math.h>

int main(void){
double xL,xR,xM;
{
int count=0;
int xL=0,xR=2.0;
while(xR-xL > 0.0000001){
print("count++ xL,xR \n");
xM=(xL+xR)/2;
if (f(xL)*f(xM)<0)
xR=xM;
else{
xL=xM;
count ++;
}
    

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • rkhs

    2017/01/03 21:59

    どのようなエラーですか?

    キャンセル

  • seastar3

    2017/01/05 00:19

    これはC#ではなく、Cのコードです。

    キャンセル

回答 4

checkベストアンサー

+1

出力に関しては、seastar3さんのコメントをご参照ください。
その他気になった点を以下修正しています。
1.countについては、xM~xRの間に解がある時だけカウントするのではなく、iterationの数をカウントしたいのではないかと思います。この場合、count ++の場所が間違っています。(以下コードをご参照ください)
2.f(x)の定義がありません。2分法で求める対象の関数がなにか分かりませんが、例えばf(x)=(x-1)^3を対象の関数とした場合の例を、以下に記載します。
3.初期値が悪く、xL~xRの間に解が無い場合やtoleranceがシビアな場合、収束計算が終了しませんので、繰り返し計算の回数に上限を設けた方が良いかと思います。以下のコードでは50回を上限にしています。

#include <stdio.h>
#include <math.h>
double f(double x)
{
    return x*x*x-3*x*x+3*x-1;
}

int main(void){ 
double xM; 
xM=0.1;
int count=0; 
double xL=0,xR=1.9; 
while((xR-xL > 0.0000001) && (count < 50)){ 
printf("count++=%3d xL=%2.9f , xR= %2.9f, f1(xM)=%2.9f \n",count,xL,xR,xM,f(xM)); 
xM=(xL+xR)/2; 
if (f(xL)*f(xM)<0) 
{xR=xM;
}
else{ 
xL=xM; 
} 
count ++; 
}
}


以下上記の実行結果ですが、解析解のx=1付近で収束計算が終了しています。
イメージ説明

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

+1

そもそもコンパイルが成立すると思えないです。作製途中か何かだと思います。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

+1

f(0)*f(2)<0であり、かつfが[0,2]で単調であれば、提示したアルゴリズムは動くはずです。
エラーが出る原因は、プログラミング作法に問題があるか、使用しているfが仮定を満たしていないかのいずれかでしょう。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

0

これはC#ではなく、Cのコードですね。とりあえず動くように手直ししてみました。

#include <stdio.h>
int main(void){ 
 double xL,xR,xM; 
 int count=0; 
 xL=0.0,xR=2.0; 
 while(xR-xL > 0.0000001){ 
   printf("%4d  %10.9f  %10.9f\n",count,xL,xR); 
   xM=(xL+xR)/2.0; 
   if ((float)xL * (float)xM < 0)
     xR=xM; 
   else { 
     xL=xM; 
     count++; 
   }
 }
}


printfの使い方、xL,xRの二重宣言、(float)キャストの間違い、++直前の余分なスペース、}の不整合など、いかにも初心者のコードミスが散見されました。
私も最初は訳の分からない間違いを沢山して、直すのもお手あげで先輩に見てもらったものです。コードを一字一句忠実に読み解き、たくさんエラーを直して、デバッグのセンスを養っていきましょう。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2017/01/05 08:51

    自分でコーディングしておいてなんですが、(float)xL の浮動小数点型へのキャストは不要でしたね。
    元コードの f(xL)がキャスト操作かと推測して当てはめましたが、もしかしたら別のf()関数を定義するつもりなのかも知れません。
    その点の補足をお願いします。

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.83%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

  • 解決済

    2重ループのプログラム

    現在forの2重ループで下のような図形ができるように下のプログラムを書いたのですが、うまくいきません。なぜ下のプログラムではいけないのかは大体わかるのですが、どのように直せばうまく

  • 解決済

    1からNまでの二乗和を求める

    1からNまでの二乗和を求めるプログラムを作れ、という課題をやっているのですが、完成したと思いコンパイルし実行してみると、=0 という答えが返ってきました。 コンパイルした時、警

  • 受付中

    c言語のコードに一部誤りがあるようですがそれがどこなのかわかりません

    ある自然数nを素数の足し算に分割し、その組み合わせの数を求めるプログラム(例えば5なら5、2+3。8なら、2+2+2+2、2+3+3、3+5といったような分け方。n=0なら1が結果

  • 解決済

    gifの繰り返し回数の指定がわからない

    gif_libの画像処理に関する質問です。 https://github.com/katahiromz/imaio_bmp/blob/master/imaio_win.c この

  • 解決済

    while文のなかにif文

    その名の通りwhile文のなかにif文が入っているコードがあったのですが なんだか違和感を感じました。while文の中にif文をいれることってけっこうありますか? 詳しいかた教えて

  • 解決済

    コードの問題点 c

    このコードに疑問があります。 prime[500]と500個の配列をつくって、 prime[ptr++]=2としたら prime[1]=2となってしまいませんか? それとも[ptr

  • 解決済

    これは何故エラーになるのでしょうか?

    int *h = {1, 2, 3, 4, 5, 6};     printf("%d", *(h));          int h_2[] = {1, 2, 3, 4, 5, 

  • 解決済

    google spreadsheet Importxmlについて質問です。

    <?xml version="1.0" encoding="UTF-8"?> <countset xmlns="http://www.sitemaps.org/schemas/si

同じタグがついた質問を見る

  • C#

    6012questions

    C#はマルチパラダイムプログラミング言語の1つで、命令形・宣言型・関数型・ジェネリック型・コンポーネント指向・オブジェクティブ指向のプログラミング開発すべてに対応しています。