2个回答
展开全部
以下是二分法求已知区间的根,自己将其修改一下,增加一个区间确定就可以了
#include <math.h>
#include <stdio.h>
double f(double a, double b, double c, double x) // 求函数值
{ return a * x * x + b * x + c; }
int main()
{
double x, u, v, f1, f2, fm, a, b, c;
printf("请输入一元二次方程的系数:");
scanf("%lf%lf%lf", &a, &b, &c);
/* 这里可以增加是否合法二次方程的判断*/
printf("请输入区间端点:");
scanf("%lf%lf", &u, &v);
f1 = f(a, b, c, u);
f2 = f(a, b, c, v);
if (f1 * f2 > 0)
{
printf("区间内无解!\n");
return 1;
}
do
{
x = (u + v) / 2;
fm = f(a, b, c, x);
if (fm == 0)
break;
if (f1 * fm > 0.0) /* 同号时修改端点*/
{
u = x;
f1 = f(a, b, c, u);
}
else
v = x;
} while (fabs(fm) >= 1e-2);/*这里还可以在条件中增加循环次数控制*/
printf("根为:%8.3lf\n", x);
return 0;
}
#include <math.h>
#include <stdio.h>
double f(double a, double b, double c, double x) // 求函数值
{ return a * x * x + b * x + c; }
int main()
{
double x, u, v, f1, f2, fm, a, b, c;
printf("请输入一元二次方程的系数:");
scanf("%lf%lf%lf", &a, &b, &c);
/* 这里可以增加是否合法二次方程的判断*/
printf("请输入区间端点:");
scanf("%lf%lf", &u, &v);
f1 = f(a, b, c, u);
f2 = f(a, b, c, v);
if (f1 * f2 > 0)
{
printf("区间内无解!\n");
return 1;
}
do
{
x = (u + v) / 2;
fm = f(a, b, c, x);
if (fm == 0)
break;
if (f1 * fm > 0.0) /* 同号时修改端点*/
{
u = x;
f1 = f(a, b, c, u);
}
else
v = x;
} while (fabs(fm) >= 1e-2);/*这里还可以在条件中增加循环次数控制*/
printf("根为:%8.3lf\n", x);
return 0;
}
展开全部
double F(double a, double b, double c, double d, double x)//函数表达式
{
return (((a * x + b) * x) * x + d) / c;
}
double Function(double a, double b, double c, double d, double low, double high, double e)
{
double mid = (low + high) / 2;
if (F(a, b, c, d, mid) == 0)
return mid;
while ((high-low) >= e)
{
mid = (low + high) / 2;
if (F(a, b, c, d, mid) == 0)
return mid;
if (F(a, b, c, d, low)*F(a, b, c, d, mid) < 0)
high = mid;
else
low = mid;
}
return low;
}
你试试看,这是针对一般方程的。
{
return (((a * x + b) * x) * x + d) / c;
}
double Function(double a, double b, double c, double d, double low, double high, double e)
{
double mid = (low + high) / 2;
if (F(a, b, c, d, mid) == 0)
return mid;
while ((high-low) >= e)
{
mid = (low + high) / 2;
if (F(a, b, c, d, mid) == 0)
return mid;
if (F(a, b, c, d, low)*F(a, b, c, d, mid) < 0)
high = mid;
else
low = mid;
}
return low;
}
你试试看,这是针对一般方程的。
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询