関数内で値の更新がうまくできていないと思われます。

Question

前提・実現したいこと

二分法というものを使って x - cos(x) = 0 の近似解を求めようとしています。
作ったプログラミングでエラーは発生していないのですが、出力される値が想定しているものではありません。

二分法の原理はとある値a, b と c = (a + b) / 2 を用いてf(a)とf(c)を計算し、f(a)とf(c)が異符号だったらaとcの間に解があるので新たにb = cとし、同符号だったらaとcの間に解はないので新たにa = cとする、というものを繰り返して設定した誤差の値より小さくなったら終わりという形になります。

発生している問題・エラーメッセージ

プログラムを実行すると実行結果が2.50をずっと繰り返しているだけになります。
誤差の設定がうまくいっていないのか、a,b,cそれぞれの値をうまく更新できないないのかのどちらかだと思われます。

該当のソースコード

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

/*与えられた式の右辺を計算する関数*/
double func(double hensu)
{
  double migi;
 migi = hensu - cos(hensu);
  
  return migi;
}


/*二分法を計算する関数*/
double bisec(double a, double b, double c)
{
  double kekka, gosa;

  c = (a + b) / 2;
  kekka = func(a) * func(c);


  if (kekka < 0){
    b = c;
  } else {
    a = c;
  }
  
  return kekka;
}

/*初期値はa = 0, b = 5 , gosa = 1.0*10^-8 */
int main(void)
{
  double fx, kekka, seido, hensu, atai, gosa;
  
  seido = pow(1.0, -8);
  
  do {
    double a = 0.0, b = 5.0, c = (a + b) / 2;

    bisec(a, b, c);

    if (a > b){
      gosa = a - b;
    } else {
      gosa = b - a;
    }

    printf("%lf\n", c);

  } while (gosa > seido);
  
  return 0;
  
}

補足情報（FW/ツールのバージョンなど）

https://replit.com/languages/c
このサイトでコードを書いて実行しています。

Comments

[m.ts10806]

質問は編集できます。間違えたのであれば新しく投稿するのではなく質問を修正してください。
回答つく前に削除依頼を。
またコードはマークダウンのcode機能にてご提示ください

[fana]

> x - cos(x)の近似解

という言葉で一般的に通用するのでしょうか？
解こうとしている方程式は，
x - cos(x) = 0
ですか？

Answer 1

関数 bisec は何一つ有益な仕事をしていません．

※また，仮に関数 bisec が何かしらの「更新」をしてくれたとしても，
main関数側のdo～while内で毎回 a や b の値を初期化しているので，やはりうまいこと動かないでしょう．

問題に取り組むよりもまえに，まずはC言語の基礎的な知識を習得する必要があるでしょう．

---

やさしさに満ち溢れたヒント：

とりあえずこのあたりの事柄から復習なりすれば良いのではないかな，と．

void f( int a )
{  a = 800;  }  //←この処理は，main関数側の変数の値に全く影響を与えない．

int main()
{
  int a = 1;  //aの値は1である
  f( a );  //関数fにaの値を引数として渡す
  printf( "%d", a );  //←これは 1 と表示される．800 ではない．
}

Answer 2

まず、このプログラムがどう動くのかを考えてみてください。
「実行するとまず〇〇行が動いて次に△△行が動いてこの変数がXXになって次に・・・」
という感じです。
そして、各行の前か後に printf で "A""B""C" とか "あ""い""う" とか(他の printf と区別するための文字)と共に 変化するはずの変数の値 を表示する行を追加して実行してみて、
想定した通りに動作しているかを確認(「A通ってB通って変数がこうなってCのif文で…あれ？」と行った感じ)しては如何でしょう。

関数で値の更新がうまく出来ていないと思われるのなら、関数を呼ぶ前後で変数の値を表示してみたり、関数内の最初と最後で表示してみる等で確認できると思います。

全体を動かしてから考えるよりも、少しずつ動かして確認していくほうが良い場合が多いかと思います。