程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> 關於C語言 >> C語言 三次方程求法 代碼

C語言 三次方程求法 代碼

編輯:關於C語言
 

#include <iostream.h>
#include <math.h>
#include <iomanip>
#include <fstream.h>
#include <stdlib.h>

#define f(x) (1*x*x*x-3*x*x-6*x+8)

double bisection(double x,double y,double z)
{
double a=x,b=y,p;
while(((b-a)/2)>=z)
{ p=(a+b)/2;
if(f(p)*f(b)<=0)
a=p;
else b=p;
}
return p;

}

void main()
{
double a[5];
char filename[256];
cout.setf(ios::fixed);
cout<<"請輸入誤差容限:";
cin>>a[0];
cout<<"輸入目的文件全路徑:";
cin>>filename;
ofstream outfile(filename);
a[2]=(3-sqrt(21))/2;
a[3]=(3+sqrt(21))/2;
a[1]=a[2];
a[4]=a[3];
for(int i=1; ;i++)
{
if(f(a[1])>=0) a[1]--;

else if(f(a[4])<=0) a[4]++;
else break;
}
for(i=1;i<4;i++)

outfile<<"函數f(x)在區間("<<a[i]<<’,’<<a[i+1]<<")上的根為"<<bisection(a[i],a[i+1],a[0])<<endl;

cout<<"已經將結果寫入:"<<filename<<endl;
cout<<"誤差在"<<a[0]<<"之內."<<endl;
outfile.close();
}

input=============

f:output.txt //可以是任何路徑
1e-7

output============

函數f(x)在區間(-2.79129,-0.791288)上的根為-2
函數f(x)在區間(-0.791288,3.79129)上的根為1
函數f(x)在區間(3.79129,4.79129)上的根為4

 
  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved