1. 程式人生 > 程式設計 >C++實現二分法求連續一元函式根

C++實現二分法求連續一元函式根

本文例項為大家分享了C++實現二分法求連續一元函式根的具體程式碼,供大家參考,具體內容如下

設計一個用二分法求連續一元函式根的通用函式solve
此函式有三個引數:

  • 第一個是函式指標,指向所要求根的連續函式
  • 第二、三個引數指出根的區間,且確保函式在區間的兩個端點異號

函式的返回值為求得的解

要求編寫main函式如下:

double fun(double x)
{
 double y;
 y=4*pow(x,3)-6*pow(x,2)+3*x-2;
 return y;
}

int main()
{
 cout<<"4*x^3-6*x^2+3*x-2=0在區間(1,2)的根為 x="<<solve(fun,1,2);
 return 0;
}

C++實現:

#include <iostream>
#include <cmath>

using namespace std;

double solve(double (*fun)(double x),double a,double b);

double fun(double x);

int main() {
 cout << "4*x^3-6*x^2+3*x-2=0在區間(1,2)的根為 x=" << solve(fun,2);
 return 0;
}

double solve(double (*fun)(double x),double b) {
 double i = b - a;
 double c = (a + b) / 2;
 while (i > 0.0000001) {
  i = b - a;
  if (fun(c) == 0)return c;
  if (fun(c) * fun(a) < 0) {
   b = c;
   c = (a + b) / 2;
  } else {
   a = c;
   c = (a + b) / 2;
  }
 }
 return c;
}

double fun(double x) {
 double y;
 y = 4 * pow(x,3) - 6 * pow(x,2) + 3 * x - 2;
 return y;
}

總結:

  • 函式與指標的結合
  • 注意返回的型別與要求

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支援我們。