二分法.cpp

来自「定步长辛普森公式计算.二分法.列主高斯消去法.龙贝格.牛顿迭代法.弦截法.逐次超」· C++ 代码 · 共 40 行

#include "math.h"
#include "stdio.h"

double Func(double x)
{
	return (sin(x)-(x/2)*(x/2));
}

void main()
{
	double a,b,l,eps;
	a=1.5;
	b=2;
	l=(b-a)*0.5;
	eps=1;
	for(int i=1;i<=7;i++)
		eps*=0.5;
	if(Func(a)<0)
	{
		double t;
		t=a;
		a=b;
		b=t;
	}
	double t;
	while(l>eps)
	{
		t=(a+b)/2;
		if(Func(t)>0)
			a=t;
		else
			b=t;
		l*=0.5;
	}
	printf("%.6f\n",t);
}