このプログラムがうまく作動しないのですが・・・
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
float f(float x);
float f(float x)
{
return (float)(x-sin(x)/cos(x));
}
void main()
{
float x1,x2,eps,f1,f2,xm,ff;
int i;
printf("Bisection method\n\n");
for(;scanf("%g%g%g",&x1,&x2,&eps)!=EOF;){
f1= f(x1); f2= f(x2);
if(f1>0){
xm= x1; x1= x2; x2=xm;
ff= f1; f1= f2; f2=ff;
}
printf("\nFinding a root between x1=%g and x2=%g\n", x1, x2);
printf("f(x1)=%g f(x2)=%g eps=%g\n", f1, f2, eps);
printf("\tx1\t\tx2\t\txm\t\tf(xm)\n");
if(f2<0){
printf("????\n");
continue;
}
i= 0;
for(;fabs(x1-x2)>=eps;){
xm=(x1+x2)*((float)rand()/(float)RAND_MAX)
ff= f(xm);
i++;
printf("%2d %15.6e %15.6e %15.6e %15.6e\n",i, x1, x2, xm, ff);
if(ff<0)
x1= xm;
else
x2= xm;
}
printf("A root found between %g and %g\n", x1, x2);
}
}
f(x)=x-tanx=0の解のうち３つを表示させるためにランダムサーチのプログラムを上記のように書いたつもりなんですがうまく作動しません。どこをどう直したらいいのでしょうか、教えてください。
