二分法で関数の最小値を求める

質問

題名通りのプログラムを組んだのですが、どうしてもうまくいきません。どこが間違っているのでしょうか？教えてください。

以下プログラム

#include< stdio.h >
#include< stdlib.h >

double minfx( double (*func)(double x), double a, double b )
{
if( b-a < 0.0001 )
{
if(func(a)>func(b))
{
return func(b);
}
else
{
return func(a);
}
}
else
{
if(func(a)>func(b))
{
a = (a+b)/2;
}
else
{
b = (a+b)/2;
}
minfx(func,a,b);
}
}

double f1(double x)
{
return (x-1.0)*(x-1.0);
}

double f2(double x)
{
return 3.0*x*x*x*x-4.0*x*x*x-12.0*x*x+10.0;
}

void main()
{
double ans1,ans2;

ans1 = minfx( f1, 0.0, 5.0 );
ans2 = minfx( f2, 0.0, 5.0 );

printf("関数f1の最小値：%lf \n", ans1);
printf("関数f2の最小値：%lf \n", ans2);
system("PAUSE");
}

これを実行するとf1とf2の値が0.000000といった感じで値が代入できていないと返されます。流れ的には合っていると思うのですが…なお、プログラムの形はできるだけそのままで、関数minfx内のみ変更する形でお願いします。自分としてはそこ以外に原因が考え付かないので…。

回答No.4

麻野 なぎ（@AsanoNagi） · ベストアンサー

では、ヒントだけ。

double を返す関数を抜けるとき、最後に実行されるのは、 return （式） でなければなります。

これだけでわからなければ、これを念頭に置いて、再帰で定義されたフィボナッチ数列とか、階乗とかのソースを探して眺めてみると、何かひらめくかもしれません。

あと、No.2 でも指摘されていますが、「最小値を求める」ためののだとすると、ロジック自体も間違っています。
よくある、方程式を解く（関数が０を横切る点を探す）のなら、一定の条件下では両端を見れば良いですが（両端の符号が異なれば、とりあえず少なくとも一個はゼロ点があるのは保証されるので）、最小値を探すのであれば、両端だけを見ても探索はできません。

お礼

ありがとうございます。ヒントを参考に探してみたいと思います。

回答No.3

redfox63

再帰呼び出しのminfxの戻り値を現在は破棄しています
これを捨ててしまっては 思い通りの結果にならないでしょう

補足

回答ありがとうございます。再帰呼び出しの仕方がまずいんでしょうか。再帰呼び出し時に値が破棄されているということはそういうことですよね。うーん、minfxではダメなんでしょうか…。しかし再帰ですからこれでいいと思うのですが…。

回答No.2

himajin100000

http://ideone.com/dlxJZ

まぁ確かに、27行目を修正する必要が有るのはわかる。確かにわかる。
(多分コンパイラにとって文法上そこにあることがわかるのが29行目がはじめなんだろう)

#system関数は意図的に削った。あともう一つのエラーは質問者宛に残しておく
====================
##ただ、正直出題者が何を考えて出したのかがわからない。
##例えば添付したような画像のグラフだと、最小値がいきなり右側にあることになっちゃうでしょ？(左側の垂直な直線がx = a,左側の垂直な直線がx = b,として点線がその中心とする)

補足

回答ありがとうございます。必要なのはプログラムが通るかどうかなので関数の値はあまり考慮していません。

回答No.1

Tacosan

minfx は何を返すの?