求解三次方程

原创 2007年09月17日 18:22:00
采用方法:先用二分法求一实根(三次方程必有一实根),然后化为二次方程,用解析法求得另外两根(包括实根和虚根)
#include <stdio.h>
#include <math.h>
#include <conio.h>
#define eps 1.0e-8
//三次函数
inline double pow3(double a3,double a2,double a1,double a0,double x)
{
 return ((a3*x+a2)*x+a1)*x+a0;
}
int pow3root(double a3,double a2,double a1,double a0,double r[])
{
 double left=0.0,right=0.0,half,temp,b,c;
 //不是三次方程
 if(a3<eps&&a3>-eps)
 {
  return -1;
 }
 //确定有实根的区间
 do
 {
  left-=1.0; 
  right+=1.0;
  temp=pow3(a3,a2,a1,a0,left);
  if(temp<eps&&temp>-eps)
  {
   r[0]=left;
   goto STEP2;
  }
  temp*=pow3(a3,a2,a1,a0,right);
  if(temp<eps&&temp>-eps)
  {
   r[0]=right;
   goto STEP2;
  }
 }
 while(temp>0);
 while(right-left>eps)
 {
  half=(right+left)/2.0;
  temp=pow3(a3,a2,a1,a0,half);
  if(temp<eps&&temp>-eps)
  {
   r[0]=half;
   break;
  }
  else if(pow3(a3,a2,a1,a0,left)*temp<0)
  {
   right=half;
  }
  else
  {
   left=half;
  }
 }
 r[0]=(right+left)/2.0;
STEP2:
 b=a2/a3+r[0];
 if(r[0]<eps&&r[0]>-eps)
 {
  c=a1/a3;
 }
 else
 {
  c=-a0/a3/r[0];
 }
 if(b<eps&&b>-eps)
 {
  if(c<=0)
  {
   r[1]=-sqrt(-c);
   r[2]=r[1];
   return 0;
  }
  else
  {
   r[1]=0.0;
   r[2]=sqrt(c);
   return 1;
  }
 }
 temp=b*b-4*c;
 if(temp<eps&&temp>-eps)
 {
  r[1]=r[2]=-b/2.0;
  return 0;
 }
 else if(temp>0)
 {
  temp=sqrt(temp);
  r[1]=(temp-b)/2.0;
  r[2]=-(temp+b)/2.0;
  return 0;
 }
 else
 {
  r[1]=-b/2.0;
  r[2]=sqrt(-temp)/2.0;
  return 1;
 }
}
int main()
{
 int flag=1,temp;
 double a3,a2,a1,a0,r[3];
 //pow3root(1.0,-3.0,3.0,-1.0,r);
 while(flag)
 {
  printf("请依次输入三次方程的四个系数,以空格键隔开。/n");
  printf("输入:");
  scanf("%lf %lf %lf %lf",&a3,&a2,&a1,&a0);
  switch(pow3root(a3,a2,a1,a0,r))
  {
  case -1:
   printf("这不是一个三次方程!/n");
   break;
  case  0:
   printf("结果:x1=%f  x2=%f  x3=%f/n/n",r[0],r[1],r[2]);
   break;
  case  1:
   printf("结果:x1=%f  x2=%f+%fi  x3=%f-%fi/n/n",r[0],r[1],r[2],r[1],r[2]);
   break;
  }
  printf("退出请按Q或q,继续按其他键。/n");
  temp=getch();
  if(temp==(int)('q')||temp==(int)('Q'))
  {
   flag=0;
  }
  printf("/n");
 }
 return 0;
}
 

【Openjudge:Noi】7891:一元三次方程求解 c++

【Openjudge:Noi】7891:一元三次方程求解 总时间限制:  1000ms  内存限制:  65536kB 描述 有形如:ax3+bx2+cx+d=0  这样的一...
  • u013675643
  • u013675643
  • 2016年07月01日 16:25
  • 1354

求解一元三次方程--c++

#include #include #include//调用了fabs、pow函数 using namespace std;   double f(int,int,int,int,doubl...
  • u012221917
  • u012221917
  • 2013年12月06日 23:22
  • 2116

蓝桥网 算法训练 一元三次方程求解

问题描述   有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根...
  • ShiWaiGaoRen12345
  • ShiWaiGaoRen12345
  • 2016年11月06日 12:27
  • 365

算法训练 一元三次方程求解 蓝桥杯

问题描述   有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间)...
  • sinat_35637319
  • sinat_35637319
  • 2017年02月25日 13:04
  • 739

算法训练 一元三次方程求解 二分

算法训练 一元三次方程求解   时间限制:1.0s   内存限制:256.0MB        问题描述   有形如:ax3+bx2+cx+d=0 这样的一个一元三...
  • jerans
  • jerans
  • 2017年04月26日 21:47
  • 588

一元三次方程求解c++实现

typedef  double Number; class CubicRealPolynomial { public:     static Number computeDiscriminan...
  • zg260
  • zg260
  • 2014年11月21日 15:36
  • 1518

EXCEL 单变量求解 解一元三次方程

目标单元格:是公式所在的单元格,也是方程所在的单元格,如图中的B1=3*A1^3-100*A1^2+300*A1-300 目标值:即目标单元格的值,想要让单元格的值变成多少?即要让方程的结果等于多少 ...
  • gufengaoyue
  • gufengaoyue
  • 2013年12月12日 18:24
  • 2759

一元三次方程的分治解法

/* 2011-09-21 f(x)=ax^3+bx^2+cx+d 分析:设三个根为 x0,x1,x2 有 f(x)=a(x-x0)(x-x1)(x-x2...
  • mobius_strip
  • mobius_strip
  • 2014年03月24日 02:25
  • 1829

求解一元三次方程的方法java实现

public double findFS(double a, double b, double c, double d) {   double FS = 0;   double A, B, C, ...
  • zhangxr123
  • zhangxr123
  • 2011年08月24日 17:30
  • 1510

牛顿迭代法在求解三次方程上的应用

牛顿迭代法是牛顿在17世纪提出的一种求解方程f(x)=0.多数方程不存在求根公式,从而求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。 设r是f(x)=0的根,选取x0作为r初始近似...
  • barry123
  • barry123
  • 2007年11月24日 21:26
  • 943
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:求解三次方程
举报原因:
原因补充:

(最多只允许输入30个字)