mm.cpp

来自「第一次发 不懂规矩 请大家谅解 数值分析中的进退法的VC++实现 供大家参考」· C++ 代码 · 共 35 行

#include "math.h"
#include "stdio.h"
main()
 {	int n=0; 
 float t,h,a,b,c,d,f1,f2;
 printf("t,h:"); scanf("%f,%f",&t,&h);
 f1=t*t-2*t+4; 
 if((t+h)*(t+h)-2*(t+h)+4<=f1)
 {
	 c=t+h;f2=c*c-2*c+4;
} 
 else 
 {
	 h=-h;
	 c=t+h;
	 f2=c*c-2*c+4;
} 
 while(f2<=f1) 
 {
	 h=2*h;
	 d=t;t=c;f1=f2;c=t+h;
	 f2=c*c-2*c+4;n=n+1;
} 
 if(c<=d) 
 {
	 a=c;b=d;
} else 
 {
	a=d;b=c;
}
 printf("[a,b]=[%.2f,%f]\n",a,b); 
 printf("n=%d\n",n); 
}